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SECTION  I  1 

introduction'  ,  ■  *  1 

i 

i  !  I 

The  problem  of  specifying  the  in-flight  behavior  of  an  aircraft  can  be  considered 
aa  one  of  determining  the  effect  the  dynamic  motion  of  the  craft  has  on  the  fluid  through 
which  it  is  moving  and,  hence,  on  what  forces  and  moments  are  impressed  on  ttye 
vehicle  due  to  this  motion.  The  forces  and1  moments  thus  impressed  on  the  vehicle 
together  with  thrust  and  inertial  response  completely  determine  the  aircraft  behavior. 
The  determination  of  the  forces  and  moments  experienced  by  an  aircraft  is,  therefore, 
of  importance  for  determining  aircraft  performance  and  stability  and  contrpl  require¬ 
ments.  Methods  for  predicting  the  aerodynamic  behavior  of  conventional  aircraft  have 
been  established  and  are  documented  in  a  form  suitable  for  use  by  the  preliminary  1 
designer  in  Reference  1.  V/STOL  aircraft  of  the  lift  jet,  lift,  fan,  and  vectored 
thrust  type  pose  aerodynamic  problems  which  are  significantly  different  from  the  , 
aerodynamics  of  conventional  aircraft.  New  problem  areas  arise  mainly  through 
the  interference  of  the  lifting  jet  or  far.  efflux  on  the  relative  mainstream  flow  due  to 
aircraft  forward  speed  or  natural  wind.  Also,  due  to  its  particular  form  of  propulsion, 
a  V/STOL  vehicle  is  likely  to  fly  at  large  angle  of  attack  and/or  sideslip  during  some 
part  of  its  flight  path  so  that  nonlinear  aerodynamics  become  an  important  factor  in 
determining  the  aerodynamic  forces  and  moments  on  these  vehicles. 

i  : 

To  date  no  satisfactory  methods  for  estimating  the  power  induced  aerodynamics 
of  the  lift  jet,  lift  fan  or  vectored  thrust  V/STOL  aircraft  have  been  established.  The 
subject  of  nonlinear  aerodynamics  has  received  little  attention  and  no  generally  accepted 
prediction  methods  exist  at  this  time.  This  study  has  been  concerned  with  obtaining 
the  air  induced  forces  and  moments  on  V/STOL  aircraft  of  the  lift  jet,  lift  fan  and 
vectored  thrust  type  in  a  convenient  form  to  allow  motion  studies  to  be  made.  This 
requires  that  the  magnitude1  of  the  force  and  moment  components  be  determined  as 
functions  of  a  consistent  set  of  flight  variables.  1 

*  •  i 

These  force  and  moment  components  are  treated  as  being  functions  of  the  flight 
variables,  a  consistent  set  of  which  are  a  ,  0  ,  U ,  a  ,  0  ,  p  q  ,  r  and  any  higher 
order  of  terms  such  as  ^  which  may  be  considered  pertinent.  The  basic  problem 
of  any  aerodynamic  stability  study  thus  reduces  to  obtaining  accurate  functional 


1 


I 


) 


relationships  between  the  forces  and  moments  and  the  flight  condition,  for  example, 

L  .  .  .  .} 

These  functional  relationships  then  permit  all  the  necessary  aircraft  motions  to  be 
obtained  assuming  the  thrust  and  inertial  characteristics  of  the  vehicle  are  known. 

A  complete  functional  relationship  such  as  the  above  equation  is  seldom  obtain¬ 
able  in  practice.  However,  the  effect  of  the  most  significant  parameters  on  given 
components  or  parts  of  the  aircraft  is  usually  obtainable.  These  effects  are  added 
together  and  assumed  to  hold  for  the  complete  airplane. 

1.  PURPOSE 

The  purpose  of  this  investigation  was  to  develop  analytical  engineering  methods 
for  predicting  the  static  and  dynamic  stability  and  control  derivatives  and  force  and 
moment  coefficients  of  lift  jet,  lift  fan,  and  vectored  thrust  V/STOL  aircraft  in  the 

hover  and  transition  flight  regimes.  These  methods  were  to  take  into  account  the 
strong  power  effects,  large  variations  in  angle  of  attack  and  sideslip,  and  changes  in 
aircraft  geometry  that  are  associated  with  high  disk  loaded  V/STOL  aircraft  operating 
in  the  aforementioned  flight  regimes.  Where  appropriate,  use  was  to  be  made  of  high 
speed  computers  to  obtain  solutions  having  reasonable  time  periods  for  implementa¬ 
tion.  The  aircraft  configurations  studied  were  to  have  a  conventional  wing,  fuselage 
and  empennage.  The  methods  developed  were  to  be  suitable  for  use  by  design 
personnel  during  the  preliminary  design  and  evaluation  of  V/STOL  aircraft  of  the 
type  previously  mentioned. 

2.  TECHNICAL  APPROACH 

This  investigation  has  been  concerned  with  the  development  of  theories  for 
predicting  the  forces  and  moments  on  lift  jet,  lift  fan  and  vectored  thrust  V/STOL 
aircraft.  In  general,  due  to  the  complexity  of  the  interaction  between  the  propulsive 
system  and  the  airframe,  it  has  been  found  necessary  to  use  computer  programs  to 
enable  the  forces  and  moments  to  be  calculated.  Some  of  the  methods  developed  are 
simple  enough  to  be  reduced  into  a  handbook  procedure.  Examples  of  such  methods 
are  those  based  on  slender  body  theory  and  the  viscous  cross  flow  hypothesis. 

i 

This  investigation  has  followed  the  conventional  procedure  in  a  manner  most 
suitable  for  the  treatment  of  V/STOL  aircraft. 
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A  basic  assumption,  for  the  methods  developed  in  this  investigation,  is  that 
the  power  induced  effects  can  be  treated  separately  from  the  unpowered  aerodynamic 
effects  and  the  two  sets  of  solutions  added  together  to  produce  the  total  aerodyanmic 
effect. 

A  second  assumption  is  that  the  aerodynamics  can  be  treated  component  by 
component  and  the  results  added  together  for  the  total  aircraft.  This  assumption 
holds  for  both  the  power-induced  and  unpowered  effects.  Whenever  possible, 
mutual  interference  between  components  has  been  taken  into  account  to  obtain  a  more 
accurate  solution.  This  approach  is  essentially  the  procedure  used  in  Reference  1 
and  is  considered  to  be  most  suitable  for  use  as  a  design  tool. 

The  technical  approach  used  in  developing  the  prediction  methods  was  based  on 
theoretical  and  semi-empirical  analyses.  It  was  thought  that  correlation  of  wind 
tunnel  test  data,  in  terms  of  geometric  and  flow  variables,  would  be  an  impossible 
task  without  theoretical  and  semi-empirical  methods  to  indicate  correlation  vari¬ 
ables. 


An  extensive  literature  search  was  conducted  to  identify  available  test  data 
and  areas  for  which  further  data  was  required  to  develop  and  validate  the  prediction 
methods.  As  a  result  of  this  literature  search,  two  wind  tunnel  test  programs  were 
conducted. 

In  the  first  of  these  test  programs,  a  flat  plate  model  containing  up  to  three 
nozzles  was  tested  to  assist  in  the  development  and  validation  of  the  jet  flow  field 
theory.  The  results  of  these  tests  have  been  documented  in  Reference  3.  In  the 
second  test  program,  a  complete  configuration  model  was  tested  in  the  NASA 
Langley  Research  Center  V/STOL  tunnel.  This  model  consisted  of  two  wing-mounted 
vectored  thrust  engines  and  a  lift-jet  engine  embedded  in  the  fuselage.  The  vectored 
thrust  nozzles  were  designed  for  exit  positions  at  two  longitudinal  locations.  For  the 
aft  location,  two  nozzle  exit  diameters  were  available.  Data  from  this  test  has 
been  used  to  verify  the  prediction  methods.  A  data  report  for  this  test  will  be 
published  by  NASA  at  a  later  date. 

3.  REPORT  ORGANIZATION 

The  report  consists  of  four  volumes,  entitled: 

Volume  I  -  Theoretical  Development  of  Prediction  Methods 
Volume  n  -  Application  of  Prediction  Methods 


Volume  in  -  Manual  for  Computer  Programs 
Volume  IV  -  Literature  Survey 

In  Volume  I,  the  aerodynamic  prediction  methods  are  developed  in  a  form  suit¬ 
able  for  application  to  each  aircraft  component.  The  applicable  theoretical  analyses 
or  semi-empirical  bases  are  presented.  Empirical  coefficients  are  determined, 
where  necessary,  and  then  extensive  comparisons  of  calculations  with  test  data  are 
made. 

Volume  n  gives  detailed  examples  of  the  application  of  the  prediction  methods 
to  the  determination  of  the  aerodynamic  forces,  moments,  and  in  some  cases,  surface 
pressure  distributions,  on  the  aircraft  wing,  fuselage  and  empennage.  In  each  case, 
a  sample  problem  is  given  to  illustrate  the  application  of  the  methods  and  method 
limitations  are  also  discussed.  Also  discussed  in  Volume  n  is  the  wind  tunnel 
test  program  conducted  in  the  NASA  Langley  V/STOL  tunnel  and  the  use  of  the 
resulting  data  to  validate  the  applicable  prediction  methods. 

Volume  III  contains  a  detailed  description  of  the  computer  programs  developed 
in  this  investigation. 

Volume  IV  documents  the  results  of  the  literature  survey  conducted  to  identify 
existing  test  data  and  theoretical  methods  relating  to  this  investigation. 
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SECTION  n 


JET  FLOW  FIELD  THEORY 

A  fundamental  problem  in  the  development  of  methods  for  predicting  aerodynamic 
characteristics  of  lift-jet,  vectored  thrust  and  lift-fan  V/STOL  aircraft  is  that  of  for¬ 
mulating  a  mathematical  model  to  estimate  the  effects  of  the  propulsion  system  efflux 
interaction  with  a  crossflow.  During  the  transition  flight  phase,  this  efflux  is  directed 
at  larg*,  angles  to  the  freestream  and  has  a  significant  influence  on  the  aircraft  aero¬ 
dynamic  performance  as  well  as  on  its  stability  and  control  characteristics. 

A  number  of  analytical  formulations  of  the  problem  of  a  single  jet  exhausting  into 
a  crossflow  exist,  and  details  of  the  different  approaches  may  be  found  in  Reference  2. 
An  approach  to  the  problem  of  a  single,  normally  exhausting  jet,  which  appeared  to 
offer  possibilities  of  treating  more  complex  configurations,  is  given  in  Reference  4. 

In  this  reference,  an  entrainment  model  was  developed  from  dimensional  analysis  and 
physical  considerations.  The  force  on  the  jet  boundary  as  a  result  of  the  pressure 
differential  around  the  jet  was  accounted  for  by  a  crossflow  drag.  The  geometry  of  the 
Jet  cross  section  was  represented  by  an  ellipse.  Assuming  constant  and  equal  density 
for  the  jet  and  the  crossflow,  the  continuity  and  momentum  equations  were  solved  for 
the  jet  path.  The  jet-induced  velocity  field  was  then  determined  by  replacing  the  Jet 
by  a  distribution  of  sinks  and  doublets. 

The  analytical  model  described  above  has  been  further  extended  to  treat  jets 
exhausting  into  arbitrarily  directed  freestreams  as  well  as  multiple-jet  configurations. 
Multiple-jet  configurations  are  treated  as  combinations  of  discrete  jets,  with  leading 
jets  assumed  to  develop  independently  and  downstream  jets  assumed  to  exhaust  into  a 
freestream  of  reduced  dynamic  pressure. 

To  generate  data  which  would  substantiate  these  assumptions  and  also  provide 
information  for  further  refinements  to  the  analytical  model,  a  wind  tunnel  investigation 
of  Jets  exhausting  into  a  crossflow  was  conducted.  The  test  configuration  consisted  of 
a  four-foot  diameter  circular  plate,  containing  up  to  three  circular  Jets,  which  was 
elevated  12  inches  from  the  tunnel  floor  and  aligned  with  the  air  flow.  The  plate  con¬ 
tained  pressure  taps  to  determine  the  surface  static  pressures.  Jet  centerlines  and 
decay  characteristics  were  obtained  with  a  total  head  rake.  ■  Data  from  this  investiga¬ 
tion  have  been  presented  in  Reference  3. 
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Data  from  the  wind  tunnel  investigation  have  been  analyzed  and  used  to  verify 
some  of  the  assumptions  made  in  the  development  of  the  analytical  approach.  The 
mathematical  model  has  been  refined  using  data  from  the  wind  tunnel  investigation. 

The  single  Jet  is  considered  in  Sections  II.  1  and  II.  2.  The  details  of  the  analyti¬ 
cal  model  are  presented.  Refinements  of  the  empirical  parameters,  based  on  analysis 
of  the  experimental  dataware  discussed.  Extensive  comparisons  between  theory  and 
experimental  data  are  presented. 

Two-Jet  configurations  are  treated  in  Sections  n.  3  and  n.  4.  An  expression  for 
the  effective  dynamic  pressure  which  the  downstream  Jet  "sees"  as  a  result  of  the  block 
age  of  the  crossflow  by  the  upstream  Jet  is  derived.  The  test  data  are  used  to  verify 
that  the  upstream  Jet  develops  independently  of  the  downstream  Jet.  Calculations  of 
Jet  centerlines  and  induced  surface  static  pressures  are  compared  with  test  data  for 
a  number  of  two-jet  configurations. 

Applications  of  the  analytical  model  to  more  complex  configurations  are  pre¬ 
sented  in  Section  II.5.  Computations  of  the  induced  pressure  distribution  around  a 
two-jet  configuration,  with  the  jets  at  different  thrust  levels,  are  compared  with  test 
data.  The  extension  of  the  computational  procedure  to  the  treatment  of  more  complex 
jet  arrangements  is  carried  out  for  a  three-jet  configuration.  An  empirical  relation¬ 
ship  for  the  reduction  in  crossflow  dynamic  pressure  for  the  third  jet,  as  a  result  of 
blockage  by  the  two  upstream  Jets,  is  given.  Comparison  between  theory  and  data  for 
the  variation  of  surface  static  pressure  induced  by  a  three-jet  configuration  is  shown. 
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1.  SINGLE  JET  ANALYTICAL  MODEL 
a.  Normally  Exhausting  Jet 

Consider  a  circular  Jet  exhausting  at  a  right  angle  into  a  uniform  mainstream 
as  shown  in  Figure  1.  When  a  Jet  exhausts  at  an  angle  into  a  mainstream  it  is  deflect 
ed  partly  by  viscous  entrainment  and  partly  by  the  force  on  the  jet  boundary  resulting 
from  the  pressure  differential  around  the  Jet.  It  is  assumed  that  the  flow  is  incom¬ 
pressible  and  viscous  effects  other  than  entrainment  are  neglected.  Dimensional  con¬ 
siderations  then  suggest  that  the  entrainment  of  mainstream  fluid  per  unit  length  of 
jet  may  be  written  as 


£  - 


|0  £,  Uoo  d  cos  S  + 


p  E2  (Oj  ‘~Lfoo  &IQ&)  C 
/  +  £3  Uoo  cosG/cJj 


(i) 


A  fluid  particle  approaching  the  jet  near  the  plane  of  symmetry  will  be  more 
easily  entrained  by  the  jet  than  a  particle  that  is  moving  away  from  the  Jet.  The  sus¬ 
ceptibility  of  the  approaching  particle  to  jet  entrainment  is  accounted  for  by  the  term 
/>E,Uoodcos0  in  Equation  (1).  Particles  moving  away  from  the  Jet  or  to  its  side,  with 
momentum  not  directed  toward  the  Jet, are  less  tolerant  to  jet  entrainment.  The  term 


-P-£*  (JooSi/10  )  C 

/  ■+  £3  CjOS  &/  U  j 

takes  this  into  consideration  while  satisfying  the  Ricou-Spalding  solution  for  the  freejet 
case  when  =  0.  The  net  force  acting  on  the  jet  boundary,  as  a  result  of  the  pressure 
differential  around  the  jet,  is  accounted  for  by  a  crossflow  drag.  The  force  on  a  Jet 
element  of  unit  length  is 


f~p  *  \  pL&C05z&  d  (2) 

where  CD  is  the  cross  flow  drag  coefficient  of  the  jet. 
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Actual  UqoCOS0  Assumed 


FIGURE  1.  JET  EXHAUSTING  NORMALLY  INTO  THE  FREESTREAM 
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(1)  Equations  of  Motion 


The  entrainment  of  mainstream  fluid  into  the  Jet  and  the  pressure  forces  on  the 
boundary  of  the  Jet  govern  the  equations  of  motion  which  can  now  be  written,  using  the 
coordinate  system  defined  in  Figure  1. 

The  continuity  equation  may  be  written  as 

p  d/di  (AjUj)  -  £  (3) 

The  only  tangential  forces  acting  on  the  Jet  are  the  viscous  forces  resulting  in 
entrainment,  so  that  the  momentum  equation  for  the  jet  flow  in  the  tangential  direction 
is 

p  d/dt>  (AjUj1)  —  £ Uoc  ts//?G  (4) 

The  deflection  of  the  Jet  is  due  to  the  centripetal  action  of  the  pressure  force  on 
the  Jet  boundary  and  the  entrainment  of  mainstream  fluid  with  momentum  in  that  direc¬ 
tion.  The  force  equation  governing  the  curvature  of  the  jet  may  now  be  written  as 

pAjUjV/Z  =  EUoqcqsQ  t  Cd'v  pU*cos*4d  (5) 

«  pAj L  I  ♦  IX' YJ3iz 


where  R  is  the  radius  of  curvature  of  the  jet  centerline  and  the  primes  denote  differen¬ 
tiation  with  respect  to  Z.  Replacing  d/ds  by  cos  9  (d/dz)  in  Equations  (3)  and  (4) 
reduces  the  problem  to  one  of  finding  d,  Uj  and  X  as  functions  of  Z,  the  entrainment 
parameters  ,  Eg  ,  Eg  and  the  cross  flow  drag  coefficient  CD.  However,  a  func  ¬ 
tional  relationship  between  the  cross  sectional  area  Aj,  the  circumference  C  of  the 
jet  and  the  Jet  growth  must  be  established. 

Experimental  observations  show  that  there  is  a  region  in  which  the  jet  deforms 
from  its  initial  circular  cross  section  into  a  kidney-shaped  one  (Reference  8).  Once 
this  shape  is  attained,  the  jet  cross  section  remains  relatively  similar  (Reference  7). 
Since  it  did  not  appear  possible  to  treat  the  exact  jet  shape,  a  simplified  jet  shape,  an 
ellipse,  was  chosen.  Correlation  of  data  for  a  normally  exhausting  jet  indicates  that 
the  extent  of  this  region,  in  which  the  cross-sectional  deformation  occurs,  is  a  func¬ 
tion  of  the  jet  exit  to  mainstream  velocity  ratio  U^/Uo^  The  extent  of  this  region  may 
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be  expressed  by  0* Z/d0*.  W^o/Voo-  When  Z/do*.  3Ujo/Uoo,  observations  show  that 
the  best  fit  of  the  Jet  cross  soctlon  with  an  ellipse  is  one  with  a  ratio  of  minor  to  major 
axis  of  1/4.  Experimental  data  further  indicate  that,  in  the  development  region,  the 
ratio  of  minor  to  major  axis  decreases  linearly  with  distance  from  the  Jet  orifice 
(Reference  5).  Therefore,  the  geometry  of  the  Jet  may  bo  treated  in  two  regions: 

•  A  region  in  which  the  jet  deforms  to  an  elliptical  cross  section  and  the  ratio 
of  minor  to  major  axis  decreases  linearly  with  Z  from  1  at  Z/dQ-  0  to  1/4 
atZ/dQ--  H  = 

ZD*  ratio  of  minor  to  major  axis 

-  A  (Z/dDl  •-  B 

Substituting  boundary  conditions  above,  yields 

(Z/do) 

~  /-  l &(*/&)(  iWu.)  J 

then 

C  -  rrd  1/,  U-  W(Z/d*\  1*  J 1/2  J 

-  /id  f  !*  1 1  -  ( z/dc)l (W0o) l*]'"  I 

t  2  J  \  0«Z/do<H  (G) 

4  .  rr  [/-  g-  (z/do)  ]  £  1 

=  IT  [/-  (Z/da)<UWOjo)J  d.1  J 

From  the  preceding  expressions  it  can  be  Been  that  if  Ujo*0,  as  U,*~0,  the  development 
region  becomes  OiZ/^tooand  the  Jet  retains  its  circular  cross  section,  which  is  to 
be  expected  since  the  jet  is  now  exhausting  into  quiescent  surroundings. 

•  A  region  in  which  the  jet  retains  a  similar  cross  section 


= 

i/4 

\ 

c  - 

2.?Ad 

{  Z/do  >  H  (7) 

A]  - 

rrd2/iG 

J 
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(2)  Integration  of  the  Equations  of  Motion 

The  left  hand  side  of  Equation  (3)  may  be  rewritten  as 

p  d/ds  ( Aj  Uj )  -  p  d/da  (  AjUj1) 

Uj 


(8) 


now 


p  d/da  ( AjJJjI)  a  42  d/cb  (Aj Uj2)+p (AJUJ2) d/ds (J- ) 


Oj 


whore  from  Equation  (4) 


dj 


p  d/da  (Aj  U/)  ^  BUoo  sin& 


then 


(9) 


(10) 


(O  d/cb  (.400 

Oj 


/  £(J»  <s^6  -to4,U/\  d/di  l  Oj )  (ll, 
TTj  uj1 


substituting  into  (8) 


/  £LU  smO  -pA>  d/cb  (Uj )  ^  t 

TT;  ' 


(12) 


Substituting  for  Aj,  C  and  letting  d/ds  ■  cos  9  (d/dZ) 


prr  ( /-  %  £  J£t)  d*  <*/<&  Ml)  - 


(13) 


,  p£,  Uoo  dco6  9  >  £2((Jj‘(J**nO)/rd  p  ///-■%  (*M(cWl/jo)rJ 1/2 


*  |  Lk  ' 


n 


(14) 


Dividing  by  pdU*  and  letting  Z*  =  Z/do,  Uj*  -  Uj/Uj0,  and  m  =  Uj0/Uoo 


rr/4  ( /-  s4-£)  d case  d!d2 (U4)  = 


Ei  cos  9  + 


£7  (mUi*-Ain6)  IT  [  [sM^mUL 

/  +  £3 cosQ/Ufm  L  2  J  J  l 


if  d/dZ  =  (l/'do)  d/dZ*  and  d*  =  d/d0  then  the  left-hand  side  of  (14)  becomes 

ir/4  5/2  ^)  d*  COS  &  d/dz*  (  Oj/U<*>) 


rr/4.  ( /  -  s'2  )  d  m  d/dZ*  ( Of) 


Equating  and  solving 


Ml 

dZ * 


_  (  Etcos9  +■  Ez(r7i Uj  *~&n&) iT  1  /  )  1 

i+Ejcose/u/f n  «•  2  J 


x  (  - 


<nn9~  /tjUi* 


m  (/- 2S)  dWUj*cos9 


The  equation  for  dUj*/dZ*  applicable  for  the  developed  region  is  obtained  by  substi¬ 
tuting  Equation  (7)  for  Aj  and  C  in  (12) 


ml 

di* 


rrdlml0j*c:as9 


E,cosO  4-  £1 


z  ( tndi  *-  ds/n O)  224  /jidj'] 

/*  c.os&/o£m  1 
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Substituting  Equations  (6)  and  (7)  into  Equation  (3)  and  again  letting  Z*  =  Z/dQ, 
d*  =  d/dQf  Uj*  -  Uj/Uj0  and  m  =  Uj0/Uoo  yields 


Ml 

d£* 


y 


-  ol* 
majsO 


+  S/6  ird'Uj..  %d'z(l-%£)^£ 


(17) 


and 


dr 


t:,  cos  9  y  m  U.1- s//7&)  2. 14  _  £  mcQsUd*  dU* 

_ _ /  -y  <zas&/ u,*rrr~  & _  dZ* 

( #/& )  m Uj  f  <18) 


Substituting  Equations  (6)  and  (7)  into  Equation  (5)  leads  to 


(19) 


(  {  £,  t .b£n)  cos 9  /•  Ez( rtiU, *-s//j0)/7 

i  >o-^Y 

* 

_ 

1  /+  £zcos9/ US m 

2 

• 

cos  6 


and 


d*K*  = 

dr1 


(20) 


/£co£>& 

Fn,zd*u;2 


(Et  t  .5(d)cc>$Q  a 


£i(mUf-&*0)224  ! 
/  + 1$  cos  e/Ufm  | 


Equations  (15),  (17)  and  (19)  are  applicable  to  the  development  region  of  the  jet  and 
Equations  (16),  (18)  and  (20)  are  applicable  for  the  developed  region,  that  is,  when 
Z*/m>.  3. 
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With  the  additional  substitution 


Cos  & 


/ 


suo  = 


the  preceding  equations  are  seen  to  constitute  a  set  of  differential  equations  to  be 
solved  for  Uj*,  d*,  and  X*  as  functions  of  Z*  and  the  parameters  Ej  ,  Eg  ,  Eg ,  and  Cjj. 

Initial  conditions  at  the  jet  exit  are 

Z*  =  0. ,  X*  =  0. ,  Uj*  =  1. ,  and  dX*/dZ*  =  0. 

It  has  not  been  possible  to  integrate  these  equations  in  closed  form,  but  they 
have  been  integrated  numerically  with  the  aid  of  a  digital  computer.  The  system  of 
first  order  differential  equations  is  solved  by  means  of  a  fourth  order  Adams 
predictor/corrector  method  using  a  Runge-Kutts  starting  solution. 

(3)  Determination  of  Empirical  Parameters 

The  parameter  E2  is  determined  by  considering  the  results  for  a  freejet,  since 
E2  is  the  only  parameter  remaining  in  Equation  (1)  when  Uoo  =  0.  From  the  data  of 
Reference  6,  E2  was  originally  determined  to  be  .08,  on  the  basis  of  measured  entrain¬ 
ment  rates  in  the  fully  developed  region  of  a  jet  exhausting  into  a  quiescent  environ¬ 
ment. 


The  entrainment  of  ambient  fluid  into  a  jet  in  the  development  region  will  gener¬ 
ally  be  a  function  of  the  exit  characteristics  of  the  jet  and  therefore  will  usually  vary 
from  test  to  test.  The  entrainment  characteristics  of  the  jet  of  Reference  3  are  sh<wn 
in  Figure  2.  It  is  observed  that  entrainment  of  fluid  by  the  jet  increases  in  the  devel¬ 
opment  region  and  approaches  the  asymptotic  value  of  Reference  6. 

The  variation  of  entrainment  with  distance  from  the  jet  exit  has  been  incorporated 
into  the  jet  model  by  allowing  Eg  to  vary  in  the  development  region.  Knowledge  of  the 
entrainment  characteristics  for  the  static  case  is,  therefore,  used  in  determining 
crossflow  effects.  For  the  calculations  presented  in  this  report,  the  variation  of  E2 
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FIGURE  2.  ENTRAINMENT  CHARACTERISTICS 
FOR  THE  JET  OF  REFERENCE  3 
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In  the  development  region  to  Its  asymptotic  value  of  .08  is  based  on  the  information  of 
Figure  2. 

Cp  is  also  varied  in  the  development  region  to  take  into  account  the  change  in  the 
cross-sectional  profile  of  the  jet  in  this  region.  The  relationship  between  Cp  and  the 
cress  section  of  the  jet  which  is  employed  in  toe  computations  has  been  obtained  from 
an  unrest^tion  of  the  data  presented  in  Reference  23  and  is  given  by 


where  ZO  is  the  ratio  of  minor  to  major  axis  of  the  ellipse  representing  the  cross  sec¬ 
tion  of  the  jet  in  the  development  region.  This  expression  for  Cp  was  derived  from 
subcritical  test  data  and  thus  the  range  of  Reynolds  numbers  for  which  it  can  be  con¬ 
sidered  valid  is  approximately  107s  Re  $  10  .  Variations  in  Cp  do  not  appear  to  have 
a  great  effect  on  the  predicted  values  of  jet-induced  velocities,  so  that  use  of  the  ex¬ 
pression  outside  the  Reynolds  number  range  defined  above  is  not  expected  to  lead  to 
significant  errors. 

The  parameters  Ej  and  Ey  were  chosen  at  values  giving  good  correlation  bet¬ 
ween  experimentally  and  theoretically  determined  jet  centerlines  and  jet -induced  sur¬ 
face  pressures.  The  details  on  the  test  data  utilized  are  presented  in  Section  II.  2. 
This  approach  led  to  a  final  set  of  values  of  Ej  =.  45  and  £3  =  30. 


(4)  Calculation  of  the  Induced  Velocity  Field 

To  obtain  the  velocity  field  due  to  the  jet  interference,  the  entrained  fluid  is 
represented  by  a  uniform  sink  distribution  along  an  axis  normal  to  the  mainstream 
and  the  jet  blockage  effect  by  a  doublet  distribution  along  the  jet  centerline  as  shown 
in  Figure  3.  The  strength  of  the  doublet  distribution  is  obtained  from  the  1  /Z  term 
in  the  complex  velocity  potential  expansion  W(£)  for  the  two-dimensional  flow  past  an 
ellipse.  In  effect,  by  replacing  the  jet  with  a  doublet  distribution,  the  flow  past  an 
equivalent  circular  cylinder  is  being  considered. 


Consider  an  element  of  the  jet,  length  4s,  centered  at  (X,  0,  Z).  The  sink 
strength  per  unit  distance  in  the  rj  direction  will  then  be  given  by 


m 


St  Uoo  dcosO  /  fcz  ( Uj  -Uo&5ifl9)C 

/  +  hi  COS  &/ Uj 


<5b 

d 


(21) 


lfi 


FIGURE  3.  SINK  AND  DOUBLET  DISTRIBUTIONS 


FIGURE  4.  JET-INDUCED  DOWNWASH  VARIATION 
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0  * 


4*r 


The  iadheed  iriodlf  at  a  poiat  P  (Xp,  Yp,  Zp)  dae  to  a  sink 
will  be 

V'i  =  -  a^/dr 

-  -  /ft  <afct 

4wrx 

_  _ —nidi? _ 


The  velocity  component  in  the  X  direction  is  gives  by 

Us  -  -  ^/ar  =  -^0/Dr)(  Or/Dx) 


Performing  differentiation  with  respect  to  Xp 


of  length  M 

(22) 


Ui  =  ~  m  dq  dr  (23) 

4irrz  ax 

—  /W  do  (  X~Xp) _ 

4w  [  (Z -ZP)1  +  (X-Xp>* 1-  (  Yp^*]*z 


Similarly  the  induced  velocities  in  the  direction  of  the  axes  Y  and  Z  are  obtained 


V*  *  -  _  m  dr?  (Ye+ff) _ (24) 

*/r  L  (Z-7p)*  *  (i-xP)*  +  (Yp+4)*]51* 

Ws  -  ffidq  ( Z-Zf>) _  (25) 

4,7  f  (  Z-2p\L  t-(X~Xpy+(Yp  *%)*  J*2 
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F<ysatica  (23;  represents  Ike  rekdtjr  esxsspnvtzS.  ia  Ike  X-dfcreetaae  disc  to  at  sizk 
eknaess  dr„  kse^raiioc  over  ike  range  -  ;-<£  |  mli  pre  tke  wJocity  sadaeed  Sgr  a  je* 
eJM&eai  ■ef  kro g!&  <$. 


/ft  (X'X0)dff _ 

lL2  -2pf\ (X-Xpf  +(YpH?f] 


JZL  /*')&>)  f 

^/7  J  [(Z'ZpV  +(X"Xp)z+  (Tp-?)1]3'2 


Substituting  in  Equatioc  (26;  and  evaluating  from  —  d  2  to  d/2  yields 


<3us  - 


-J2?  X-Xp  _  x 

*/7  [(2-Zp^KX'Xp^J 


; _ Yp  ~  djz _ Xprdlz _ 

}  RZ-Zpf4-(X-Xp)2+(Yp~d/z)*]w  [(Z-Zp)zt(X-Xpy2^(YpidZZ^J,te 


Equations  (24)  and  (25)  can  also  be  integrated  to  yield 


-  -  j2L  J  _ ! _ 

4,t  I  [  (z-zPr  tcx-xpV+  crp-d/zY]^ 


j.  (Z~  Zp)i+  ( X'Xp^  +■  (Yp  7  d/z^J1'2 


221  (2  ~7p) 

4fT  [(Z-4>)S(X-Xp)x] 


Yp -d/z  Yp+d/2 

UZ-Zp)2  (-(X  -XpY  f  (Yp  tyH'1  ~  [CZ-V+CX'Xp)z+(:Yptd/z)zj” 


19 


F<pfflicw  (21>,  aad  |2S|»  gjsw  She  aihced  veiccsSjr  tcgyatato  aeS  Xp,  Yp,  ass4 
Zp  dm  so  rtplanrg  a  jet  fitattn*  ct  ietegsh  is  a y  She  ssaric  tfissrs&KSsoe. 

To  shss  she  b&ockage  effect  mass  he  an&Sesit  The  strength  off  she  «hx£&et  sfestzni&a- 
Sjoe  cSilized  So  rcpreseos  she  blockage  effect  off  she  |et  es  ofocaameef  ffroa  She  oetsspiex 
velocity  poCetoEial  for  she  mo-dEEmeEfiksea!  Qo*  past  as  e&epse  By  eqaaEisg  She 
sireaglfe  off  she  doshffes  So  she  coefficEesES  off  she  £  2  Sens,  she  flow  pass  aa  eqpuraleeis 
circular  cyisader  es  being  nesdered. 

For  She  tvo-SiaeasjOQsl  flcv  past  aa  ellipse 


W(Z) 


a  -  t> 


(3©) 


for  flow  along  a,  she  major  semi -axis.  Rotating  flow  SO  degrees  makes  a  she  minor 
S6roi-2xisf  b  nisior  ~2xis,  2  -  0  . 


W(2)  =  ill  (3fb)  ^  +  ( Zz-Cz)ta 

2  3+6 


+ 


■Z-  C?-cz)lft 

3-6 


From  binomial  expansion,  the  leading  two  terms  in  a  series  of  inverse  powers 

are 

W(Z)  z  II21I  U(aibVb)(‘/z) 

the  coefficient  of  the  1/jJ  term  is  equated  to  the  strength  of  the  doublet, p 
In  the  coordinates  of  Figure  1 
U  -  Uqb  cos  9 
a  -  minor  axis  of  jet 
b  =  major  axis  of  jet 

=  ±  Ooo  (3+b)(b) 


20 


CaasuSicrmg,  the  o»©>  gtsoeatsnc  repou  oi  the  jefl  dessriM  ta  FyawM  fS>  *af 
fl>  jr»ei4» 


Z*/es>-3 


Wife  fee  noeatic®  of  Figure  1  fee  irfaced  velocity  held  at  a  given  point  P  ({,  5,  *) 
doe  to  a  doable*  of  streagfe  *  at  U,  nt  *»  -  {0,  0,  (h  may  he  evalasted 

The  uiocity  potential  of  a  doublet  is  given  by 

0  -  ? £-&S  & 

4nr* 

where  r  -  [  f 2  +  S?  t  f2}  Vz 

CCtf0  -  /r 
then 

0  -  ^  £  <32) 

-Jf/7  iy+4****]** 

with  the  induced  velocities  in  the  4,  »7  and  {  directions  being  given  by 

U,  =  -  <30/df 

V,  -  -  b0/D^ 

W,  =  -  &&/  by 
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F*yfi«2aB  t3Zr  is  diflerertaad  to  yneld 


a,  =  ifz+4z  +  fr1)**  <235 

vt  =  -£>4$/jlI nr  L£1+jzz*  Z*x>5tz 


i*f.  =  2jjl  (.  ) 

-fr^yjTV  Sr9«*  [  !  L?****^')  [ 


<3S) 


From  the  vtl^uccship  be’.a«n  asd  X,  Y,  Z  in  Figure  J,  it  can  be  seen  ifeti 

=  u,  *iu(9  +  *4  cc?5@ 

<^vft  s  -  v', 

0W5  =  U,  CQ5<?  -  *\  mrl  (? 

^  as  (Xp~x)  S/tjO  *  ( 2p  -2)cc&@ 

SI  =  -'rp 

J*  •*  fXp-X)cas6>  (Zp~2)  5/n& 


where  4U g,  4Vg,  4Wg  are  the  induced  velocities  at  P(Xp,  Yp,  Zp)  due  to  a  jet  element 
of  length  as  centered  at  (X,  0,  Z),  and  u  is  obtained  from  Equation  (31). 

The  total  interference  velocity  at  a  point  may  be  determined  by  integrating  Equa¬ 
tions  (27)  -  (29)  and  Equati  on  (36)  over  the  extent  of  the  jet,  giving  the  total  induced 
sink  and  doublet  component  velocities.  These  component  velocities  are  then  summed 
to  give  the  total  interference  velocity.  In  treating  a  jet  exhausting  from  an  infinite 
wall,  an  image  system  of  the  singularity  distributions  is  utilized  to  satisfy  the  bound¬ 
ary  condition  at  the  wall.  The  flat  plate  pressure  coefficients  shown  in  Section  II.  2 
were  computed  in  this  manner.  Boundary'  conditions  for  a  wing  or  fuselage  in  the 
presence  of  a  flow  exiting  from  a  nozzle  are  discussed  in  Sections  IV  and  VIII. 

In  the  treatment  ol  the  blockage  effect  discussed  above,  the  jet  curvature  does 
not  contribute  anyth::. g  to  the  doublet  strength.  Three-dimensional  effects  are 
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accounted  tor  ifcroogb  the  doublet  position  rod  orientate*  and  by  calculating  the 
tadoced  velocity  field  from  the  three-dimensional  counterparts  to  the  doublets. 

For  the  Icnrer  velocity  ratios  where  the  jet  ^rratcre  is  quite  small, 

this  model  can  be  expected  to  be  a  good  representation  of  the  ,'ysical  situation.  As 
the  jet  curvature  increases  with  increasing  velocity  ratio,  ts  quasi  -two  dimensional 
approach  becomes  less  valid.  Calculations  of  induced  surface  static  pressures  for  the 
higher  velocity  ratios  (C*>  X' ^  -  .250)  did  not  correlate  with  test  data  as  well  as  did 
calculations  for  the  lower  velocity  ratios 

Attempts  were  made  to  improve  the  correlation, by  changing  the  entrainment 
expression  of  Equation  (i),  without  success.  This  suggested  that  some  other 
mechanism  was  responsible  for  the  lack  of  correlation  between  calculations  and  test 
data.  It  was  thought  that  jet  curvature  had  not  been  adequately  included  in  the  mode). 
Also,  it  was  argued  that  this  curvature  would  allow  relief  for  the  crossflow-  fluid 
ahead  of  the  jet  and  might  be  accounted  for  by  applying  a  distribution  of  sources 
along  the  jet  centerline.  One  would  expect  the  strength  of  the  sources  to  be  a 
monotonicaily  increasing  function  of  the  jet  c-cr  ter  line  curvature.  Subsequently,  it 
was  learned  that  Werner  and  Chang  (Reference  50),  using  a  slender  body  matched 
asymptotic  expansion  technique,  had  deduced  similar  results.  Calculations  of  jet- 
induced  static  pressures,  using  a  source  Distribution,  displayed  markedly  improved 
correlation  with  test  data. 

The  source  distribution  has  been  made  proportional  to  jet  curvature  and  best 
correlation  with  test  data,  over  the  range  of  velocity  ratios  .  1<  .  3,  has  been 

obtained  by  taking  the  source  strength  equal  to  three  times  the  jet  curvature. 

The  formulation  for  the  induced  velocity  components  at  arbitrary  points  given  in 
Equations  (21)  to  (36)  is,  of  course,  not  valid  for  points  inside  the  jet.  Figure  4  shows 
schematically,  the  variation  of  the  downwash  component  induced  in  the  plane  of  the  jet 
exit  by  a  flow  exhausting  from  a  nozzle  located  at  Y  =  0.  The  large  variation  across  the 
jet  in  the  freestream  direction  is  typical  of  the  computed  downwash  distributions.  Thus, 
a  certain  amount  of  care  must  be  exercised  in  choosing  the  control  points  at  which  in¬ 
duced  velocities  are  to  be  evaluated  for  use  with  the  methods  of  computing  forces  and 
mon  encs  in  Sections  IV  and  VIII. 

b.  Arbitrary  Jet  Direction 

Although  the  equations  of  motion  for  the  jet  model  developed  in  the  preceding 
discussion  are  written  for  a  jet  exhausting  normally  \ito  a  cross  flow,  u.ey  are  valid 
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FIGURE  5.  GENERAL  COORDINATE  SYSTEM 


FIGURE  6.  JET-CENTERED  COORDINATE  SYSTEM 


for  any  jet  in  a  local  coordinate  system,  oriented  with  the  X'-axis  in  the  direction  of 
the  mainstream  flow  and  the  X'Z’-plane  defined  by  the  mainstream  flow  vector  and  the 
jet  exhaust  vector.  The  direction  of  the  mainstream,  the  location  of  the  center  as  well 
as  the  initial  direction  r  f  the  exhausting  jet,  and  the  position  of  the  control  points  at 
which  induced  velocity  components  are  to  be  evaluated  may  be  arbitrarily  specified  in 
a  general,  fixed  coordinate  system  as  shown  in  Figure  5. 


A  local  coordinate  system,  centered  at  X. ,  Y. ,  Z.  is  then  established  as  shown 

A  J  3  J 

in  Figure  6(a).  The  direction  cosines  of  X'  with  respect  to  the  fixed  coordinate  system 

are  a  ,  B  ,  y  and  the  direction  cosines  for  Y'  and  are  determined  from  the  follow- 
0  *0  o 

mg  .elation ships. 


Y' 

A 

2' 


X'  X 


X' 


A. 

Y' 


where  Vj  is  the  unit  vector  in  the  initial  jet  exhaust  direction,  obtained  from  4>  and  ip, 
the  jet  exhaust  angles. 

This  establishes  the  co-planar  relationship  between  the  mainstream  flow  and  the 
jet  centerline  required  by  the  equations  of  motion  formulated  for  the  normally  exhaust¬ 
ing  jet. 


Initial  conditions  at  the  jet  exit  (shown  schematically  in  Figure  6b)  are  now 


Z*  =  0. ,  X*  =  0. ,  U j*  -  1. ,  d*  -  1. ,  and  dX*/dZ*  = 
3  the  angle  <5-,  the  ang! 
plane,  are  shown  in  Figure  6(b). 


1  -  COS70r, 


cos290 


V? 


A  A 

where  cos  60  =  Vj  •  Z 1 


Both  and  the  angle  <5-,  the  angle  specifying  the  exhaust  direction  in  the  exhaust 


Observations  during  the  experimental  investigation  of  Reference  3  showed  that  a 
jet  exhausting  into  the  direction  of  the  crossflow,  Jj>90°,  appears  to  deform  from  its 
initial  circular  cross  seciion  more  rapidly  than  a  jet  exhausting  with  the  crossflow, 
ij<90°.  Consequently,  the  development  region,  in  whio  •  the  jet  deforms  to  an  ellip¬ 
tical  cross  section  with  a  minor  to  major  axis  ratio  of  1/  4,  has  been  made  a  function 
of  the  jet  exhaust  angle  tfj.  For  6 j  >90°the  extent  of  the  development  region  is  defined 
to  be  H'=  H  (cos0o)  and  for  dj<90°  the  development  region  becomes  H'“  H/cos0Q 
where  H  -  .  3  U. _  /U~  is  the  extent  of  the  development  region  for  a  normally  exhausting 

JO 

jet,  defined  previously. 


The  geometry  of  the  jet  and  the  jet-induced  velocity  components  can  now  be 
computed  in  this  local  coordinate  system,  utilizing  the  equations  developed  for  a  jet 
exhausting  normally  into  a  crossflow. 
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2.  COMPARISONS  OF  SINGLE  JET  COMPUTATIONS  WITH  TEST  DATA 
a.  Jet  Centerlines 

Figure  7  shows  a  comparison  between  computed  and  experimentally  determined 
centerlines  for  normally  exhausting  jets,  with  velocity  ratios  of  Uw/U.0  =  .125,  .250 
and  Uoo/Uj0  =  .123,  233.  The  correlations  shown  in  Figure  7  were  instruments!  in 
the  determination  of  the  entrainment  parameters  Ej  and  Eg. 

Extensive  computations  were  carried  out  for  jets  at  various  velocity  ratios  and 
exhaust  augles.  A  comparison  of  computed  jet  centerlines  with  an  empirical  equation 
for  the  centerlines  of  jets  exhausting  at  angles  other  than  90  degrees  to  the  mainstream 
is  shown  in  Figure  8(a)  through  8(d).  These  correspond  to  Figure.,  12(a),  13(a), 

14(b),  and  15(a)  of  Reference  9.  The  solid  curves  represent  centerlines  computed  for 
the  specified  conditions.  The  broken  curves  represent  the  empirical  equation 

JL  =  ~  (tWUjof  /f\3.  _Z  rio/cSj  (37) 

do  cdj  \  do  /  do 

The  curves  superimposed  on  the  photographs  of  Figures  12  through  15  of  Refer¬ 
ence  9  show  the  fit  of  this  equation.  It  should  be  emphasized  that  the  broken  curves 
reoresent  an  empirical  rather  than  experimental  determination  of  the  centerlines. 

They  are  included  in  Figure  8  primarily  to  aid  in  orienting  the  computed  centerlines 
with  respect  to  the  photographs  of  the  jet  waxes  in  Reference  9. 

Figure  9  shows  a  comparison  of  experimental  and  theoretical  jet  centerlines  for 
a  jet  of  velocity  ratio  Uoo/Uj0  =  .120  exhausting  at  various  angles  into  the  mainstream. 
Since,  in  this  figure,  the  agreement  between  theory  and  experiment  is  not  a,-?  close  as 
might  be  desired,  a  comment  regarding  the  differences  is  in  order.  It  is  noted  that 
the  theoretical  curves  in  Figure  9  employ  values  of  entrainment  coefficients  based 
partly  on  correlation  with  the  experimental  data  of  Reference  7  (Figure  7).  Compari¬ 
son  of  the  experimental  centerline  tor  the  normally  exhausting  jet  of  Figure  9  with 
the  test  data  of  Figure  7(a)  reveals  differences  which  appear  to  be  inconsistent  with 
the  small  change  m  velocity  ratio.  In  view  of  these  experimental  differences  for  the 
normally  exhausting  jet,  the  correlation  for  the  jet  centerlines  shown  in  Figure  9  is 
considered  to  be  adequate. 

A  further  comparison  of  theory  with  experiment  for  jets  exnausting  at  an  angle 
other  than  90  degrees  is  presented  in  Figure  10.  Good  agreement  between  theory  and 
experiment  can  be  observed. 
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FIGURE  7.  CENTERLINES  OF  JETS  EXHAUSTING  NORMALLY  INTO  THE  FREESTREAM 
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FIGURE  8.  CENTERLINES  OF  JETS  EXHAUSTING  INTO  THE  FREESTREAM 

AT  V  ARIOUS  ANGLES 
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Penetration  into  Freestream,  Z/do 


FIGURE 


Distance  in  Freestream  Direction,  X/dQ 


9.  CENTERLINES  OF  A  JET  EXHAUSTING  INTO  THE  FREESTREAM 

AT  VARIOUS  ANGLES 
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FIGURE  10.  CENTERLINES  OF  JETS  EXHAUSTING  INTO  THE  FREESTREAM 

AT  AN  ANGLE  <J j  =  60° 
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b.  Induced  Pressures 


Extensive  computations  to  determine  induced  surface  static  pressure  distribu¬ 
tions  due  to  single  jets  exhausting  into  a  crossflow  weie  carried  out.  Comparisons 
between  theory  and  test  data  from  References  3,  12,  13,  and  14,  are  shown  in  Figures 
12  through  18  for  a  number  of  jet  configurations.  Figure  11  is  a  schematic  of  the 
single  jet  arrangement,  with  surface  pressure  taps  indicated,  employed  in  the  tests  of 
Reference  3. 

The  surface  static  pressure  variations  with  X/do  at  Y/do  =  0  and  1. 5  are  shown 
in  Figure  12  for  a  jet  with  a  velocity  ratio  l'ac/l'j0  =  •  125.  Test  data  from  References 
3,  12,  and  13  are  plotted  for  comparison  with  the  calculations.  At  Y/do  =  0,  the  dif¬ 
ferences  between  the  test  data  of  References  3  and  12  are  within  experimental  accuracy 
ahead  of  the  jet,  whereas  the  data  from  Reference  13  show  significant  differences. 

Behind  the  jet  there  are  discernible  differences  in  test  data,  due  in  part  to  the  unsteady 

nature  of  the  flow  in  this  region.  At  Y/do  =  1.5,  the  test  data  of  References  3  and  12  are 
seen  to  be  in  good  agreement.  There  is  some  discrepancy  between  these  data  and  the 
data  of  Reference  13  in  the  region  behind  the  jet. 

The  same  static  pressure  variations  for  a  jet  with  velocity  ratio  L„/rjo  -  .250 
are  shown  in  Figure  13.  At  Y/do  -  0,  the  lata  of  References  3  and  12  are  in  good 
agreement  ahead  of  the  jet,  while  significant  differences  are  noted  for  the  data  of 
Reference  13.  Differences  in  the  test  data  behind  the  jet  are  again  evident.  All  three  sets 
of  test  data  show  good  agreement  at  Y/do  =  1.5. 

The  purpi  oc  cf  making  these  comparisons  was  to  determine  whether  there  were 
discrepancies  among  data  from  different  tests  and,  if  so,  the  reasons  for  these  dis¬ 
crepancies.  The  data  of  Reference  3  are  in  good  agreement  with  the  test  data  of 
Reference  12  except  in  the  wake  region  behind  the  jet,  where  discrepancies  might  be 
expected  because  of  the  unsteady  nature  of  the  flow  in  this  area.  The  differences 
between  the  data  of  Reference  13  and  the  data  from  the  other  tests  could  be  due  to 
boundary  layer  effects.  The  plate  boundary  layer  was  thicker  for  tests  of  Reference 
13  than  for  the  others.  Boundary  layer  thickness  has  been  shown  to  have  a  significant 
effect  on  the  surface  static  pressure  distribution  (see  Figure  5  of  Reference  13)..  The 
correlation  among  the  sets  of  test  data  also  gives  confidence  in  the  data  of  Reference 
3  for  inclined  jets  and  multiple  jet  configurations  for  which  there  are  no  other  data 
available  for  comparison. 
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(a)  Station  Y/do  -  0 

FIGURE  12.  INDUCED  PRESSURE  VARIATION  DUE  TO  A  SINGLE  JET 
OF  VELOCITY  RATIO  Uoo/Uj0  -  .  125 


* 
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-2  0  2  4  6 

Distance  in  I-  reestream  Direction,  X/do 


(b)  Station  Y/do  =  1. 5 
FIGURE  12.  (Concluded) 
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Proatiuro  Coufflclunt, 


FIGURE  13.  INDUCED  PRESSURE  VARIATION  DUE  TO  A  SINGLE  JET 
OF  VELOCITY  RATIO  Ua,/Uj0  =  .  250 


-2  0  2  4  6 

Distance  in  Freestream  Direction,  X/do 


(b)  Station  Y/d0  =  1.5 
FIGURE  13.  (Concluded) 
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The  theoretical  pressure  varialJoess  show  good  agreement  with  test  data  is 
Figures  12  and  13  for  (he  rtgisss  ahead  of  the  jet  exit.  !t  might  be  looted  that  the 
agreement  between  theory  and  test  data  for  the  velocity  ratio  of  .  25©  was  greatly 
improved  by  utiiizatioc  of  the  source  dislribctffon  to  account  for  the  curvature  of  the 
jet,  which  was  discussed  oe  page  23. 

There  are  significant  differences  between  theory  and  test  data  in  the  wake 
regions  behind  the  jet.  The  singularity  representation  cf  the  jet  used  to  compute  the 
jet-induced  pressure  distribution  assumes  potential  flow  in  the  areas  external  to  the 
jet.  Physically  meaningful  results  can  only  be  expected  in  regions  where  the  flow- 
outside  the  boundary*  layer  is  potential.  It  is  possible  that  a  distribution  of  sources 
along  the  jet  boundary  may  be  used  to  represent  the  wake  effects  and  to  achieve  better 
correlation  with  test  pressures  in  the  wake  region.  Also,  it  should  be  noted  that  the 
wake  region  is  usually  small  with  respect  to  the  total  area  being  investigated  (see 
Figure  15  for  example). 

Figures  14,  15  and  16  show  the  pressure  distributions  in  the  plane  of  the  jet 
exits  around  jets  exhausting  at  angles  of  ij  =  120:,  90'  and  60  into  the  freest  ream. 

When  the  jet  is  directed  against  the  crossflow  t  ij  =  120  )  extensive  regions  of  negative 
pressure  result  primarily  ahead  of  the  jet.  As  the  jet  xhaust  angle  changes  to  90 
degrees  and  then  to  60  degrees,  these  regions  of  negative  pressure  are  observed  to 
move  downstream  (Figures  15  and  16).  Agreement  between  theory'  and  test  data  is 
noted  to  be  better  for  the  6j  =  90;  and  ij  =  60:  configurations  than  for  the  ij  -  120- 
configuration.  In  all  three  cases,  discrepancies  between  theory  and  test  data  exist  in 
the  wake  regions  of  the  jets  for  reasons  discussed  prt  viously. 

Additional  comparisons  between  computed  and  e  cperimental  induced  pressure 
distributions  around  a  single  jet  exhausting  normally  into  the  freestream  are  shown  in 
Figures  17  and  18,  for  velocity  ratios  of  .125  and  .089.  The  experimental  data  of 
Figure  17  provide  additional  confirmation  for  the  correlation  between  theory  and  experi¬ 
ment  observed  in  Figure  15.  They  also  indicate  that  the  theoretical  results  lie  within 
experimental  differences.  Very  good  agreement  between  theory  and  experiment  is 
evident  in  Figure  18. 


FIGURE  14.  PRESSURE  DISTRIBUTION  AROUND  A  SINGLE  JET 
EXHAUSTING  AT  AN  ANGLE  d<  =  120°  INTO  THE 
FREESTREAM  (U^/Ujo  =  .  125) 


FIGURE  16.  PRESSURE  DISTRIBUTION  AROUND  A  SINGLE  JET 
EXHAUSTING  AT  AN  ANGLE  6i  ^60°  INTO  THE 
FREESTREAM  (Uoo/Uu  =  .  125) ' 


—  T§st  Data  (Ref  13)  ' 
— Test  Data  (Ref  14) 

_ 1 _ 1 _ 

FIGURE  17.  PRESSURE  DISTRIBUTION  AROUND  A  SINGLE  JET 
EXHAUSTING  AT  AN  ANGLE  <5  \  =  90°  INTO  THE 
FREESTREAM  (Uoo/Uj0  = .  125) 
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FIGURE  18.  PRESSURE  DISTRIBUTION  AROUND  A  SINGLE  JET 
EXHAUSTING  AT  AN  ANGLE  <5  j  =  90°  INTO  THE 
FREESTREAM  (Uoo/Uj0  = . 089) 
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3.  MULTIPLE  JET  ANALYTICAL  MODEL 


The  basic  single  jet  analytical  model  has  been  applied  to  the  computation  of  the 
interaction  flow  field  due  tc  multiple  exhausting  jets.  A  multiple  jet  configuration  is 
treated  as  a  combination  of  discrete  jets,  with  each  jet  (including  jets  resulting  from 
the  coalescence  of  other  jets)  being  replaced  by  its  representative  singularity  distri¬ 
bution  to  obtain  the  induced  velocity  field.  In  the  development  of  the  two-jet  model 
described  in  detail  in  this  section  two  assumptions  were  made 

•  The  leading  (or  upstream)  jet  develops  independently  of  the  downstream  jet 

•  The  downstream  jet  exhausts  into  a  freestream  of  reduced  dynamic  pres¬ 
sure  which  it  sees  as  the  result  of  blockage  by  the  upstream  jet. 

Data  from  the  wind  tunnel  investigation  of  Reference  3  have  been  analyzed  to 
substantiate  these  assumptions  and  to  establish  quantitatively  the  empirical  relation¬ 
ships  utilized  in  the  two-jet  computations.  Details  of  this  analysis  ° re  also  given  in 
this  section. 

a.  Two-Jet  Computations 

Figure  19  shows  the  planform  of  three  jet  configurations  in  relation  to  the  main¬ 
stream  flow.  Arrangements  (a)  and  (c)  represent  limiting  cases.  The  arrangement  of 
(a)  allows  each  iet  to  develop  independently  to  the  point  where  the  growth  of  the  jets  in 
the  direction  normal  to  the  flow  causes  them  to  intersect.  The  configuration  shown  in 
(c)  places  the  downstream  jet  entirely  in  the  zone  of  influence  of  the  upstream  jet.  The 
downstream  jet  is  located  partially  in  the  zone  of  influence  of  the  upstream  jet  in  the 
arrangement  depicted  in  (b). 

Although  Figure  19  shows  the  relationship  of  the  jets  in  the  plane  of  the  jet  exits, 
as  an  illustration,  the  determination  of  the  degree  of  influence  of  the  upstream  jet 
(JET1)  on  the  downstream  jet  (JET2)  can  be  carried  out  for  each  element  of  JET1, 
as  shown  in  the  general  case  of  Figure  20. 

Plane  L,  defined  by  the  mainstream  flow  vector  and  the  diameter  of  JET1,  iden¬ 
tifies  the  element  of  JET2  influenced  by  a  given  element  of  JET1.  Plane  M,  which 
contains  the  mainstream  flow  vector  and  the  local  jet  velocity  vector  of  the  given 
element  of  JET1,  is  then  establ'shed.  The  intersection  of  the  plane  M  with  the  jef 
diameter  of  the  appropriate  element  of  JET2  is  computed  to  determine  the  extent  to 
which  the  two  jets  overlap. 
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FIGURE  19.  JET  EXIT  CONFIGURATIONS 


J 


FIGURE  20.  SCHEMATIC  OF  JET  INFLUENCE 
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An  expression  for  the  effective  dynamic  pressure  which  the  downstream  jet 
"sees"  as  a  result  of  the  blockage  by  the  upstream  jet,  based  on  analysis  of  experi¬ 
mental  data  will  be  derived  in  the  following  subsection.  Configurations  with  complete 
■overlap  between  the  two  jets  (jets  aligned  in  crossflow  direction)  and  partial  overlap 
will  be  considered,  and  the  extent  of  overlap  between  the  two  jets  will  be  a  principal 
parameter.  Thus  the  influence  of  the  upstream  jet  on  the  downstream  jet  is  intro¬ 
duced  into  the  computations  as  a  reduced  freestream  velocity,  Ue/U«>=  (Qc/flool^in  the 
continuity,  momentum  and  force  equations  governing  the  development  of  the  down¬ 
stream  jet. 

When  calculations  of  the  distance  bewteen  the  two  jet  centerlines  indicate  that 
the  jets  have  intersected,  initial  conditions  for  the  merged  jet  which  results  are  deter¬ 
mined  from  the  continuity  and  momentum  considerations: 


A  Uj,  t  Az  Ujj. 
(AiL/j,)LJj, x  + (AzUji)Uj2^ 
l 'A,Uj,)L Ij'x  HAzUJ2)Ujz y 
(AiUj,)Oj,£  (Azl/jz)  Uj2  2 

where 


Az  Uj3 
(AiCJjz) 

(  A  djSjj 
(AiL/js)  (Jjzz- 


(38) 


A.,  A,  A  =  area  of  JET1,  JET2  and  resulting  merged  jet,  respectively 
1^0 

Uji,  Uj2,  Uj3  =  jet  velocity  of  JET1,  -IET2  and  resulting  merged  jet,  respectively. 

Initial  conditions  for  the  coalesced  jet,  in  the  local  coordinate  system  centered 
at  the  point  of  intersection  of  the  two  jets  (taken  to  be  an  average  of  the  coordinates  of 
the  centerlines  of  JET1  and  JET2  at  intersection),  are  now 

Z*  -  0. ,  X*  =  0. ,  Uj*  =  1.,  d*  =  1.  and  dX*/dZ*  =  Uj3x»/Uj3zi 
with  the  velocity  ratio  being  given  by  Uoo/Uj0  =  Uoo/Uj3. 


Subscripts  x'  and  z'  denote  components  in  this  local  coordinate  system. 

These  initial  conditions  are  employed  in  integrating  the  set  of  differential 
equations  forlH.  d*  and  X*. 

The  cross-sectional  area  of  this  coalesced  jet  is  not  taken  to  be  circular  as  is 
normally  specified  for  an  exhausting  jet  at  Z*  -  0.  Experimental  observations  and  cor¬ 
relation  with  test  data  were  instrumental  in  establishing  the  geometrical  features  of  the 
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initial  cross  section  of  the  coalesced  iet  for  different  jet  configurations.  Details  of 
the  initial  geometry  of  the  merged  jet  will  be  discussed  in  the  following  subsection. 

The  velocity  field  induced  by  a  two-jet  configuration  can  now  be  determined  by 
replacing  each  jet  (including  the  coalesced  jet)  by  its  representative  singularity  distri¬ 
bution.  The  induced  velocity  components  due  to  each  singularity  distribution  are  addi¬ 
tive  at  every  control  point. 

b.  Determination  of  Empirical  Parameters 

Analysis  to  substantiate  the  assumptions  made  for  the  two-jet  computations  and 
to  establish  quantitatively  the  empirical  relationships  for  the  computational  procedure 
described  above  was  based  on  the  experimental  data  from  Reference  3. 

(1)  Effective  Dynamic  Pressure  for  Downstream  Jets 

The  blockage  effect  of  the  leading  jet  in  a  two-jet  configuration  may  be  deter¬ 
mined  from  the  test  data.  The  centerline  of  a  single  jet  exhausting  into  a  crossflow 
at  a  velocity  ratio  U^/U^  =  .125  together  with  test  data  from  References  5  and  12  are 
shown  in  Figure  21.  The  data  from  the  three  tests  are  seen  to  be  in  close  agreement. 
This  is  significant  since  all  three  tests  were  carried  out  using  different  values  for  the 
jet-exit  dynamic  pressure. 

Data  for  the  two-jet  configurations,  aligned  in  the  crossflow  direction,  are 
plotted  in  Figures  22,  23  and  24  for  the  three  different  spacings  between  the  jets. 
These  data,  together  with  the  single  jet  data,  may  be  used  to  derive  "effective  velo¬ 
city  ratios"  for  the  downstream  jets,  making  use  of  a  similarity  relationship  derived 
in  Reference  15.  It  was  determined  in  Reference  15  that  centerlines  of  single  jets 
with  different  velocity  ratios  are  similar  and  that  for  a  given  displacement  in  the 
crossflow  direction,  Ax/d0,  the  penetration  into  the  crossflow  normal  to  the  jet 
exit  is  inversely  proportional  to  the  velocity  ratio  (Uqo/U^)  .  By  applying  this 
similarity  relationship  to  a  number  of  displacements,  Ax/d0,  of  the  downstream 
jet,  it  is  shown  in  Figure  25  that  the  centerline  of  the  downstream  jet  in  a  two-jet 
configuration  is  the  same  as  the  centerline  of  a  single  jet  exhausting  into  a  crossflow 
of  reduced  dynamic  pressure. 

The  information  in  Figure  25  has  been  used  to  obtain  an  empirical  relationship 
for  the  dynamic  pressure  qe  which  the  downstream  jet  "sees,  "  as  the  result  of 
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crossflow  blockage  by  the  upstream  jet,  in  terms  of  the  crossflow  dynamic  pressure, 
qgp  and  spacing  between  the  two  jets,  s.  It  is  given  by 


5/dc.  - 1 
S/do  +  -  73' 


{  d/do  >  I  } 


(39) 


The  expression  for  fnn/ar>l'/j  given  in  Equation  (39)  is  now  used  as  a  limiting  value 
when  the  two  jets  are  aligned  in  the  crossflow  direction,  when  complete  overlap 
between  the  two  jets  exists. 

The  above  expression  has  been  derived  from  test  a  da  for  which  the  Reynolds 
number,  based  on  freestream  velocity  and  jet  exit  diameter,  was  approximately  2xl03. 
Comparison  with  jet  trajectories  from  References  7  and  12,  in  Figure  21,  indicates 
that  there  is  no  measurable  change  in  jet  trajectory  in  the  Reynolds  number  range 
103s  Re«2xl03.  However,  it  is  possible  that  a  large  change  in  Reynolds  number 
could  affect  the  jet  path  and  consequently  influence  the  formulation  for  effective  dynamic 
pressure  given  in  Equation  (39). 

When  the  two  jets  are  not  aligned  in  the  crossflow  direction,  an  effective  cross- 
flow  dynamic  pressure,  qa,  which  is  a  weighted  mean  of  q  and  q®,  is  utilized.  The 

p  e 

weighting  of  the  dynamic  pressure  is  determined  from  the  degree  of  overlap  between 
the  upstream  jet  and  the  downstream  jet  (see  sketch  below) 


Thus 

-K  -  x  f  -*)  jg™ 

cXz 


(40) 
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(2)  Jet  Decay  Characteristics 

The  decay  characteristics  of  the  single  jet  and  the  three  two-jet  configurations 
are  shown  in  Figures  26  through  29.  The  information,  obtained  by  analyzing  the  data 
of  Reference  3,  is  presented  in  terms  of  the  jet  dynamic  pressure  rather  than  the 
impact  pressure  which  has  customarily  been  used.  The  reason  for  this  is  that  the 
decay  of  jet  momentum  is  due  to  entrainment  of  ambient  or  crossflow  fluid.  Thus, 
for  jets  of  differing  pressure  ratios,  the  jet  momentum  flow  parameter  should  be 
invariant  rather  than  the  impact  pressure  parameter.  This  observation  is  supported 
by  the  data  presented  in  Reference  16. 

Figure  26  illustrates  that,  in  the  presence  of  a  crossflow,  a  jet  decays  at  a 
greater  rate  than  when  it  exhausts  into  a  quiescent  environment.  Figures  27,  28  and 
29  show  the  dynamic  pressure  decay  for  the  two-jet  configurations  at  three  different 
jet  spacings.  The  dynamic  pressure  has  been  obtained  from  the  impact  pressure  by 
assuming  that  the  static  pressure  throughout  the  jet  is  constant  and  equal  to  ambient 
pressure.  Small  variations  in  static  pressure  are  expected  but  they  should  not  affect 
the  general  conclusions.  The  distance,  s,  along  the  jet  centerline  is  measured  from 
the  center  of  the  jet  orifice  with  the  downstream  jet  exit  used  as  the  origin  for  the  jet 
resulting  from  the  coalescence  of  the  upstream  and  downstream  jet. 

From  the  data  in  Figures  27  through  29  it  is  deduced  that  the  rate  at  which  down¬ 
stream  jet  dynamic  pressure  decreases  with  distance  along  the  jet  centerline,  increases 
with  increased  spacing  between  the  two  jets. 

The  decay  characteristics  of  the  upstream  jet  in  each  of  the  three  two-jet  con¬ 
figurations,  together  with  data  for  the  single  jet,  are  shown  in  Figure  30.  It  is  deduced 
that  the  decay  of  dynamic  pressure  for  the  leading  jet  in  a  two-jet  configuration  is  inde¬ 
pendent  of  the  spacing  between  the  two  jets.  Since  the  centerline  of  the  leading  jet  is 
also  independent  of  the  jet  spacing  (see  Figures  22  through  24)  and  identical  to  that  for 
the  single  jet,  it  is  concluded  that  a  leading  jet  in  a  two-jet  configuration  develops 
independently  of  the  downstream  jet. 

This  result,  along  with  the  finding  that  the  downstream  jet  behaves  as  if  it  were 
exhausting  into  a  crossflow  of  reduced  dynamic  pressure,  is  important  in  the  applica¬ 
tion  of  the  analytical  model  to  the  multiple  jet  problem  since  it  permits  a  two-jet 
configuration  to  be  treated  as  a  combination  of  three  discrete  jets. 
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Inverse  Effective  Velooity  Ratio,  •^/qj0/qe 


(a)  Jet  facing  5  do 


d 

OS  10. 


(b)  Jet  Spacing  7. 5  do 


FIGURE  25.  EFFECTIVE  VELOCITY  RATIOS  OF  DOWNSTREAM  JETS 
IN  TWO-JET  CONFIGURATIONS  flWUio  =  .  125) 


\  A  Quiescent  Environment 

\  O  Crossflow,  Uoo/Uj0  =  .125 

A  I  I 


Distance  along  Jet  Centerline,  s/dG 

FIGURE  26.  JET  DYNAMIC  PRESSURE  DECAY 
FOR  A  SINGLE  JET 
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Dynamic  Pressure/Dynamic  Pressure  at  Exit,  qj/qj 


0  4  8  12  16  20 

Distance  along  Jet  Centerline,  s/do 


FIGURE  28.  JET  DYNAMIC  PRESSURE  DECAY  FOR  TWO  JETS 
AT  A  SPACING  OF  5  DIAMETERS 
(Uoo/Ujo  =  .  125) 
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amic  Pressure/Dynamic  Pressure  at  Exit,  qj/qj 
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FIGURE  30.  DECAY  CHARACTERISTICS  OF  LEADING  JETS 
IN  TWO-JET  CONFIGURATIONS 
(Uoo/Ujo  =  .  125) 
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FIGURE  31.  WATER  VAPOR  VISUALIZATION  FOR  TWO  JETS 
AT  A  SPACING  OF  2.  5  DIAMETERS  (U^/U^  -  .  125) 


Best  Ava\\ab\e  COPY 
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(3)  Geometry  of  Merged  Jet 

Test  observations  seemed  to  indicate  that  the  downstream  jet,  being  the  stronger 
jet  at  the  point  of  intersection,  due  to  shielding  by  the  upstream  jet  is  the  dominant 
jet. 

This  feature  of  the  two-jet  interaction  is  readily  discernible  in  Figure  31.  The 
dominant  influence  of  the  downstream  jet  on  the  characteristics  of  the  coalesced  jet 
results  in  the  high  degree  of  penetration  exhibited  by  the  coalesced  jet.  This  same 
feature  may  be  noted  in  Figures  22  through  24.  As  a  result  of  these  observations, 
geometrical  features  of  the  initial  cross  section  of  the  coalesced  jet  are  treated  as  a 
function  of  the  jet  orientation  angle,  £2  .  The  jet  orientation  angle  is  defined  as  the 
enclosed  angle  between  the  freestream  vector  and  the  line  joining  the  centers  of  the 
two  jet  exits. 

For  the  jet  spacings  of  the  experimental  investigation,  it  appears  that  best 
correlation  between  theory  and  test  data  is  obtained  by  employing  an  initial  circular 
cross  section  for  the  coalesced  jet  for  jet  orientation  angles  up  to  20  degrees.  For  a 
jet  o rotation  angle  £2  =  90°  (the  jet  configuration  shown  in  Figure  36)  the  geometry 
of  lie  jet  cross  section  of  the  coalesced  jet  is  initially  represented  by  an  ellipse, 
with  a  minoi  to  major  axis  ratio  of  1/2. 

At  this  time  the  functional  relationship  of  the  coalesced  jet  cross  st  tion  with 
orientation  an?rle  has  not  been  established  for  the  range  20°  <  £2  <  90°.  It  should  also  be 
noted  that,  while  the  above  discussion  is  limited  to  jet  configurations  with  a  maximrm 
spacing  of  7.5  jet  diameters,  the  importance  of  establishing  a  relationship  between 
orie  ation  angle  and  coalesced  jet  cross  section  decreases  for  large  jet  spacings,  since 
the  .ect  of  the  merged  jet  on  induced  pressures  in  the  plane  of  the  jet  exit  diminishes. 
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4.  COMPARISON  OF  TWO  JET  CALCULATIONS  WITH  TEST  DATA 

Computations  for  a  range  of  two-jet  configurations  were  effected.  Comparisons 
between  theory  and  test  data  from  Reference  3  are  presented  in  Figures  32  through 
42. 

‘  '  I 

Comparisons  between  predicted  and  experimentally  determined  pressure  distri¬ 
butions  in  the  plane  of  the  jet  exits  around  two-jet  configurations  with  spaeings  of 
2.5,  5  and  7.5  jet  diameters  are  shown  in  Figures  32  through  34.  In  the  computations 
shown  in  Figures  32  through  34',  "effective  velocity  ratios"  of  .058,  .087  and  .098 
were  utilized  for  the  downstream  jets  at  the  respective  spaeings1  of  2.5  dQ, '  5  dQ  and 
7.5  d  .  These  are  based  on  an  actual  velocity  ratio  of  .  125  and  the  relationship 
between  crossflow  dynamic  pressure,  q^,  and  the  reduced  dynamic  pressure  experi¬ 
enced  by  the  downstream  jet  as  given  by  Equation  (39).  For  the  spacing  of  2. 5  diam¬ 
eters  it  may  be  noted  that  the  pressure  distribution  is  similar  to  that  induced  by  a 
single  jet  exhausting  into  a  crossflow,  although  the  regions  of  negative  pressure  are, 
more  extensive.  As  the  jet  spacing  increases,  the  downstream  jet  is  seen  to  have  its 
own,  discrete  effect  on  the  pressure  distribution.  'Isobars  of  C^.=  -.3  are  noted 
around  the  leading  jet  and  the  downstream  jet  in  Figure  33;  and  in  Figure  34,  which 
shows  the  larger  jet  Spacing,  isobars  of  Cp  =  -.  2,  3,  are  noted  around  each  jet. 

Generally  good  agreement  between  theory  and  test  data  is  discernible  for  all 
three  jet  configurations,  It  should  be  pointed  out  that  relatively  small  differerices  in 
measured  pressure  coefficients  at  stations  of  constant  Y/dQ,  when  the  variation  with 
X/dQ  is  small  (i.e.,  when  the  curve  is  very  flat)  can  lead  to  what  appear  to  be  much 

larger  discrepancies  in  the  position  of  the  isobars  shown  in  the  figures. 

1 

Theory  and  test  data  for  the  centerlines  of  the  three  two-jet  configurations  above 
are  shown  in  Figures  22  through  24.  Generally,  the  prediction  for  the  centerline  of 
the  coalesced  jet  does  not  achieve  the  degree  of  penetration  exhibited  by  the  experi¬ 
mental  data. 

Figure  35  shows  a  sketch  of  a  two-jet  configuration,  indicating  the  location  of 
pressure  taps,  and  defines  the  sideslip  angle  fi.  The  variation  of  induced  pressure 
with  angle  of  sideslip  for  this  configuration  is  shown  in  Figures  37  through  39. 

Figure  37  shows  comparison  between  theory  and  experimental  data  for  induced 
pressure  variation  with  X/dQ  at  Y/dQ  -1.5  and  3.0  for  zero  sideslip.  Figure  38  shows 
the  same  comparison  for  a  sideslip  angle  &  -  20'.  The  effects  of  positive  sideslip 
are  seen  to  be  an  increase  in  the  magnitudes  of  the  induced  pressures  in  the  region 
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ahead  of  the  leading  jet  exit  and  a  decrease  in  the  region  immediately  behind  it. 

Figure  39  displays  the  effects  of  a  sideslip  angle  0-  -20°  on  the  above  pressure  dis¬ 
tributions  .  Negative  sideslip  is  seen  to  reverse  the  two  effects  discussed.  It  results 
in  a  decrease  in  the  magnitudes  of  induced  pressures  in  the  region  ahead  of  the  leading 
jet  exit  and  an  increase  in  the  region  immediately  behind  it.  These  trends  in  the  test 
data  are  also  evident  in  the  computed  results.  For  both  positive  and  negative  sideslip, 
the  downstream  jet  is  discernible  as  a  more  discrete  influence  on  the  pressure  distri- 

I 

button  when  compared  to  the  zero  sideslip  configuration.  This  feature  of  the  induced 
pressure  distribution  is  predicted  quite  well  by  theory. 

( The  correlation  between  theory  and  test  data  in  Figures  37  through  39  and  the 
effects  of  positive  and  negative  sideslip  are  representative  of  results  obtained  for  a 
two-jet  configuration  with  a  jet  spacing  of  7. 5  diameters  also.  Figures  40  and  41  show 
the  pressure  variation  with  X/dQ  at  zero  sideslip  and  at  0  =  20°  for  a  two-jet  configu¬ 
ration  with  a  spacing  of  2.5  diameters.  The  correlation  between  theory  and  test  data 
is  seen  to  be  very  good  for  the  zero  sideslip  condition  of  Figure  40.  Noticeable  differ¬ 
ences  between  theory  and  test  data  exist  for  the  0  -  20°  case  of  Figure  41. 

At  zero  sideslip,  (qc/qoolyi  =  -46  for  the  downstream  jet,  as  computed  from 
Eq.  (39)  with  s  =  2. 5  dQ.  This  means  the  downstream  jet  "sees"  a  low  crossflow 
dynamic  pressure  and,  consequently,  does  not  exert  a  strong  influence  on  the  induced 
flow, field.  The  assumption  that  the  upstream  jet  develops  independently  of  the  down¬ 
stream  jet  is  therefore  justified,  despite  the  close  jet  spacing.  The  good  agreement 
between  theory  and  test  data  in  Figure  40  supports  this  conclusion. 

At  0  =  20°,  [qc/qaJ%  =  .  93  for  the  downstream  jet ,  as  given  by  Equation  (40). 
The  jet  now  has  a  stronger  influence  on  the  induced  flow  field.  This  strong  influence, 
together  with  the  close  jet  spacing,  makes  the  assumption  that  the  upstream  jet 
develops  independently  of  the  downstream  jet  no  longer  representative  of  the  physical 
situation.  Further  mutual  interference  effects  between  the  two  jets  must  be  included 

i 

to  .improve  correlation  between  theory  and  test  data. 

Figure  36  shows  a  schematic  representation  of  a  two-jet  configuration  with  a 
spanwise  jet  spacing  of  7.5  diameters.  Theoretical  and  experimental  pressure  varia¬ 
tions  with  Y/dQ  are  shown  for  five  stations  of  constant  X/dQ  in  Figure  42.  As  expected, 
correlatioh  between  computed  answers  and  test  data  is  better  at  stations  ahead  of  the 
jets  than  at  stations  immediately  behind  the  jets.  The  wake  region  of  the  two  jets  is 
qleariy  defined  at  station  X/dQ  =  1.5  in  Figure  42(d).  Correlation  again  improves 
farther  behind  the  two  jets  where  the  importance  of  wake  effects  tends  to  diminish. 
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FIGURE  33.  PRESSURE  DISTRIBUTION  AROUND  TWO  JETS 
AT  A  SPACING  OF  5  DIAMETERS 
(Uoo/Ujo  =  .  125) 
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FIGURE  34.  PRESSURE  DISTRIBUTION  AROUND  TWO  JETS 

AT  A  SPACING  OF  7.  5  DIAMETERS 
(U®/Uj0  =  .  125) 
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FIGURE  35.  TWO-JET  CONFIGURATION  WITH  SIDESLIP 


FIGURE  36.  SPANWISE  TWO- JET  CONFIGURATION 
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(a)  Station  Y/do  =  1.5 


Distance  in  Freestream  Direction,  X/do 

(b)  Station  Y/do  =  3 

FIGURE  37.  INDUCED  PRESSURE  VARIATION  FOR  A  TWO-JET  CONFIGURATION 

AT  ZERO  SIDESLIP  (Spacing  =  5  do,  Uoo/Uj0  =  .125) 
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Pressure  Coefficient,  CD  Pressure  Coefficient, 


Pressure  Coefficie.'  Cp  Pressure  Coefficient 


(a)  Station  Y/do  =  1.5 


Distance,  X/d0 
(b)  Station  Y/d0  -  3 

FIGURE  30  INDUCED  PRESSURE  VARIATION  FOR  A  TWO-JET  CONFIGURATION 
AT  SIDESLIP  fi  -  -20°  (Spacing  5  do,  U«»/r,0  .123) 
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Distance  in  Freestream  Direction,  X/do 

(a)  Station  Y/d0  =  1.5 


Distance  in  Freestream  Direction,  X/iIq 


(b)  Station  Y/d0  -  3 


FIGURE  40.  INDUCED  PRESSURE  VARIATION'  FOR  A  TWO-JET  CONFIGURATION 
AT  ZERO  SIDESLIP  (Spacing  -  2.5do,  LVl'jo  -  .125) 


(a)  Station  Y/d0  =  1.5 

FIGURE  41.  INDUCED  PRESSURE  VARIATION  FOR  A  TWO-JET  CONFIGURATION 
AT  SIDESLIP  0  =  20°  (Spacing  =  2. 5 do,  Uoo/Uj0  =  .  125) 
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(b)  Station  Y/do  =  3 
FIGURE  41.  (Concluded) 
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5.  APPLICATIONS  TO  MORE  COMPLEX  CONFIGURATIONS 


As  outlined  previously,  the  computational  procedure  treats  each  jet  in  a  multi jet 
arrangement  as  a  single  jet,  with  shielded  downstream  jets  exhausting  into  a  cross- 
flow  of  reduced  dynamic  pressure. 

This  approach  does  not  restrict  jets  in  a  multi  jet  arrangement  to  have  the  same 
velocity  ratios.  Figure  43  shows  the  computed  pressure  distribution  around  a  tvc  ,« t 
arrangement  with  the  leading  jet  at  a  velocity  ratio  of  .  215  and  the  downstream  jet  at 
a  velocity  ratio  of  .  150.  Test  data  are  also  shown  for  comparison. 

The  extension  of  the  computational  procedure  to  the  treatment  of  more  cor  pi  ex 
jet  configurations  was  effected  for  a  three-jet  arrangement.  In  the  three-jet  configu¬ 
ration  sketched  below,  s^  is  the  distance  between  Jet  No.  1  and  Jet  No.  3  and  s9  is 
the  distance  between  Jet  No.  2  and  Jet  No.  3. 

Jet#l  Jet  #2  Jet  #3 


The  effective  freestream  dynamic  pressure  into  which  Jet  No.  2  exhausts  is 
obtained  from  the  relationship  for  the  two-jet  configuration  given  by  Equation  (39),  where 
s  is  the  spacing  between  Jet  No.  1  and  Jet  No.  2  above.  The  crossflow  dynamic  pres¬ 
sure  which  Jet  No.  3  "sees"  is  given  by 


^/do  -  1  Sddo  -1  [  [Si/do  7'] 

3t/cfc,'T.7?_  SVd  t  +  ,75i 


(41) 


Equation  (41)  accounts  for  shielding  by  both  of  the  leading  jets  on  Jet  No.  3.  If 
Sg  is  very  large,  there  is  effectively  no  shielding  and  ~  1-  If  Sj  is  very 

large,  only  Jet  No.  2  provides  shielding,  and 


lE*  «  _  - 1 

'  5.i/dc  f  -  75 

Both  jets  contribute  to  the  reducbon  in  dynamic  pressure  when  Jet  No  1  and  Jet  No. 
2  are  close  to  Jet  No.  3. 
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FIGURE  43.  PRESSURE  DISTRIBUTION  AROUND  TWO  JETS 
OF  DIFFERENT  THRUST  LEVELS 
(Upstream  Jet,  Uoo/Uj0  =  .  215) 
(Downstream  Jet,  U<x>/Uj0  =  .  150) 
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FIGURE  44.  TIIREE-JET  CONFIGURATION 
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Figure  44  shows  a  sketch  of  the  three-jet  configuration  tested  in  Reference  3, 
The  computed  centerlines  for  this  jet  arrangement  are  shown  in  Figure  45.  A  com¬ 
parison  between  predicted  and  experimental  pressure  variations  at  spanwise  stations 
of  1. 5  and  3.0  jet  diameters  is  presented  in  Figure  46.  In  the  computations  shown  in 
Figures  45  and  46,  "effective  velocity  ratios"  of  .058  and  .069  were  utilized  for  the 
two  downstream  jets  spaced  at  2.5  dQ  and  7.5  dQ  from  the  leading  jet,  respectively. 
These  are  based  on  an  actual  velocity  ratio  of  .  125  and  the  relationships  for  reduced 
dynamic  pressure  of  Equations  (  J9)  and  (41).  Good  correlation  between  theory  and 
data  is  noted. 
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SECTION  m 


MAPPING  METHOD  FOR  ARBITRARY  CROSS  SECTION 

To  use  the  transformation  method  and  the  nonlinear  body  aerodynamics  method 
which  are  described  in  Sections  IV  and  VI,  it  is  necessary  to  have  a  method  of  mapping 
an  arbitrary  wing  or  body  cross  section  ir.lo  a  unit  circle.  It  is  the  purpose  of  this 
section  to  describe  the  mjthod  developed  during  this  investigation  to  obtain  such  a 
mapping. 

There  are  numerous  methods  available  to  map  specific  shapes  into  a  unit,  circle 
such  as  the  Joukowski  transformation  for  wings.  It  is  also  possible  to  map  polygons 
by  the  Schwarz-Christoffel  tran^^rmation,  but  this  transformation  maps  the  polygon 
into  a  straight  line,  r,ot  into  a  circle. 

The  only  suitable  transformation  for  mapping  an  arbitrary  section  into  a  ,nit 
circle  is  that  developed  by  Skulsky  in  Reference  46.  This  method  uses  a  procedure 
similar  to  the  Schwarz-Christoffel  transformation.  This  method  would  be  suitable  for 
the  mapping  required  by  the  programs  developed  during  this  study  except  for  the 
following  restrictions  on  the  method.  First,  the  method  replaces  the  cross  section  by 
one  made  up  of  straight  line  segments  giving  only  an  approximation  to  the  desired 
shape.  Second,  the  mapping  procedure  is  an  iterative  one  and  the  question  of  con¬ 
vergence  in  the  iteration  arises.  Third,  a  large  number  of  terms  in  the  mapping  are 
required. 

These  restrictions  are  not  so  severe  as  to  prevent  the  use  of  the  Skulsky  method 
for  the  present  study  program  methods,  but  fhe  method  is  restricted  to  symmetrical 
shapes  and  this  restriction  means  that  general  airfoils  could  not  be  treated.  For  this 
reason,  it  was  decided  to  develop  a  mapping  function  which  would  be  more  suitable 
for  the  present  study. 

The  mapping  method  developed  in  this  section  presents  an  alternative  approach 
to  the  computation  of  a  mapping  function.  This  method  is  not  restricted  to  symmetri¬ 
cal  sections  and  does  not  require  any  iterative  procedure.  The  mapping  is  obtained 
by  finding  the  potential  for  ar.  arbitrary  section  which  corresponds  to  a  simple  vortex 
flow  about  the  corresponding  circle.  The  potential  is  equated  to  that  for  a  circle 
and  equivalent  points  are  located  on  the  mapping  circle.  Knowing  the  point-to-point 
correspondence,  the  mapping  function  is  then  obtained. 
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1.  COMPUTATION  OF  COMPLEX  POTENTIAL  FOR  VORTEX  FLOW 

The  method  used  to  obtain  a  mapping  function  is  to  develop  a  potential  function 
about  an  arbitrary  shape  which  satisfies  the  boundary  condition  of  no  flow  through  the 
body  and  behaves  as  a  simple  vortex  flow  at  large  distances  from  the  body.  This 
expression  is  written  in  body  plane  coordinates  (see  Figure  47)  and  is  compared  to 
a  similar  flow  about  a  circle  expressed  in  circle  plane  coordinates,  which  is  a  well- 
known  function. 


The  complex  potential  of  a  vortex  flow  about  a  circle  is 


W  *  +  itfl 


(42) 


The  velocity  in  the  circle  plane  is 


dW 

<y? 


~uc  +  L  Vc 


A 


(43) 


To  obtain  the  corresponding  velocity  about  an  arbitrary  section  it  is  necessary 
to  obtain  an  analytic  function  of  Z  which  satisfies  the  boundary  conditions  at  the  body 
and  behaves  like  i/Z  at  large  distances  from  the  body.  It  is  possible  to  obtain  such  a 
function  by  writing  the  velocity  in  the  body  plane 

(44) 


where  f(Z)  is  some  function  whici'  approaches  j  as  Z  approaches  infinity  and  which  is 
chosen  to  satisfy  the  boundary  conditions  at  the  body.  By  writing  the  velocity  in  terms 
of  polar  coordinates  i.e., 


dW 

dZ 


(45) 


satisfying  the  boundary  condition  at  the  body  reduces  to  finding  a  function  f(Z)  which 
matches  a  specified  variation  of  a  around  the  section. 


Assuming  temporarily  that  the  section  to  be  mapped  is  smooth,  that  is  has  no 
corners,  the  log  of  Equation  (45)  can  be  taken  to  give 


Jtn>  la  -  -  hn  E  t  Jm  -f  ( Z) 


(46) 


and  In  ard  a  are  conjugate  functions, 
series  wYh  no  nonnegative  powers,  i.e., 


di  *  it  i  <3;  +  ibj 


Thus,  if  In  f(Z)  is  taken  to  be  a  Laurent 


('7) 
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the  coefficients  of  this  expansion  can  be  chosen  to  obtain  any  desired  variation  of  a 
about  the  body. 


This  givea  the  expression  for  the  complex  velocity  ir  '  ie  section  plane  as 


i/be  (48) 

When  the  section  to  be  mapped  has  comers  a  will  be  discontinuous  and  it  is  necessary 

to  obtain  an  expression  which  has  similar  discontinuities.  Assume  that  a  corner  exists 

\k 


on  the  section  to  be  mapped  at  Z  =  Z^.  Then,  the  factor 


P#) 


m 


where 

JL  m - 45£L_ 

W+Aam 

will  have  a  singularity  at  Zm  wi  ,h  a  discontinuity  in  angle  of  This  can  be  shown 

by  comparing  the  limits  on  the  t”"'  ^ides  of  the  corner.  Letting  Z  =  Z^  -  4  Z 


(/-  !=)*”-, 

/  Aseia)*i»_ 

[As 

1  B  / 

{  -2/  12  )  l 

■  »t  J 

rj*m  (a.-i t-w) 

Again,  letting  Z  =  +  4  Z 


subtracting  these  two  expressions  gives 

i  Am  (&+•  =  ~i  dam 


(49) 


which  is  the  result  sought.  The  magnitude  of  this  factor  gives  an  infinite  velocity  for 
an  outside  comer  and  zero  velocity  for  an  inside  comer  which  is  to  be  expected. 
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The  complex  velocity  for  an  arbitrary  section  can  thus  be  written 


(50) 


m-f 


where  M  is  the  number  of  corners.  This  expression,  being  written  wholly  in  terms  of 
the  section  coordinate  Z,  can  be  evaluated  directly. 


Defining, 


the  log  of  this  function  can  be  written 

M  /M  \ 

InVfy-ia  =  Rb  +  i  (H  Am  Pm- Tl-to)+q(2)  (51) 

m-l  vn*l  ' 

The  imaginary  part  is 

M 

a-Tr-O)*  Z  Am  =  (g(Z ))  (52) 

/»■/ 

The  left  hand  side  of  this  equation  is  a  continuous  and  periodic  function  of  s  and  can  be 
evaluated  at  any  point  on  the  body.  This  makes  it  possible  to  evaluate  the  coefficients 
of  g(Z)  to  satisfy  che  boundary  conditions  to  any  degree  of  accuracy  desired. 


Letting 


M 


t-a-7r-a)  *Z  Ampm+tfn(g(zj) 

(!•/ 


where 


a  f  _  V  t>n  COS  nu>-dn  Sin  niv 
(9  (Z))  -  2.  „ 

Rb 

and  setting 


E  represents  the  least  square  error  in  fitting  a .  Differentiating  with  resoect  to  each 
of  the  a's  and  b's  and  setting  the  results  to  zero 


dbj 


gives  a  set  of  linear  equr tjons  solvable  for  the  a's  and  b's. 
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Once  these  coefficients  are  determined,  the  real  part  of  Equation  (51)  can  be 
evaluated  and  obtained.  Then,  comparing  the  potential  functions  in  the  circle  plane 
with  that  about  the  section  gives 


G 


27rjvbds 

fvbds 


(53) 


where  6  is  the  angular  distance  about  the  circle.  This  gives  a  relationship  between  the 
points  on  the  section  and  points  on  the  mapping  circle  which  permits  a  mapping  to  be 
obtained. 


2.  COMPUTATION  OF  DERIVATIVE  OF  MAPPING  FUNCTION  WITH  CORNERS 


Once  it  is  known  which  point  on  the  section  maps  into  which  point  on  the  circle 
it  is  relatively  simple  to  obtain  the  mapping  function,  and  this  relationship  has  been 
obtained  above. 


as 


7  he  derivative  of  the  mapping  function  can  be  written  as  a  Laurent's  series  in  t 


T* 


♦ 


-*■ 


dn 

<n 


(54) 


the  l/t  terra  being  zero  for  a  closed  section.  It  can  also  be  shown  that 

dZ  _  Js_ 
c/C  "  M© 


(55) 


when  the  mapping  is  evaluated  at  the  mapping  circle.  Then  the  coefficients  of  the 
mapping  can  be  evaluated  as  follow.*.  The  radius  of  the  mapping  circle  can  be  deter¬ 
mined  by  setting, 


or 


271 i  —  4 


(a-Q-Tr/t) 


ds 


or 


i(a-9-rr/2)ds 


(56) 
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similarly 


/< 


f*"' If  <*« -*»*•<, 


zirid„*f,c-<e  i(n")e~iee ifr-e-Wf  ,c  € iode 
dn  »  ^j~fe  i^-Oe-Jr^ 


(57) 


Since  it  has  been  shown  above  how  to  evaluate  6  as  a  function  of  s  it  is  possible  to 
obtain  the  derivative  of  the  mapping  function. 

It  is  possible  at  this  point  to  integrate  dZ/'ds  directly  and  to  obtain  the  mapping 
function.  When  the  section  to  be  mapped  has  comers,  however,  the  region  in  the 
vicinity  of  the  comers  is  not  mapped  accurately  unless  an  extrem  »ly  large  number  of 
terms  in  the  mapping  function  is  retained.  For  this  reason,  it  is  desirable  to  rewrite 
dZ/dC  in  such  a  manner  that  the.  comer  singularities  are  expressed  explicitly.  This 
requires  that  dZ/d£  be  written  as 


It  is  now  possible  to  expand  the  product  factors  of  Equation  58  by  binomial  expan¬ 
sion,  complete  the  expansion  by  further  multiplication  and  the  terms  equated  with  like 
terms  of  Equation  -54).  This  results  in  a  set  of  expressions  which  can  be  solved  for 
the  A’s.  Equation  (58)  allows  a  satisfactory  mapping  near  the  corner  points  with  a 
limited  number  of  terms  but  cannot  be  integrated  directly  to  give  the  mapping  function 
analytically.  It  is  possible  to  integrate  (58)  numerically,  howevr  r. 


3.  COMPUTATION  OF  MAPPING  FUNCTION 


The  computer  programs  written  for  the  transformation  method  and  the  nonlinear 
body  aerodynamic'’  method  do  not  permit  the  use  of  the  derivative  of  the  mapping  func¬ 
tion  but  require  that  the  nmpning  function  itself  be  known  and  be  written  in  the  form. 
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+ . .  ±n. 

Kn 


(59) 


To  obtain  the  mapping  function  in  this  form,  the  expression  of  Equation  «5g)  is  expanded 
to  obtain  the  derivative  of  the  mapping  by  expanding  each  of  the  terms 


by  a  binomial  expansion.  These  expressions  are  then  multiplied  together  and  also 

multiplied  by  the  summation  term  in  Equation  (58).  This  results  in  an  expression  for 

dZ/dC  which  can  immediately  be  integrated  analytically  into  the  form  of  Equation  (59). 

This  integration  allows  calculation  of  all  the  coefficients  of  Equation  (59)  except  for  cq. 

To  obtain  cq  it  is  necessary  to  compare  the  original  and  the  mapped  section  and  shift 

the  mapped  curve  to  best  reproduce  the  original  cross  section.  The  amount  of  shift 

renuired  is  then  the  value  of  c  . 

o 

4.  SAMPLE  CALCULATIONS  OF  MAPPINGS 

This  method  has  been  used  to  map  several  sections,  figure  48  foi  example  is  a 
mapping  of  an  extremely  thick  Joukouski  airfoil  which  was  mapped  since  the  exact  point- 

to-point  mapping  correspondence  could  be  determined.  Shown  on  the  figure  are  both  a 
complete  mapping  from  the  physical  coordinates  of  die  airfoil  giving  an  excellent  repro¬ 
duction  of  the  airfoil  with  the  exception  that  some  error  is  introduced  near  the  trailing 
edge,  and  a  mapping  starting  /ith  the  true  6  distribution  which  shows  even  better 
agreement,  hi  the  complete  mapping  twenty  terms  of  the  expansion  for  g(Z)  were  used. 

Figure  -»9  shows  a  mapping  of  a  65-010  airfoil  using  the  above  method.  In  this 
case  only  ten  terms  of  the  expansion  were  used  so  that  the  agreement  in  this  case  is 
not  quite  so  good,  although  still  satisfactory.  To  obtain  this  mapping  it  was  necessary 
to  place  a  singularity  within  the  airfoil  at  the  center  of  the  leading  edge  radius.  With¬ 
out  this  additional  singularity  an  unsatisfactory  mapping  was  obtained. 

Figure  50  shows  a  mapping  of  a  muc>  more  difficult  section,  the  section  being 
that  of  a  T-38  section  in  the  region  where  the  ;et  inlets  are  adjacent  to  the  body.  This 
section,  which  is  symmetrical,  has  five  discontinuities  on  each  side.  The  mapping 
was  made  using  twenty  terms  in  the  expansion  for  g(Z).  The  resultant  agreement 
between  the  true  section  and  the  mapped  section  is  very  good. 
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SECTION  IV 


TRANSFORMATION  METHOD 

1.  INTRODUCTORY  REMARKS 

Under  the  assumption  of  inviscid  flow,  the  force  and  moment  on  an  immersed 
body  arise  from  the  pressure  alone.  Thus,  the  problem  is  that  of  determining  the 
pressure  distribution  on  a  wing  or  fuselage  with  one  or  more  jets  exhaust¬ 
ing  into  a  uniform  crossflow.  The  flow  is  assumed  to  be  in  viscid,  incompressible 
and  irrotational.  The  wing  or  body  configuration  can  be  fairly  general  as  long  as  its 
cross  sectional  profile  can  be  mapped  into  a  circle  on  a  two-dimensional  space.  The 
governing  differential  equation  for  the  motion  is  the  three-dimensional  Laplace 
equation  (viscosity  has  been  implicitly  included  in  the  jet  model).  The  boundary 
condition  is  the  usual  flow  tacgeacy  condition. 

In  addition  to  these  assumptions,  we  further  assume  that  the  presence  of  the 

body  has  negligible  influence  on  the  jet.  Thus,  the  jet  flow  field,  ?  ,  mav  be  computed 

v 

independently  and  used  in  the  formulation  of  the  boundary’  conditions  for  the  present 
p .  oblem. 

The  total  velocity  vector,  *|  ,  may  th°n  be  written  as- 


}  =  *  ?;  -  %  *  % 


(60) 


where  i^is  the  freestream  velocity,  the  perturbation  velocity  due  to  the  presence 
of  the  body  in  the  freestream.  and  ^is  the  perturbation  velocity  due  to  the  body 
presence  in  the  jet  induced  flow  field.  The  problem  of  determining  ^  is  here  termed 

-J  -a 

the  power-or.  problem,  that  of  determing  LL*  ,  the  power-off  problem  and  that  of 
determining  the  power-effect  problem..  Emphasis  will  be  placed  on  deter- 

nrning  the  power  effect  since  we  assume  that  methods  already  exist  for  determining 
the  power-off  pressure  distribution. 


i  c-s,  the  problem  is  reduced  to  that  of  solving  for  the  perturbation  velocity 
,  in  the  presence  of  a  body  and  a  given  jet  induced  flow  field,  ^ .  The  governing 
differential  elation,  as  nr.mtioeed  above,  is  the  three-dimensional  Laplace  equation 


dz‘  *V 


(61) 


wbere  x,  y  and  z  are  body -oriented  coordinates  (Figure  51)  and  <P  is  the  perturbation 
velocity  potential.  The  perturbation  velocity  components  are  iL,  =  -  —$■ .  VL  =  -  , 

'  ji  - 

40  0  -  .  The  perturbation  velocity  vector  may  then  be  written  as: 


^  =  Aji  r  ^  j  r  (62) 

The  problem  is  then  to  find  a  solution  of  Equation  {61),subject  to  the  boundary' 
conditions  that 


(  fj  -  •  n  —  O  on  the  body  surface. 


(63) 


and  grad  <p  vanishes  at  infinity.  The  vector  n  refers  to  the  body'  surface  outward 
normal. 

It  is  possible  to  obtain  a  numerical  solution  to  the  problem  by  applying  Green’s 
Theorem  using  distributed  singularities  on  the  body  surface.  The  basic  formulation 
of  the  problem  would  be  similar  to  that  given  by  Hess  and  Smith  in  Reference  17. 
This  approach,  although  giving  a  numerically  exact  solution,  is  not  convenient  for 
the  present  investigation  for  two  reasons.  The  first  reason  is  that  a  large  amount  of 
computing  time  is  required  to  generate  solutions  for  practical  problems.  The  second 
reason  is  that  a  large  number  of  control  points  are  required  to  describe  the  body 
geometry. 

Thus,  in  this  investigation,  a  quasi-two-dimensiona!  method  has  been  developed 
which  is  ar  approximate  solution  to  Equation  61  and  satisfies  the  boundary  conditions 
in  an  approximate  way.  The  technique  consists  of  essentially  two  parts  which  corre¬ 
spond  to  two  steps  in  the  computation  and  are  railed  the  segment  method  and  the 
three-dimensional  modification.  Step  fi).  or  the  segment  method,  consists  of 
obtaining  an  approximate  three-dimensional  velocity  potential  by  de5er mining  exact 
solutions  of  the  two-dimensional  Laplace  Equation  at  a  number  of  1.  >dy  stations.  In 
Step  (2),  this  velocity  potential  is  modified  to  better  satisfy  the  three-dimensional 
Laplace  Equation. 


‘•9 


Step  (1),  or  Jie  segment  method,  is  an  application  of  slender  body  theory  in 
which  the  x -derivatives  are  neglected  in  favor  of  derivatives  b  the  other  two  direc¬ 


tions. 


Thus,  the  governing  differencial  equation  reduces  to 


=  0 


(64) 


The  boundary  conditions  are  then 


ft  y 


on  each  v  -  z  section 


and  grad^p^  vanishes  at  infinity,  where  V  is  an  outward  normal  to  each  y-z  section. 
This  boundary  value  problem  may  new  be  solved  by  first  mapping  the  contour  of 
every  y-z  section  imo  a  circle  and  then  using  the  method  described  in  the  following 
subsection.  The  potential  so  determined  involves  x  as  a  parameter  and  may  be 
written  symbolically  as  <p  (y,z.x).  These  potentials  could,  however,  be  consolidated 
to  form  a  three-dimensional  velocity  potential.  This  latter  potential  will  not,  in 
general,  satisfy  the  three-dimensional  Laplace  equation.  To  remedy-  this,  a  techni¬ 
que  is  developed  as  Step  (2). 

Step  (2)  is,  in  effect,  an  alternative  method  of  including  some  three-dimensional 
effects  in  the  problem  in  place  of  the  usual  pi^ -cedurc  e«i«pi  jyed  in  the  slender-body 
theory.  (A  general  description  of  the  slender-body  theory  may  be  found  in 
Reference  18.)  The  present  technique  is  based  on  the  (act  that  failure  of  the  consolidated 
potential  to  satisfy  the  three- dimensional  Laplace  Equation  implies  local  violation  of 
•r  aserv^tion  of  mass.  To  reme-jv  this  deficiency,  a  distribution  of  residual  sources 
aid  sinks  are  included. 

Inclusion  of  these  residual  sources  and  sinks  supplements  the  old  velocity  field, 

9j  -  grad9D  <p(y,  z;x),  and  is  denoted  by  The  existence  of  the  body  will  also 
have  to  be  recognized  in  this  supplementary  field.  To  do  this,  we  again  use  the 
segment  method.  The  new  flow  field  after  Step  (2)  may  then  be  written  as: 


f'=  %  -  <p (66) 

This  process  could,  in  principle,  be  repeated.  Each  time  when  the  process  of  three- 
dimensional  modification  is  involved,  two  more  terms,  similar  to  ^  and  grad9p  <ps  , 
would  be  added  to  the  above  expression. 
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2.  computation  of  fuselage  pressure  distribution 


As  stated  above,  the  present  technique  comprises  two  steps,  the  segment 
method  and  the  three-dimensional  modification.  The  former  can  be  used  inde¬ 
pendently,  whereas  the  latter,  being  a  modification  procedure,  must  be  used  in  con¬ 
junction  «:t!>  the  first  step.  Some  details  of  these  two  steps  will  now  be  discussed. 


a.  Segment  Method 

(1)  The  Boundary  Function 

The  method  is,  in  effect,  a  numerical  extension  of  the  procedure  described 
in  Reference  19.  The  problem  considered  in  Reference  19 is  that  of  determining 
the  complex  velocity  potential  for  a  cylinder  moving  two-dimensionally  in 
infinite  fluid  at  rest  at  infinity.  The  motion  of  the  cylinder  is  described  by  a 
velocity  of  translation  U  and  an  angular  velocity  w.  Let  C  be  the  contour  of 
the  cylinder  cross-section  and  O  the  origin  of  coordinates  about  which  the 
cylinder  is  rotating 


Z -plane 


£ -plane 


PHYSICAL  PLANE 


MAPPED  PLANE 
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Suppose  the  region  outside  C  in  the  Z-pIane  is  mapped  conformally  on  the  outside  of 
the  unit  circle  1^1  =  1  in  the  5 -plane  by  the  relation 

Z  =  f(S)  <«> 

with  the  points  at  infinity  in  the  z  and  $  -plane  corresponding.  Again,  following 
Reference  19,  the  normal  velocity  of  the  cylinder  along  the  outward  normal  is 

(litas*  -  gat) S/nO' ~  (Usinoc  rx.to^cos6/ 


in  which  if  is  the  inclination  of  the  tangent  to  the  x-axis.  Thus, 
-fjr  =  (Uc*sa-$u>)&  -  (Usina  +xu>)  jff 

Integrating  along  the  cylinder  boundary,  we  get 

iff  =  Uxs.n<x  -  Ug&s*  t-  f)  +  6 


(68) 


(69) 


where  B  is  an  arbitrary  constant  which,  without  loss  of  generality,  may  be  set 
equal  to  zero. 

Thuo,  y  is  equal  to  the  imaginary  part  of  the  function 

{(Iff-)  *  -U2e 


so  that  finally,  on  the  boundary  of  the  cylinder,  the  stream  function  <i>  is  such  that 

2i{f--Uie  tU?e  i-  iujff  (70) 

Let  d-  e°  denote  a  general  point  on  the  unit  circle  in  the  £-plane.  Then,  on 
the  unit  circle.  Equation  (70'  gives 

2i  =  &(if)  -  -  Uf(f)e~l'a+  uT(j)e*  *■  j)  (71) 

This  'unction,  3(d),  is  callrJ  the  boundary  function  by  Milne-Thoir.son. 
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If  the  boundary  function  is  expanded  in  a  power  series  in  <f  ,  we  can  write 
6(<f)  ~  B,(J)  +  &J<t) 

where  6,-fd’)  ontains  the  negative  powers  of  d  only  so  that  is  holo- 

morphic  (  is  finite,  single  valued  and  has  a  finite  single  valued  differential 

coefficient)  outside  the  unit  circle  and  vanishes  at  infinity. 

Thus,  in  terms  of  the  complex  velocity  potential  function  W(<f ) ,  Equation  (71) 
may  be  written  in  the  form 


U(<f)  -  w(%)  =  B,(<r)  +  BZU)  (72) 


Let  i  be  the  circumference  of  the  unit  circle.  Then  Equation  72  may  be  used 
to  obtain 

_L  f  1  [  H(7T)d*  _  JL  f  B,WWrf  ±  \ 

ivi)  <*-$  %ni jf  iviJ  6~S  J  "V-5 

—  < 

Now,  W($)  and  B((5)  are  holomorphic  outside  /,  while  W(  "j  )  and  B2(f) 

are  holomorphic  inside  1 .  Therefore,  if  $  is  outside  i  ,  application 
of  Cauchy's  Theorem  shows  that  the  second  and  fourth  integrals  vanish  and 
the  first  and  third  give 

V($)  =  (73) 

Thus,  Equation  (73)  gives  the  complex  velocity  potential  as  a  function  of  jj  . 
Elimination  of  5  between  Equations  67  and  73  yields  the  complex  potential 
as  a  function  of  Z  . 

In  summary,  to  determine  the  complex  potential  W(Z) ,  it  is  necessary  to 
know  the  mapping  function  of  Equation  67  and  the  stream  function  on  the 
boundary  of  the  cylinder  (Equation  (70)).  The  boundary  function  in  the  circle 
plane  may  then  be  determined  and,  hence,  that  part  which  contains  only 
negative  powers  of  £  .  The  complex  potential  is  then  immediately  deter¬ 
mined. 
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(2)  Application  of  Boundary  Function  to  the  Segment  Method 

Although  a  cylinder  under  translation  and  rotation  was  given  as  an  example 
of  the  boundary  function,  this  approach  may  be  used  to  solve  more  general  problems 
in  which  $  may  be  determined  around  a  contour  *or  which  the  mapping  into  a  unit 
circle  is  known. 

Consider  an  elongated  body  with  its  longitudinal  axis  aligned  in  the  direction 
of  the  free3tream  (Figure  51).  Let  the  velocity  components  in  the  x,  y  and  z 
direction  be  denoted  by  u,  •/,  and  w  respectively.  The  body  is  divided  by  cuts  as 
shown  and  the  volume  adjacent  to  each  cut  represents  a  segment. 

From  the  computation  of  the  jet  flow  field  (see  Section  II),  we  know  what  the 
jet  induced  velocity  components  of  u,  v  and  w  along  the  cross  sectional  boundary  of 
each  cut  would  be  if  the  body  were  not  present.  Due  to  the  presence  of  the  body 
segment,  we  must  satisfy  the  boundary  condition  that  no  fluid  penetrates  the  surface 
of  the  body  segment.  In  the  method  of  segments,  this  boundary  conditioi  is  not,  in 
general,  satisfied  exactly.  It  is  satisfied  approximately  by  moving  the  cross  section 
of  each  cut  in  such  a  way  that  the  velocity  component  normal  to  the  cut  is  equal  and 
opposite  to  that  induced  by  the  jet.  Thus,  application  of  the  boundary  condition  does 
not  account  for  growth  of  the  body  in  the  x-direction. 

The  method  developed  in  Section  III  may  be  used  to  determine  the  mapping 
function  of  Equation  67  so  that  the  problem  of  obtaining  the  complex  potential  and, 
hence,  velocity  and  pressure  distribution,  is  reduced  to  one  of  determining  the 
stream  function  on  the  given  boundary. 

Follov/ing  the  procedure  for  determining  a  boundary  function,  we  first  transform 
the  cross  sectional  area  of  each  cut  to  a  circle  by  the  mapping  function 


Z  =  t,  *  K 


bn 

r 


> 


where  Z  =  y  +  iz  (physical  plane), 


re 


id 


(mapped  plane). 
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FIGURE  51.  REPRESENTATIVE  CONFIGURATIONS 


The  complex  constants  b^,  b^, .  .  -  bQ,  which  signify  the  geometry  of  the  cut, 
are  determined  by  the  procedure  discussed  in  Section  III  of  this  volume.  The  complex 
potential,  on  the  mapped  plane  due  to  the  motion  of  each  cut  is  assumed  to  be 

W  =  4>  T  i  'F  =  a%&n%  +  J1-  +  Jy  +  '  •  ‘  +  '  (74) 

To  compute  the  boundary  function  for  every  cut,  we  write 


When  we  perform  the  integration  of  this  formula  along  a  circle  where  the  variable 
r  remains  unchanged,  it  reduces  to 


0 


Putting  ^  ;n  a  form  convenient  for  computation  yields 


where  v.  and  w, 
J  j 


are  the  jet  induced  velocity  components. 


(75) 


If  we  set  the  upper  limit  to  2  n  ,  this  integral  gives  the  net  flow  of  fluid 
through  the  boundary.  This  may  not  be  zero,  due  to  the  fact  that  only  two  velocity 
components  have  been  accounted  for  instead  of  three.  This  source  or  sink  inside  the 
boundary  is  accounted  for  by  the  logarithmic  term  in  Equation  74  . 

We  first  determine  the  coefficients  a^.  .  .  an  in  Equation  (74)  by  consider¬ 
ing  the  periodic  part  of  Equation  (75).  Thus,  we  form 


y  -=  Y  + 


_0_ 

o?7T 


2  TT 


„ri£ 

~%Je 


)d0 
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The  boandary  function  ^  can  wnr  be  expanded  In  the  form  of  a  Fourier  aeries 
by  Ibe  usual  method 


%  =  ^  ^rtCcS/?e 

n=c 


Sinnd . 


The  boundary  function  may  now  be  used  to  determine  the  coefficients  a^,  .  .  . 
an  in  Equation  (74-  and,  hence,  determine  the  potential  for  the  flow  which  will 
cancel  the  normal  component  of  velocity  calculated  by  the  jet  program. 

The  first  term  in  Equation  (74)  represents  the  source  and  circulation  potential,  and 
will  be  considered  later.  The  coefficients,  a^,  a2,  .  .  .  an,  of  the  remaining  terms 
may  be  determined  by  equating  them  to  A  and  Bn-  Setting  a^=  a  j  +  i  a.’’ ,  &2  =  aj>  * 
i  a",  .  .  .  ,  we  obtain 


rb  ’  "  ft 


A  -  3* 

ht  rx  > 


(76) 


or 


a,  =  ~  rb  8,  +  1  Yb  > 

o-z  -  -rb2B2  +  t 


where  is  the  radius  of  the  mapped  circle. 

To  evaluate  the  coefficient  of  the  logarithmic  term  in  Equation  (74),  we  write 

~  &o  Q+t  ?  +■  iO-oO  . 


For  a  body  without  circulation,  aQ  is  real  and  is  given  by 


ITT 


J 


0 
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Thai  far,  the  complex  coefficients  in  Equation(74)hare  all  been  determined, 
which  sets  the  flow  through  every  cross  sectioc  of  the  cut  to  zero  in  a  predetermined 
external  flow  field.  The  complete  complex  potential  is  then  the  summation  of  Equa¬ 
tion  74  and  the  potertial  induced  by  the  exhausting  jets.  Once  this  potential 
is  known,  it  is  a  relatively  simple  matter  to  write  a  computer  program  to  calculate 
the  pressure  distributions  along  the  boundary  of  every  cut,  which  result  in  the  power 
effect  by  the  segment  method.  Sample  calculations  discussed  in  Volume  n. 

b.  Three-Dimensional  Modification 


By  establishing  the  real  part  of  die  complex  potential  in  the  physical  plane,  we 
obtain  the  velocity  potential  at  every  point  exterior  to  the  body  on  the  planes  passing 
through  the  cuts.  These  velocity  potentials  all  exhibit  the  same  form  except  that 
their  coefficients  are  different  from  one  plane  to  another.  In  other  words,  every 
coefficient  is  actually  a  function  of  x.  To  show  this  parameter  x  explicitly  in 
Equation  74 ,  we  write  the  velocity  potential  at  station  x,  with  the  help  of  Equation  (76) 


♦  1771 


in  which  r  and  6  are  functions  of  y  and  z  given  by  the  mapping  function.  By  consoli¬ 
dating  these  expressions  from  every  station,  we  could  construct  a  three-dimensional 
velocity  potential.  This  potential  does  not,  however,  satisfy  the  three-dimensional 
Laplace  equation.  To  remedy  this  deficiency,  we  proceed  to  Step  (2). 

In  Step  (2),  we  first  divide  the  region  of  the  flov,  field  into  a  network  of  small 
parallelepipeds  whose  "centers"  are  situated  on  the  planes  passing  through  each  cut 
(Figure  52)  and  then  compute  the  flux  through  the  surfaces  of  each  parallelepiped  by 
using  the  available  velocity  potential.  To  account  for  this  flux  approximately,  we 
have 


Qi 


=  ( 
\ 


in  which  the  partial  derivatives  are  r-’aluated  at  the  center  of  the  parallelepiped  and 
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A  V  refers  so  its  soBaease.  Cws  iemuegsg  sfee  poeeseiiBs  at  sMs  amt  Bss  jarsi^fcwrasg 
sfaftioas  to  be  ,  3ca£  tfee  velocity  cggspogieggs  to  be  %■=  ~ 

aad  .  weobcaiffl 


f  <fi  .  t- 

_  .9  I -  -  —  - - -  4- ; -  Hi 


The  quantities  in  ti  e  bracket  represent  an  approxiroatioo  for  tk  partial  deriniire 
j  These  redace  to  tbe  central  difference  formula  for  stations  of  equal 

intervals.  The  sons  is  always  equal  to  zero,  since  is  the  potential 

obtained  uy  »k  segment  method.  However,  the  net  flux  Q_,  computed  in  this  master,  is 
not  in  general  equal  to  ze>-c>.  This  implies  that  fluid  has  been  created  or  destroyed 
within  the  parallelepiped.  In  other  words,  the  Laplace  equation  has  a  residual  term, 
since  the  continuity  equation  and  the  Laplace  equation  are  equivalent.  To  counter¬ 
balance  this  residual  term,  a  source  of  the  same  strength  but  of  opposite  sense  is 
placed  at  the  center  of  the  parallelepiped,  so  that  the  three-dimensional  Laplace 
equation  is  approximately  satisifed  in  its  neighborhood.  Repeating  the  process  for 
each  parallelepiped  yields  a  network  of  residual  sources  and  sinks  in  the  space.  These 
sources  and  slides  create  a  new  field  which  will  modify  the  old  one.. 

A  schematic  diagram  for  the  parallelepiped  network  is  shown  in  Figure  52. 

Though  there  are  only  two  layers  and  few  parallelepipeds  indicated,  it  is  understood 
that  for  computational  purposes  we  need  at  least  four  or  five  layers  to  cover  a  space 
extending  several  body  diameters  outward.  To  illustrate  the  effects  of  residual 
sources  and  sinks,  let  A  be  an  arbitrary  point  on  the  body*  surface  (Figure  52)  with 
the  coordinates  xq,  vq  and  zq.  Let  the  centers  of  the  surrounding  parallelepipeds 
be  x.,  y.,  and  z^.  The  residua!  velocity  components  at  r’oint  A  are  then  given  by 
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in  which  w  :  -Q.  and  M  refers  to  the  total  number  of  the  residual  singularities  whose 
presence  may  affect  the  velocity  dist  ibutions  at  Foir.t  A.  A  schematic  arrangement 
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T  COMPUTATION  OF  *ATNG  PRESSURE  DISTRIBUTION 

A  finite  wing  :s  also  a  asree-ffinjeESioas!  both-,  bu*  st  has  a  distinct  property  of 
generating  lift  in  a  uniform  flw.  Consequently,  the  above  procedu re  is  applicable 
with  the  additional  requirement  that  the  Kutta-Joukowsky  cowiitwr.  be  satisfied  at 
the  trailing  edge.  To  accomplish  this,  consider  the  following: 
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If  it  is  assumed  that  the  trailing  edge  occurs  at  9  •>.  it  become.1  :  '-e=sar;  ;o  set 

the  tangential  velociU  c-omponen*  equal  to  7cm  here,  i .e..  to  set  the  imagmar>  part 
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to  compote  die  dmnrasfa  components  at  eien  section.  Tfte  synsxri  b  denotes  use 
span  length  or  that  portion  of  the  span  length  over  which  the  power  effect  is  significant. 
Then  a  new-  flow  field  on  the  fictitious  wing  is  const  ratted  by  subtracting  '*  from  the 
modified  w  components  is  step  (2).  This  constitutes  the  modification  to  account  for 
the  s pairwise  variation  of  the  circulation  around  the  wing  after  inclusion  of  the  residual 
source  distribution. 


It  is  noted  here  that  there  are  two  arbitrary  element-:  olved  in  the  above 

procedure.  These  are: 


(a>  The  factor  1/3  in  Equation  78  is  an  added  empirical  factor  to  the 
original  equation. 

(b)  Since  only  the  lifting  line  theory  is  invoked,  all  the  component6  at  any  given 
section,  regardless  of  the  chordvise  position,  are  reduced  by  an  equni 
amount  equal  to  w'(y). 
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4.  COM PTCATHBS  A3fi>  BERT T5 
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f?S  £25223"  tittsi-  hscsizr}  irsmr-rss  is  F :trxr  series. 

(gj  Derer-rsiae  ie  roeffanesss  of  the-  c-j-sspfex  pecesnaai  s»  each  seetioc  irocr 
are  Focrier  ooefScaescs. 

(fe>  Compete  she  pressure  cEs-rribcsioe  frees  rbe  compeex  poseKial. 

fi)  Dfiercne  aje  bret  aiti  aesect  by  irsegrasing  tat  pressure  distribeatiou. 

If  the  Sir«-fi!EeE5=<a2:  nwificattia  is  needed,  we  sfcp  saep  fs)  and  isrlcde  ah e  iofitnr- 
;iej  steps: 

(f»  Form  a  network  of  parallelepipeds  >n  the  exterior  space. 

(k)  Compute  the  net  flux  in  each  parallelepiped. 

(l)  Place  a  source  of  the  same  strength  but  of  opposite  sense  at  the  center  of 
each  parallelepiped. 

(ml  Determine  the  velocity  distributions  along  the  sectional  boundaries  induced 
by  the  residual  sources  and  sinks. 

(n)  P.epeat  steps  fe)  through  (i)  to  obtain  the  force  and  moment. 

Extensive  calculations  have  been  carried  out  for  two  wings  and  two  bodies 
with  lifting  jets  in  a  uniform  crossflow.  The  ’parameters  involved  in  these  computa¬ 
tions  are:  freestream  to  the  jet  velocity  ratio  .  angle  of  attack,  sidcsl'p  angle, 
location  of  the  jets  and  the  number  of  lifting  jets.  The  calculated  results  have 
been  compared  with  wind  tunnel  test  data.  The  agreement  has  been  fairlv 
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SECTION  V 


INLET  FLOW  FIELD 

Islet  flow  models  are  developed  winch  allow  lift,  drag,  and  moment  estimation 
for  bodges  containing  inlets  normal  to  die  free  stream.  A  propulsion  model  for  lift¬ 
ing  fans  is  developed  relating  internal  and  external  floss. 

1.  ANALYSIS 

Tbe  basis  of  die  inlet  flow  field  analysis  is  the  assumption  that  the  flow  about 
an  inlet  normal  to  die  free  stream  mar  be  approximated  by  die  potential  flow  produced 
by  a  concentrated  sink  embedded  in  an  infinite  plane  above  which  passes  a  uniform 
flow.  The  force  produced  cn  .be  surface  containing  the  inlet  is  assumed  to  be  die 
force  acting  on  that  pom  zn  of  the  infinite  plane  obtained  by  projecting  the  planform  of 
die  surface  onto  the  infinite  plane.  For  purposes  of  calculation,  the  force  upon  the 
affected  portion  of  the  infinite  plane  is  divided  into  two  portions:  a  lip  force  acting  in 
the  immediate  vicinity  of  the  inlet,  and  a  surface  force  acting  on  *be  remaining  plan- 
form  area. 

a.  Lip  Forces 

Hie  lip  forces  and  moments  are  calculated  by  establishing  a  control  surface  of 
radius  Rj  about  the  inlet  centroid  and  applying  the  law  of  conservation  of  momentum 
(Figure  53). 

The  iip  force  in  the  free  stream  direction  is  expressed  as  shown  below. 
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The  pressure  and  velocity  at  the  hemispherical  control  surface  may  be  expressed  in 

terms  of  the  sink  induced  velocity,  u  . 

s 

PR,  =  PTa~  /I  (0®*2l4,  UsS +  Uj  ) 

U„  =  Us  *•  14  Sill/  cos  0 
u*  =  ^4>  +  us  s/n^cos0 
U2  -  -  Us  cos  ^ 

The  velocity  u  is  obtained  from  the  sink  strength  or  fan  flow  rate. 
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2 

LI  r  ~  W./4-TfR, 

v  I 

m  =  2  Af  l ^ 


The  conditions  at  the  lower  control  surface  are  shown  below.  The  constant  K  is 
defined  as  shown  and  the  constant  tj  is  a  dynamic  head  recovery  factor  for  the  inlet. 
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The  drag  is  then  found. 

0L  *  tff  14,  (  I-K) 

As  the  free  stream  velocity  component  is  reduced  to  zero  or  near  zero  in  passing 
through  the  fan,  we  may  calculate  the  total  drag  by  assuming  the  turning  to  take  place 
in  the  inlet. 
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The  moment  about  the  centroid,  and  normal  to  the  free  stream  is  found  as  shown 
below. 


X  r)u„dS  f  jP(nAr]clS 

Muy  Af  14U*  (R,/2  *  KA*) 

MLyAfUa,Of  ft,/ 2 

Mml_  =  Mu  coS/a 
Mn  =  '  Ml  sm/3 

As  in  the  drag  calculation,  the  parameter K  has  been  set  equal  to  zero  for  consistency. 
The  lift  may  be  calculated  from  the  expression 

Thus, 

L,=/|  Af [ofO-tyM?)  f  l42(rK")  ] 

The  parameter,  K,  is  retained  in  the  expression  as  the  change  flow  direction 
in  passage  through  the  fan  does  not  affect  the  pressure  at  the  lower  control 
surface,  the  source  of  the  term.  In  the  equations  for  drag  and  pitching  moment, 
the  parameter  K  reflects  an  unremoved  freestream  momentum  flux  and  later 
elimination  of  this  component,  in  the  fan  or  noizle,  will  cancel  the  terms 
containing  K. 
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b.  Surface  Forces 


The  forces  and  moments  on  the  remaining  effective  area  are  calculated  by 
integration  of  the  pressure  distribution  created  by  the  sink.  The  pressure  coefficient 
at  any  point  is  described  by  the  expression  shown  below. 


Of  _Af 

(4»  nr  * 


f  Of  Af  \ 
'  jrz  / 


The  lift  can  then  be  expressed  as  the  following  integral. 
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Integrating  with  respect  to  r  yields  the  intermediate  result  which  may  be  used  in 
conjunction  with  a  body  description  to  obtain  the  force. 
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In  a  similar  in?.mer  the  pitching  and  rolling  moments  are  found. 
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It  must  be  noted  that  the  effects  of  the  body  edges  on  the  pressure  distribution  created 
by  the  inlet  have  been  neglected,  and  that  the  forces  and  moments  due  to  edge  effects 
are  not  included  in  the  analysis. 

2.  LIFT-FAN  PROPULSION  MODEL 

A  simple  expression  has  been  obtained  relating  the  flow  parameters  to  propul¬ 
sion  parameters  for  lift-  fan  units.  The  model  is  based  upon  the  assumption  that  the 
pressure  difference  across  the  fan  may  be  represented  by  the  following  relation. 

AP*fyctutz 


If  the  exit  static  pressure  is  assumed  to  be  that  of  the  free  stream,  and  if 
allowance  is  nrnde  for  loss  of  dynamic  head  at  the  inlet  and  for  nonaxial  flow  at  the 
fan  entrance,  the  pressure  differential  may  be  related  to  flow  conditions. 

AP  *  k/>  [uf  *-(  .] 

The  fan  flow  rate  may  then  be  described  in  terms  of  free  stream  conditions  and  fan 
speed. 

.,2  2  /  Z\  ,Z 

Uf  -  Ctut  f  (l- K  J14 
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The  fan  flow  rate  is  often  expressed  in  terms  of  flow  ratio  and  tip  speed  ratio. 
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The  thrust  carried  by  the  fan  blades  may  be  obtained  in  terms  of  fan  flow  or  tip  speed 
parameters. 

Tf  =  AfAP 


ctutz 


The  effect  of  the  presence  of  a  fan  centerbody  may  be  considered  by  assuming  it  to 
produce  a  base  drag  force,  or  momentum  loss. 
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3.  COMPARISON  WITH  PREVIOUS  RESULTS 

If  the  preceding  relations  are  applied  to  ♦V  case  of  a  simple  fan  in  an  infinite 
plane,  having  an  inlet  which  turns  th*»  by  ninety  degrees  with  no  flow  losses,  the 
results  of  the  present  expressions  may  be  compared  with  previous  work  (20,  21). 

The  lip  and  surface  forces  become: 

LLrS  =  \Uf  fU<*  ] 

The  fan  thrust  reduces  to: 

Tf  =/f  Af  [U/-Uj-  ] 

And  the  total  lift  may  be  expressed  as: 

'  Ll>S  ' 

=  fAfV/' 

a  result  in  agreement  with  Theodorsen  (20) . 

The  expression  for  the  fan  flow  rate  is  then  found  to  be 


Z  Z 

Z  2  . 

Uf  t-  U&  for  Uj.  =  constant 
<4r° 


This  result  is  in  agreement  wi*h  Gruhame  (21),  who  derived  the  expression  by  other 
means. 
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4.  DEVELOPMENT  OF  HANDBOOK  METHOD 

The  simplicity  of  the  enpncd  model  outlined  abote  eaMe?  estimation  of  in¬ 
let  effects  without  resoit  to  sophisticated  computes  techniques.  The  lip  force  and  fan 
flow  equations  are  applied  in  the  form  given.  The  integrals  used  to  describe  the  sur¬ 
face  forces  may  be  easily  integrated  by  finite  summation  using  a  work  sheet  as  sk«n 
in  Table  1.  Because  the  induced  forces  are  concentrated  in  the  immediate  vicinity'  of 
the  inlet,  large  angular  increments  may  often  be  used  with  little  loss  of  accuracy. 
Semigraphical  integration  is  often  a  convenient  alternate. 
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TABLE  1 

INLET- SURFACE  FORCE  AND  MOMENT  CALCULATIONS 
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SECTION  VI 


NONLINEAR  BODY  AERODYNAMICS 

To  evaluate  the  body  contribution  to  the  aerodynamic  derivatives  of  V  STOI.  air¬ 
craft  it  is  necessary  to  inrlotle  the  range  of  flight  conditions  which  is  likely  to  be 
encountered  in  operation.  This  requires  the  aircraft  to  operate  at  extensive  angle  of 
attack  and  sideslip  since  the  flight  aititade  near  hovering  flight  is  not  determined  by 
aerodynamic  forces  but  rather  from  the  application  of  power.  There  is,  however,  an 
area  of  flight  where  the  aerodynamic  and  newer  effects  will  tend  to  be  of  the  same 
order  of  magnitude  and  where  neither  component  cap  be  safely  neglected. 

It  is  to  treat  the  aircraft  behavior  in  the  region  where  both  aerodynamic  forces 
and  power  are  im  jo  riant  in  supporting  the  aircraft  weight  that  nonlinear  aerodynamics 
is  important. 

To  obtain  the  total  aerodynamics  of  bodies  at  large  angles  of  attack  and  sideslip 
it  is  necessary  to  treat  both  the  potential  flow-  contribution  and  the  viscous  effects. 
This  study  has  utilized  slender  body  theory  ior  the  potential  contributions,  while  two 
methods  of  treating  the  viscccr*  ejects  iiave  been  examined.  The  first  method  using 
vortex  tracking  was  found  to  be  unreliable  within  the  present  state  of  the  art  and  to  be 
too  unwieldy  for  use  as  a  Handbook  procedure.  The  second  method,  using  the  viscous 
crossflow  concept  appears  to  be  quite  suited  for  a  handbook  procedure. 

t.  SLENDER  BODY  THEORY 

To  completely  represent  the  aerodynamic  forces  and  moments  on  a  body  it  is 

necessary  to  combine  both  the  potential  and  the  viscous  contributions.  To  obtain  me 

potential  contribution,  slender  body  theory  is  used.  The  development  of  the  equations 

(22) 

representing  the  slender  body  solutions  is  an  adaptation  of  the  methods  of  Sacks 
The  basic  equations  have  been  modified  to  permit  large  angles  of  attack  and  sideslip 
to  be  accounted  for.  The  ability  to  handle  arbitrary  body  shapes  has  been  attained  by 
introducing  «  method  of  obtaining  a  mapping  function  for  an  arbitrary  body  shape. 
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A  caesgszer  program  btti  nrina  atack  f**yel*«  five  coeagomessS.  force  and 
nooMaft  coeffiaeas  as  functions  a  .-esneltamB  aagse  of  an  act:,  rofii  angle,  yaw  veSeccty 
and  petchaqg  velocity.  Free  uteie  ot&aciats  the  aerodynamic  derivatives  can  be 
easily  obtained. 

Tfe  eigaataom  for  the  forifs  and  evm«U  pres  by  Reference  22  wiser:  racdLkesf 
to  eliminate  tie  roiliag  velocity,  the  effect  of  which  can  be  considered  negfigjisle  for  a 
body,  aad  the  time  derivatives,  can  be  am*s  is  complex  form: 
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The  evaluation  of  these  equations  requires  that  the  Comdex  potent"  al 


F-BMjLlZ+Y-  An [x)/z.r,-ej) (x) 

n*i  ' 


(82) 


be  known.  When  the  residue  of  the  potential  A((x)  is  known,  the  only  difficulty  in  eval¬ 
uating  the  total  forces  and  moments  arises  in  evaluating  the  first  term  of  the  moment 
equation  above: 
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rest  of  the  terms  are 
Sacks^^  has  shown 


readily  evaluated, 
that 
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Ii  rrmtmi  to eraiwe  the coefficients  ate  the  corffawl  PjO  is  act  required. 

The  cod&oca!5  at  be  evnlaased  by  re-expresssng  Eqpatjon  is  terms 
of  the  mainiiitg  circle  plane  coontete  Z,  i_e.. 
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Reference  22,  .Appendix  B,  shews  frat  this  expressive  can  be  era  hated  from  the 
coefficients  of  the  mapping  from  the  circle  plane. 

and  is  expressible  as 

F~rY.  *„/s”-r  r*jt  *T„(0  m 

II*/  * 

where 

<»ret0  (87) 

aod  V  *)  is  the  portion  of  the  potential  required  to  satisw  the  boundary  conditions  of 
the  body  and  is  equivalent  to  the  expression  of  Equation  (2i>)  of  Reference  23: 

/  Jit 

TNtt)-rc*U0{-f_o(a,n-2n  j£*n)£j 

(88) 
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This  expression  for  H  has  been  modified  from  the  definition  of  reference  22  to  utilize 
resultant  angle  of  attack  and  roil  angle  instead  of  angle  of  attack  and  sideslip. 

All  of  the  re'essary  information  is  now'  available  to  solve  Equations  (79)  and  (90) 
if  it  is  assumed  that  the  coefficients  of  the  mapping  function  (Equation  (85))  are  known. 
The  expression 

Rfepi) 

can  be  evaluated  by  integration  by  parts  to  be 

.  i(27r+u>J 

R$rj(il) = R [fzz\  E" 

-r[£?w] 


lis 


(93) 


(94) 


dace  the  first  term  of  the  iategrafioe  bjr  ports  Is  intagicary.  Os  dm  sarfice  of 
body  2  caa  be  represented  by  the  expression 


aad  the  expression  for  F  can  be  differentiated  with  respect  to  >  and  the  residue  of  the 
integral  can  be  obtained. 

Thus,  assuming  that  a  mapping  function  is  available,  it  is  possible  to  solve 
Equations  (79),  (SO)  and  (51)  for  Y,  L,  X,  M  and  L*  as  functions  of  o,  d,  q  and  r.  Also 
by  direct  differentiation  or  by  perturbing  either  o,  q,  q  or  r,  it  is  possible  to  obtain 
the  derivatives  of  the  coefficients  as  well.  The  above  analysis  can  also  be  readily 
extended  to  include  either  the  rolling  velocity  or  the  time-dependent  behavior  of  the 
coefficients  should  these  be  desired.  In  the  treatment  of  these  equations  it  has  been 
assumed  that  the  forward  velocity  can  be  resolved  into  an  axial  and  a  crossflow  com¬ 
ponent  which  permits  the  equations  to  be  solved  at  angles  of  attack  up  to  90  degrees. 

a.  Simplified  Method  of  Obtaining  Major  Mapping  Coefficients 

Section  III  has  described  the  method  which  has  been  developed  under  this  study 
for  determining  the  coefficients  of  the  mapping  function.  This  solution  makes  the 
problem  of  solving  for  the  slender  body  aerodynamics  solvable.  It  is  desirable,  how¬ 
ever,  to  si  mplify  the  procedure  for  treating  the  body  aerodynamics  by  eliminating  the 
need  for  mapping  the  body  since  this  requires  a  great  amount  of  work. 

To  obtain  a  method  of  using  slender  body  theory  for  a  handbook  method,  it  is 
necessary  to  find  an  approximate  method  for  estimating  the  most  important  coefficients 
of  the  mapping  function  in  a  simplified  fashion.  An  examination  of  Equations  (79)  and 
(80)  shows  that  except  for  the  term  Ax  all  the  variables  are  independent  of  the  mapping, 
i.e.,  they  can  be  obtained  directly  from  the  body  geometry  and  flight  attitude.  In  addi¬ 
tion,  it  can  generally  be  assumed  that  the  rolling  moment  of  the  body  can  be  considered 
negligible.  It  can  also  be  shown  that  the  derivatives  of  Y,  L,  N  and  M  with  respect  to  q 
and  r  are  dependent  only  on  the  radius  of  the  mapping  circle  and  the  coefficient  a,.  The 
other  coefficients  enter  into  the  derivatives  with  respect  to  angle  of  attack  and  sideslip 
(or  resultant  angle  of  attack)  through  the  cos  a  terms  appearing  in  Equation  (90).  It 
appears  from  this  expression  that  coefficients  other  than  r,,  aQ,  and  a,  may  be 
neglected  in  treating  the  body.  This  would  imply  the  simplification  of  '.'quation  (90)  to 
the  form 
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A,  (x)  <W-7TT  £ 


(95) 


This  approximation  and  a  suitable  method  of  estimating  r,,  and  a,  would  per¬ 
mit  a  relatively  simple  procedure  to  be  outlined  for  treating  the  potential  contribution 
to  the  body  aerodynamics.  The  most  suitable  method  of  approximating  these  three 
coefficients  is  to  treat  a  as  if  it  were  equivalent  to  the  centroid  of  the  cross  section, 
and  to  estimate  r  anti  a,  from  the  maximum  vertical  and  lateral  dimensions  of  the 
body.  If  2a  is  the  maximum  lateral  dimension  and  2b  is  the  maximum  vertical  dimen¬ 
sion,  then  by  letting 


_A*b 

•c  ~r 


(96) 


2 


(97) 


a  relatively  good  approximation  for  these  coefficients  is  obtained.  These  approxima¬ 
tions  have  been  made  for  the  T-38  body  and  their  values  are  compared  with  those 
obtained  from  an  analog  mapping  method  and  shown  in  Figures  54  through  56.  Consid¬ 
ering  the  complexity  of  the  chosen  body,  it  is  considered  that  the  agreement  obtained 
is  quite  good.  The  effects  on  the  coefficients  of  using  these  approximate  values  and 
neglecting  all  other  coefficients  are  shown  in  Figure  57. 

It  is  felt  that  this  is  an  adequate  demonstration  of  the  capaoilities  of  predicting 
this  portion  of  the  body  aerodynamics. 
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FIGURE  54.  EXACT  AND  APPROXIMATE  CALCULATION  OF 
MAPPING  CIRCLE  RADIUS 
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FIGURE  56.  EXACT  AND  APPROXIMATE  CALCULATION  OF 
MAPPING  COEFFICIENT 


2.  VORTEX  TRACKING 


One  method  of  obtaining  the  nonlinear  characteristics  of  bodies  is  to  treat  the 
viscous  flow  as  if  the  vorticity  were  collecting  in  a  single  vortex  pair.  Under  this 
assumption  the  paths  of  the  concentrated  vortices  are  traced  and  the  effects  of  these 
vortices  on  the  body  are  computed. 

This  computational  model  requires  several  simplifications  to  the  true  flow  pat¬ 
tern  about  a  body  and  the  degree  of  success  in  predicting  the  aerodynamics  is  depen¬ 
dent  upon  obtaining  a  satisfactory  set  of  assumptions.  Reference  23  presents  in  detail 

a  procedure  which  has  been  used  in  an  attempt  to  obtain  a  valid  model  for  arbitrary 

(24) 

bodies  at  angles  of  attack.  The  basic  model  used  is  that  of  Bryson  which  balances 
the  forces  between  the  concentrated  vortex  and  a  "feeding  sheet"  to  obtain  the  equation 

$  +  (98) 


which  relates  the  vortex  path  to  the  local  potential  and  to  the  rate  of  growth  of  the  vor- 
tex  strength.  This  equation  differs  from  that  for  a  free  vortex  (( ,  =  w,'  in  that  an  addi- 

r 

tional  velocity  ({,  -  fQ)  j: ,  is  imposed  on  the  vortex  due  to  the  presence  of  the  feeding 
sheet.. 

The  Bryson  model  also  postulates  the  existence  of  a  stagnation  line  on  the  body 
which  in  effect  permits  the  strength  of  the  vortex  to  be  computed  when  used  in  conjunc¬ 
tion  with  Equation  (98). 

The  method  of  Reference  23  is  a  generalization  of  the  Bryson  model  which  per¬ 
mits  the  introduction  of  an  arbitrary  body  with  unsymmetrical  flow  conditions  in  the 
crossflow  plane.  Thio  formulation  retains  the  major  features  of  the  Bryson  model 
but  generalizes  the  method  by  relating  the  flow  in  the  body  plane  back  to  a  similar 
flow  about  a  body  of  revolution  by  means  of  a  mapping  function.  Other  modifications 
are  made  as  required  to  permit  a  more  general  flow  pattern  and  to  permit  more 
flexibility  in  specifying  the  stagnation  line. 

The  method  was  used  on  the  T-38  body,  f  d  doing  so  showed  several  limitations 
in  the  model,  the  most  serious  of  which  was  caused  by  the  feeding  sheet  force.  This 
force  at  times  became  exceedingly  large  which,  in  turn,  produced  an  unrealistic  path 
for  the  vortex. 

During  the  present  study,  modifications  were  made  to  the  above  model  with  the 
intent  of  improving  the  accuracy  of  the  method.  The  mcst  satisfactory  modification 
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consisted  of  two  changes  in  the  model.  First,  the  effect  of  the  feeding  sheet  was 
neglected.  That  is,  the  vertices  were  assumed  to  move  freely  with  the  potential  flow 
field  including  the  effects  of  'he  other  vortex  df  the  pair  and  both  image  vortices  but 
assuming  no  feeding  sheet  force.  The  second  modification  was  to  relax  the  stagnation 
line  requirement.  Originally,  the  selection  of  a  pair  of  separation  lines  together  with 
the  requirement  that  the  cross  flow  velocity  went  to  zero  there,  permitted  the  specifi¬ 
cation  of  the  vortex  strengths.  This  condition  was  replaced  by  the  assumption  that  the 
rate  at  which  the  vortex  was  being  fed  was  equal  to  the  crossflow  velocity  at  the  sepa¬ 
ration  line.  This  condition,  which  may  be  stated  as 

%‘Uf>tUr  <89> 

is  more  realistic  as  a  condition  from  a  consideration  of  the  buildup  of  vorticity  in  the 
boundary  layer.  It  also  has  the  advantage  of  permitting  a  gradual  change  of  vortex 
strength  buildup  regardless  of  the  location  of  the  separation  line.  With  the  original 
assumption,  an  abrupt  change  in  the  separation  line  would  cause  an  abrupt  change  in 
the  vortex  strength. 

t 

Figure  58  shows  a  comparison  of  the  results  obtained  with  this  modified  model 
and  test  data.  The  agreement  at  angles  of  attack  up  to  40  degrees  is  quite  good  but  the 
theoretical  results  begin  to  act  erratically  at  higher  angles.  Further  attempts  to 
extend  the  results  to  higher  angles  of  attack  were  unsuccessful.  Additional  calcula¬ 
tions  were  made  with  this  model  in  sideslip.  In  sideslip  the  flow  pattern  in  the  cross- 
flow  plane  is  unsymmetrical  and  the  path  of  each  vortex  is  different.  Calculations 
made  under  these  circumstances  showed  a  very  erratic  behavior  and  did  not  agree  at 
all  well  with  available  test  data.  The  apparent  reason  for  this  was  the  unsymmetrical 
vortex  location  which  tended  to  exaggerate  the  asymmetry. 

This  model  has  several  drawbacks  due  to  the  simplifications  necessary  to  per-, 
mit  a  solution  to  be  obtained.  Among  these  are  separation,  a  quasi  two-dimensional 
flow  field  and  the  assumption  that  the  vorticity  behaves  as  a  rectilinear  vortex  with 
the  rotation  confined  to  the  crossflow  plane. 

It  was  decided  that  it  would  not  be  possible  to  refine  the  model  sufficiently  to 
allow  realistic  aerodynamic  loads  to  be  computed  within  the  program  time  limitations 
so  the  alternative  viscous  crossflow  model  was  studied  and  utilized  to  account  for  the 
viscous  effects  instead. 
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3.  VISCOUS  CROSSFLOW 


The  viscous  crossflow  concept  assumes  that  the  effects  of  an  angle  of  attack  or 
sideslip  can  be  estimated  by  treating  only  the  component  of  dynamic  pressure  in  the 
crossflow  plane,  and  that  the  forces  exerted  at  each  segment  of  the  body  are  indepen¬ 
dent  of  those  sections  forward  or  aft  of  this  segment.  This  treats  the  body  by  what  is 
essentially  a  strip  theory  based  on  local  two-dimensional  values  of  the  body  drag. 

Using  this  concept,  the  incremental  force  on  any  segment  can  be  written  in  the 

form 

f/vv'.  (100: 

where  (Cn  c  )  and  (Cn  c  )  are  local  values  of  the  crossflow  drag  areas  per  unit 

uy  y  uz  2 

length  in  the  yaw  and  pitch  directions  respectively. 

To  include  the  effects  of  pitching  and  yawing  velocity,  it  is  necessary  to  define 
-  the  crossflow  velocity  components  VQ  and  WQ  as 

V0  =  U^sina  Sin  0  -r(x-XC9) 

W0  =  “  U*  sin  a  s in  0  -  - XCg ) 

With  these  definitions  and  Equation  (100),  the  viscous  contributions  to  tho  forces 
and  moments  can  be  written 


(101) 

(102) 


t 
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i 

[n-im]  —  J  Jh-xc,)[cD  C  14 1  v0  \  -  i  c0  ce  W0 1 W,  I  Lx 
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(103) 


(104) 


The  above  formulations  permit  the  viscous  contributions  to  Y,  L,  N  and  M  to  be 
computed  as  functions  of  a,  <f>,  q  and  r.  No  attempt  has  been  made  to  include  a  viscous 
contribution  to  rolling  moment  since  it  does  not  seem  to  be  a  significant  contribution 
to  the  overall  aircraft  aerodynamics.  Nor  is  it  possible  to  formulate  a  model  which 
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would  include  acceleration  effects  of  the  body.  It  is  believed  that  the  above  formula¬ 
tions  give  the  most  'significant  effects  due  to  viscosity. 

The  values  for  the  crossflow  drag  coefficients  can  be  obtained  by  considering 
two-dimensional  section  data  for  the  shape  under  consideration  or  a  similar  shape. 
Calculations  made  thus  far  have  shown  that  treating  the  body  as  an  equivalent  ellipse 
will  give  reasonable  values  for  the  viscous  forces  and  moments. 

Calculations  were  made  for  the  T-38  body  using  for  the  crossflow  drag  coeffi¬ 
cient  Equation  (20)  of  Reference  25,  Chapter  ID. 


with 

Cj  = .0075 


(105) 


ag  and  bg  which  in  Equation  (105)  are  the  ellipse  semi-axes  parallel  and  perpendicular 
to  the  crossflow  respectively  were  assumed  to  be  the  maximum  dimensions  of  the 
local  T-38  cross  section. 

The  results  obtained  with  this  equation,  together  with  the  slender  body  values, 
are  shown  in  Figure  59  as  compared  with  test  data. 


SECTION  vn 


NONLINEAR  WING  AERODYNAMICS 

When  a  V/STOL  aircraft  is  operating  in  the  hover  or  transition  speed  regime, 
there  will  be  occasions  when  either  the  aircraft  angle  of  attack  or  sideslip  angle  are 
of  such  a  magnitude  as  t?  necessitate  an  investigation  of  the  nonlinear  aerodynamics 
of  the  vehicle.  This  section  deals  with  the  problem  of  determining  the  nonlinear  aero¬ 
dynamic  characteristics  of  a  wing. 

The  method  which  has  been  developed  is  a  modification  and  extension  of  the 

Weissinger  method.  Weissinger^26^  replaced  the  wring  by  a  concentrated  vortex  at 

the  quarter  chord  position,  and,  by  determining  the  downwash  due  to  this  vortex  and 

associated  trailing  vortex  system  at  a  number  of  control  points  on  the  three  quarter 

chord  line,  was  able  to  determine  the  strength  of  this  vortex.  To  do  this  he  balanced 

the  downwash  due  to  the  vorticity  with  the  upwash  due  to  the  attitude  of  the  wing  at  the 

control  points  on  the  wing.  This  procedure  is  effectively  one  of  taking  theoretical 

section  data  and  including  aspect  ratio  effects  by  accounting  for  the  trailing  vorticity. 

Choice  of  the  quarter  chord  line  for  the  vortex  implies  C  =0  and  satisfying  the 

ml/4 

downwash  boundary  condition  at  the  three  quarter  chord  line  results  in  a  limit  of  2  n 
per  radian  for  the  lift  curve  slope  for  rectangular  wings  as  the  aspect  ratio  is  increased 
indefinitely.  Or  equivalently  the  section  lift  curve  slope  is  2*r  per  radian. 

The  model  for  the  nonlinear  wing  problem  is  similar  in  that  nonlinear  section  data 
are  used  to  generate  information  for  finite  wings.  The  variation  of  the  section  lift  curve 
slope  throughout  the  angle-of-attack  range  is  accounted  for  by  locating  two  lifting  lines 
at  appropriate  locations  on  the  section  chord  and  the  relative  strengths  of  the  vortices 
are  determined  from  test  or  estimated  section  characteristics.  The  problem  of  pre¬ 
dicting  wing  lift  and  pitching  moment  characteristics  is  then  one  of  solving  for  the 
induced  downwash  angle  associated  with  the  planform  under  consideration  and  of  evaluat¬ 
ing  the  lift  and  moment  contributions  at  an  effective  angle  of  attack  as  determined  from 
the  induced  downwash. 
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1.  NONLINEAR  SECTION  MODEL 

The  analytical  model  used  to  represent  nonlinear  wing  section  data  assumes  that  at 
an  angle  cf  attack  of  90  degrees  the  normal  force  on  the  sect  ion  is  entirely  Aie  to  a 
viscous  drag,  arid  that  there  is  a  contribution  to  the  normal  force  at  an  arbitrary  angle 
cf  attack  c,  proportional  to  sin”a.  The  circulation  effect  of  the  section  is  represented 
by  two  vortices  as  illustrated  in  the  sketch  below.  The  boundary  condition  for  no  flow 
through  the  airfoil  surface  is  set  at  .  75  of  the  chord. 


t  Jisfying  this  boundary  condition,  we  obtain 


J7  +  JL  _  IL.-ioc  <106> 

»’rf, 

The  normal  foico.  contribution  due  to  circulation  effects  is  given  by 


r 

'~Nlr*D 


(C~L 


90 


.  A 

St,rt  x 


(107) 


in  which  (C.J  _  OA  is  the  normal  force  at  a  -  90°. 
n  a  —  yu 
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The  contraction  to  the  normal  force  coefficient  due  to  circulation  effects  is  assumed 
to  b?  represented  by 


(108) 


Equations  (106),  (107),  and  (108)  may  be  solved  for  Tj/Cc,  /*2  /Uc  from  given  section 
data  for  any  combination  of  and  As  is  the  case  of  the  linearized  model  of 

Weissinger^26^,  the  positions  of  the  two  vortices  are  found  by  considering  the  section 
pitching  moment  data. 


The  pitching  moment  abort  the  1/4  chord  position  C  is 


m 


i/4 


(109) 


in  which  x  is  distance  on  the  airfoil  chord  through  which  the  visoous  force  acts.  Sub¬ 
substituting  for  7^  /Uc ,  r  / Uc  from  Equations  (106),  (107),  and  (10 S)  we  obtain 

C«H  --  (c  +•  T  -■5)Ct,  JHD  sin* 

4-  (C„\  Sin *  (.25  -  Vc')  (U®> 

The  positions  of  the  vortices  represented  by  hj  and  in  Equation  (110)  are  chosen  to 
give  good  correlation  between  the  pitching  moment  calculated  from  Equation  (110)  and 
the  section  test  data.  The  strength  of  the  two  vortices  may  then  be  obtained  from 
Equation^  (106),  (107),  and  (108)  using  the  section  normal  force  data. 

This  procedure  has  been  checked  against  data  from  a  NACA  0012  airfoil (27\ 

With  the  two  vortices  positioned  so  that  Yi^c  -  .546  and  h2/4  =  .053,  calculations  of 
the  pitching  moment  against  angle  of  attack  are  compared  with  test  data  in  Figure  60. 
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FIGURE  60.  SECTION  MOMENT  COEFFICIENT  TEST  AND  CALCULATED  DATA 


The  strengths  of  ■*  two  vortices  r  and  r  which  will  give  the  required  section 

1  2 

normal  force  (given  in  Figure  61}  are  shown  in  Figure  62. 

The  capability  of  the  section  model  to  fit  the  normal  force  and  pitching  moment 
data  very  closely  suggests  that  this  section  model  may  be  used  to  calculate  wing 
characteristics . 


FIGURE  61.  SECTION  NORMAL  FORCE  COEFFICIENT  FOR  NACA  0012  AIRFOIL 


FIGURE  62.  CIRCULATION  STRENGTHS  OF  TWO  LIFTING  LINES  FOR  NACA  0012  AIRFOIL 


2.  APPLICATION  OF  THE  SECTION  MODEL  TO  WING  CALCULATIONS 


The  method  used  to  determine  finite  wing  characteristics  involves  determining 
the  effective  angle  of  attack  at  a  number  of  reference  stations  on  the  wing,  using  the 
section  modeL  The  weighting  of  the  circulation  at  these  sections  may  then  be  deduced, 
and  from  this  weighting  of  circulation,  the  total  circulation  around  the  sections  deter¬ 
mined  by  satisfying  the  no-flow  boundary  condition  at  appropriate  points  on  the  wing. 
Since  the  circulation  is  not  a  linear  function  of  angle  of  attack,  an  iterative  procedure  is 
necessary,  first  an  effective  angle-of-attack  distribution  is  assumed.  The  weighting 
of  the  circulation  between  the  two  vortices  is  then  determined.  The  circulation  across 
the  wing  span  may  then  be  calculated  and  a  new  effective  angle-of-attack  distribution 
determined,  and  this  procedure  is  repeated  until  convergence  of  effective  angle  of 
attack  is  achieved. 

a.  Down  wash  Due  to  the  Bound  Vorticity 

It  is  necessary  to  determine  the  down  wash  at  the  control  points  due  to  bound  and 
trailing  vorticity.  We  first  consider  the  downwash  due  to  the  bourn,  vorticity. 

Consider  a  bound  vortex  line  of  strength  Hr))  extending  from  (£i,  *?.)  to 
*7.+1).  Let  the  downwash  control  point  be  at  P(x,y,o, ).  This  arrangement  is  shown 
in  Figure  63. 


FIGURE  63.  GEOMETRY  FOR  BOUND  VORTICITY 
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Let  VT  be  the  velocity  Iniiced  by  the  bound  vorticity  frcm  ((.,  T) .)  to  (4.  ,  tj  )  at  P 

i-i  11  1+1  1+1 
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The  downwash  due  to  this  piece  of  the  lifting  line  w  follows  as 
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Equation  (114)  enables  the  downwash  to  be  calculated  in  terms  of  ^(rj)  for  a  jjfiven 
downwash  control  point  and  wing  geometry. 

b.  Downwash  Due  to  the  Trailing  Vorticity 

Let  V  be  the  velocity  induced  by  the  trailing  vortex  system  at  the  point  P(x,y,  z). 
Then  with  the  aid  of  Figure  64  we  obtain 
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The  lower  limit  in  the  integral,  LL,  represents  the  point  cm  the  bound  vortex  from  which 
the  vorticity  dT has  been  shed, 


wnen 


and 


T  =  (x-x,)l  +(£  j.  +(i  ri,){ 
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(116) 


(117) 


y.yj.1? 


FIGURE  64.  GEOMETRY  OF  TRAILING  VORliCITY 
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If  we  assume  that  the  trailing  vorticity  follows  the  effective  angle  of  attack  and  sideslip 
directions  we  obtain 


cLxt  =  Cos^t  cos  J Uxf  *  djf  +  At* 

dj,  =  —  sin- fit,  Jdx*  +  dtf  ir  dif 

cCt,  =  COS  fa  SL7i*t  Jdyf  +  d#*  f  dif 

Thus  equation  (117)  gives 

(it  =  £  l  -  J  ~f~  dx,  (118) 

Substituting  for  y^,  and  z^,  in  Equation  (116)  we  get 

r  =(x-x,)t+  (U9) 

Equations  (118)  and  (119)  may  be  used  to  deduce 
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in  which 


ft  =  ftx-’OVr*-?  -I-*9®*  +  1 ^  Xi) 
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From  Equation  (122)  we  may  obtain  RJin  the  form 
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in  which 
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in  Equation  (121)  may  be 


integrated  in  closed  form  and  an  integral  expression  deduced  for  w^.  The  analysis, 
given  in  Appendix  II,  results  in  the  following  expression  for  wt 
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Equation  (114)  with  Equation  (124)  are  the  basic  expressions  for  determining  the  down- 
wash  at  a  control  point  for  a  given  wing  geometry  and  flight  condition. 


3.  METHOD  FOR  CALCULATION  OF  WING  CHARACTERISTICS 


The  downwash  at  a  control  point  P(x,y,  z)  due  to  the  bound  voricity  In  the  interval 
(7^,,  17. +1)  together  with  the  associated  trailing  vorticity  is  WL  +  wt  which  is  given  by 
the  integrals  in  Equations  (114)  and  (124).  By  determining  similar  integral  expressions 
for  the  other  bound  and  trailing  vortex  segments,  it  is  possible  to  determine  a  series 
of  integrals  whose  sum  will  represent  the  downwash  at  a  general  control  point.  The 
problem  we  are  faced  with  is  the  indirect  one  of  being  given  the  downwash  at  the  con¬ 
trol  points,  and  being  required  to  determine  che  r(T))  distribution  which  will  give  this 
specified  downwash.  Since  r(TJ)  is  to  be  integrated  to  determine  the  downwash  we  have 
an  integral  equation  for  r(rj).  The  procedure  which  has  been  taken  to  solve  this  equation 

is  to  choose  a  number  of.  circulation  control  points,  say  rj . ,  and  to 

’"rite  r(T))  as  a  Lagrangian  polynomial  in  terms  of  the  circulation  at  these  control  points. 


ThttS  rii)  (i "lt-h . . (nJ  Ptu  (125) 

K) . (7.-7.)  (?«-?.) . (7.-1,.) 

This  representation  for  r(rj)  permits  the  integral  for  and  Wt  to  be  evaluated 

in  terms  of  the  coefficients  T(^),  i  =  0,  1, . n.  The  total  circulation  at  a  given 

wing  station  is  distributed  between  the  two  lifting  lines  by  weighting  the  vorticity.  This 
weighting  of  vorticity,  which  is  a  function  of  effective  angle  of  attack,  is  determined 
from  the  section  model.  An  example  of  the  weighting  function  calculated  for  the  NACA 
0012  airfoil  is  shown  in  Figure  65.  The  expression  for  the  downwash  may  be  written 
in  matrix  form  as  follows: 


AW  +  B(l  -w)]r  =  D 


(126) 


in  which  A,B  are  the  matrices  for  the  leading  and  aft  lifting  lines  representing  the 
coefficients  of  the  circulation  column  vector  I"?  Wand  I  -Ware  matrices  which  give 
the  weighting  of  the  vorticity  between  the  leading  and  aft  lifting  lines.  D  is  the 
downwash  column  vector. 

The  procedure  has  been  first  checked  for  the  linear  angle-of-attack  range. 
Calculations  of  the  spanwise  loading  for  an  aspect  ratio  5  rectangular  wing  with  sweep 
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back  angles  of  0°  and  45°  are  shown  in  figure  66.  These  calculations  compare 
favorably  with  the  calculation  of  Reference26and  test  data  of  Reference  28.  Further 
calculations  for  nonlinear  a  will  be  given  in  Volume  IL 
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FIGURE  66t  SPANWISE  LOADING  FOR  ASPECT  RATIO  5  RECTANGULAR  WING 


It  is  possible,  using  this  approach,  to  make  calculations  for  a  wing  with  a  flap. 
The  section  lift  and  pitching  moment  would  be  ch&nged  depending  on  flap  deflection 
angle,  flap  chord,  etc.  The  change  in  section  characteristics  would,  in  tun.  * T^nge 
the  circulation  weighting  function.  No  calculations  have  been  made  for  wings  with 
flaps.  It  may  be  necessary  to  include  two  more  lifting  lines  to  represent  the  flap 
for  more  satisfactory  correlation  with  test  data. 
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SECTION  Vffl 


EMPIRICAL  METHODS  USING  JET  FLOW  FIELD  THEORY 

The  computer  program  which  calculates  the  induced  velocity  field  due  to  single 
or  multiple  jets  exhausting  into  an  arbitrarily  directed  mainstream,  developed  in 
Section  II,  may  be  utilized  in  the  development  of  handbook  methods  for  estimating 
induced  forces  and  moments  on  supporting  structures  adjacent  to  the  exhausting  jets. 

In  this  section,  methods  for  the  wing  and  fuselage  are  presented  to  Illustrate  the 
utility  of  the  jet  flow  field  program  for  the  development  of  handbook  methods. 

1.  EMPIRICAL  METHODS  FOR  THE  WING 

Three  different  approaches,  each  utilizing  the  jet  flow  field  program,  have  been 
used  to  obtain  empirical  methods  for  the  wing.  The  first  approach  taken  was  to  substi¬ 
tute  an  equivalent  area  circle  for  a  wing  to  evaluate  the  induced  forces  and  moments. 
The  second  approach  considered  the  wing  planform  geometry  ( A©  method).  In  the 
third  method  an  equivalent  plain  wing  is  constructed  so  that  classical  methods  for 
calculating  wing  force  and  moment  may  be  used.  A  discussion  of  the  methods  and 
results  follows. 


a.  Equivalent  Area  Circle 


The  major  part  of  the  lift  loss  experienced  by  a  VTOL  vehicle  is  due  to  adverse 
pressure  distributions  on  the  lower  surface  of  the  vehicle  This  suggests 

that  a  good  estimate  of  theselift  losses  night  be  obtained  by  determining  the  lift 
losses  on  an  equivalent  circular  area  of  an  infinite  flat  plate  from  which  the  jet  or 
jets  are  exhausting. 

Pressure  coefficients  computed  by  the  Jet  Flow  Field  Program  were  integrated 
numerically  to  obtain  total  induced  forces  in  the  plane  of  the  jet  exit.  Figure  67  shows 
the  computed  variation  of  the  induced  force  on  a  flat  plate  around  a  normally  exhaust¬ 
ing  jet  with  velocity  ratio  for  two  values  of  the  parameter  ,  where  is  the  jet 

exit  area  and  S  the  portion  of  the  infinite  flat  plate  area  over  which  the  induced  force 
is  evaluated.  The  lift  increment  due  to  jet  operation,  AL,  is  defined  as  AL  =  F  +  T 
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FIGURE  68.  INDUCED  FORCES  ON  FLAT  PLATE  AND  RECTANGULAR  WING 
x  /c  =  0.5,  EQUIVALENT  AREA  CIRCLE 


where  F  is  the  jet-induced  force  and  T  the  static  thrust  of  the  jet.  Good  correlation 
with  the  experimental  data  of  Reference  26  is  discernible.  Experimental  data  on 
induced  forces  for  delta  wings  with  similar  values  indicate  the  utility  of  the 

jet  model  and  the  associated  computer  program  in  estimating  forces  induced  by 
exhausting  jets  on  finite  wings.  The  variation  of  induced  force  with  velocity  ratio 
for  a  range  of  the  parameter  is  shown  in  Figure  68.  Figure  69  shows  the 
induced  moments  computed  about  the  axis  of  the  jet  as  a  function  of  velocity  ratio  for 
a  range  of  the  parameter  ^5  . 

Experimental  data  points  obtained  for  a  rectangular  wing  with  an  aspect  ratio 
of  3.0  are  shown  in  Figure  68  for  comparison  with  computations  for  a  flat  plate  with 
the  same  value.  Very  good  agreement  for  the  low  velocity  ratio  data  points 

is  evident.  The  e>*ror  incurred  in  the  computed  values,  by  letting  a  circle  of  equiva¬ 
lent  area  represent  the  rectangular  wing,  becomes  significant  for  the  high  velocity 


ratio  data  point. 

This  substitution  of  an  equivalent  area  circle  for  the  purpose  of  evaluating 
induced  forces,  as  shown  schematically  in  Figure  70(a),  disregards  induced  forces 
on  the  shaded  portions  of  the  rectangular  wing  and  replaces  them  with  the  induced 
forces  on  the  cross-hatched  portions  of  the  circle,  For  the  rectangular  wing  of 


Reference  2,  this  affects  approximately  36  percent  of  the  total  area  of  the  wing.  By 
comparison,  Figure  70(b)  shows  the  relationship  between  a  delta  wing,  sweep  back 
angle  of  60°,  and  the  circle  of  equivalent  area.  The  shaded  area  which  is  replaced 
by  the  cross-hatched  area  now  represents  approximately  16  percent  of  the  total  wing 
area.  The  shaded  and  cross-hatched  areas  are  also  seen  to  constitute  a  better  sub¬ 


stitution,  in  terms  of  respective  distances  from  the  centrally  located  jet  exit,  for  the 
delta  wing  configuration  as  compared  to  the  rectangular  wing. 


FIGURE  70.  EQUIVALENT  CIRCULAR  AREAS  FOR  RECTANGULAR 
AND  DELTA  WINGS 
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b.  A9  Method 


Due  to  the  need  for  improved  accuracy  a  method  was  developed  which  includes 
the  wing  planform  in  greater  detail.  The  contribution  to  the 


FIGURE  71.,  SCHEMATIC  OF  SAMPLE  WING 


ratio  of  total  induced  force  to  installed  thrust  on  a  wing  of  the  shaded  segment  of 

Figure  71  can  be  evaluated  if  ^d0 

I Cp  (r/dj  d  (r/dj 

I/-2 

is  known  as  a  function  of  0  .  for  a  given  velocity  ratio  and  deflection  angle  of  the 
exhausting  jet.  Then  the  ratio  of  induced  force  to  thrust  ratio  on  the  segment, 


£  = 
T 


A 


7f 


ts  +0*3/ a 


f 

-  A&/2 


rm 


1/2. 


(r/do)  d  (r/«U  da 


(127) 


£.  = 
T 


R/d, 

A©  /  cP  (r/dj  d  (r/d 0) 

i/2 


(128) 


150 


By  employing  nume~ical  integration  of  the  ressure  coefficients  computed  by  the  jet 
flow  field  program  in  the  plane  of  the  jet  exit,  curves  o I 

ft/d!. 

fcF  (r/d.)  <i(r/do) 

i/2 

versus  R/c  at  constant  9  can  be  generated  for  a  range  of  velocity  ratios  and  et 
deflection  angles.  Figures  III—  1  through  Ifl-6  in  Appendix  HI  are  examples  of  sets  of 
curves  that  can  be  utilized  in  the  method. 

A  wing  of  arbitrary  plan  form  can  be  broken  down  into  a  number  oi  segments  as 
illustrated  in  Figure  72.  The  r2tio  of  induced  force  to  thrust  is 

(^,r2,a  oz) 

(#3.  R3.A*3> 

FIGURE  72.  SEGMENTATION  OF  SAMPLE  WING. 

R/c?, 

evaluated  by  obtaining  the  appropriate  value  of  J" Cp(r  a  ( *  / d») 

1/2 

from  the  cuives  and  multiplying  by  the  appropriate  AO  end  coefficient.  The 
cont-ibutions  arc  then  added  *o  retain  the  v>tal  mdecec  foiec  t'>  thrust  rati"  °n  the 

•A  "!j,. 

A  c  iTPans  ei  of  exp<  rimen'a!  and  calru  at  i  values  "f  surface  (■'  rcc  on  a 
fiat  riat*  :s  she  m  F:c;r»  7:-..  lie  resj‘t<-  ■■{  a  delta  wing  ard  a  re-tar.gular  wmg 
are  sir  at  jr  Figure  "4  a:  ",  Tab  • 


O  Experimental  data  for  U^/U^  =  8,  Ref  12. 

—  Calculated  values  by  jet  flow  field  program  for  U./U^  = 


FIGURE  73.  DISTRIBUTION  OF  SURFACE  FORCE  AROUND  THE  JET 
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FIGURE  74.  INDUCED  FORCES  CN  A  DELTA  WING 
A  ft  METHOD 


c.  Equivalent  Plain  Wing 


i 


A  jet  exhausting  into  a  freestream  will  induce  a  downwash  component  of  velocity 
over  the  wing  planform.  The  surface  boundary  condition  will  be  satisfied,  in  a 
linearized  form,  when  this  downwash  velocity  component,  is  cancelled  by  some  mecha¬ 
nism.  A  convenient  way  of  producing  a  downwash  to  cancel  that  produced  by  the  jet 
is  to  construct  an  equivalent  plain  wing  with  an  appropriate  camber  distribution. 

I 

,  Classical  methods  may  then  be  used  to  calculate  the  power  induced  force  and  moment 
on  the  wing. 

Power  induced  section  characteristics  may  be  determined  using  linearized 

i 

thin  airfoil  theory  as  follows. 

,  Consider  the  problem  of  determining  the  lift  and  pitching  moment  on  a  wing 

section  in  the  presence  of  a  jet  exhausting  into  a  uniform  freestream  U  directed  along 

the  positive  x-axis.  Let  the  jet  induced  downwash  over  the  wing  section,  0  4  Y  4  C. 

be  w(x).  Let  z  be  the  displacement  normal  to  the  x-axis  of  the  equivalent  airfoil 

section  as  defined  above.  Then,  ,  and  so 

dx  Uoo 

x 

2  (x)  =  fi**  (129) 

J-  UcO 

;Thus,  writing  ^ 

dl.  =  H  =  X-  cos n.<t>  (130) 

dx  Uo  0 


in  which  p  is  given  by  x  =  c/2  (1  +  cos  p  )  , 

•s' 

we  obtain 
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R  =  —  I  cosn.&  cl(t  for  a  i 

n  T  J  ’ 


The  section  lift  and  pitching  moment  coefficients  may  be  determined  as  functions 
of  Eq,  Bj,  and  B2  (details  may  be  found  in  Reference  47) .  Thus, 


i  S,) 

and 


(1**3) 


Cm--  It  (S,  .  3,)  -  4  Ct  (134, 

Finite  wing  lift  and  pitching  moment  coefficients  may  be  determined  using  .he 
methods  described  in  References  48  and  49  . 

Examples  of  equivalent  section  camber  lines  are  shown  in  Figures  75,  70  and  77. 
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ER  CALCULATION 


FIGURE  77  .  CURVES  FOR  EFFECTIVE  ANGLE  OF  ATTACK  AND  CAMBER  CALCULATION 


d.  Application  to  Calculations  of  Downwash  at  Tall  Surface 


The  jet  flow  field  program  may  be  used  to  compute  the  velocity  components  in¬ 
duced  by  exhausting  jets  at  arbitrary  points  in  space.  Thus,  the  program  has  the 
capability  of  determining  the  jet  induced  downwash  in  the  vicinity  of  the  horizontal 
tail.  The  power  induced  contribution  to  the  pitching  moment  due  to  the  presence  of 
the  horizontal  tail  may  then  be  deduced. 

The  results  of  a  typical  calculation  are  shown  in  Figure  78.  The  jet  induced 
downwash,  €  ,  is  shown  for  two  tail  stations. 


y/d0 

13  Jet  Diameters  Aft,  S’.  =  90°,  U./U„  =  5 

J  J  “ 

FIGURE  78  .  JET  EXHAUST  EFFECTS  ON  TAILPLANE 


2.  BODY  METHOD 

a.  Application  of  Lagally’s  Theorem  in  Relation  to  the  Present  Problem. 

The  most  fundamental  information  in  the  present  study  is  the  knowledge  of  the 
induced  velocity  field  due  to  the  powered  lift  system.  Based  on  this,  alternative 
procedures  may  be  found  to  predict  the  power  induced  forces  and  moments  on  an 
immersed  body.  Such  alternative  methods  provide  not  only  competing  approaches 
but  often  an  understanding  of  the  problem  derived  from  taking  a  different  viewpoint. 

The  end  result  of  the  present  jet  model  is  a  distribution  of  sinks  and  doublets 
along  the  jet  path.  A  source,  a  sink,  or  a  doublet  will  all  induce  a  force  to  attract 
the  submerged  body  towards  the  singularity  on  the  basis  of  incompressible  inviscid 
flow  theory.  This  explains  in  part  the  loss  in  lift  of  a  V/STOL  aircraft  in  the  presence 
of  a  downward  jet.  If  instead  of  a  single  source,  a  system  of  sources  (sinks  are  con¬ 
sidered  to  be  negative  sources)  and  doublets  is  introduced,  the  total  force  and  moment 
are  the  additive  forces  and  moments  induced  by  individual  singularities.  This  is,  in 
effect,  Lagally's  theorem  and  holds  true  in  two-dimensional  as  well  as  in  three- 
dimensional  cases  (References  30  and  31).  It  applies,  in  principle,  to  a  submerged 
body  of  arbitrary  shape  in  a  uniform  stream  with  arbitrary  number  of  sources,  dotblets, 
and  higher  order  singularities.  Using  this  theorem  for  the  present  problem,  the  task 
is  ther  to  determine  the  velocities  and  their  gradients  "induced"  by  the  fuselage  and  the 
image  systems  for  singularities  outside  the  immersed  body,  since  these  velocities  and 
gradients  represent  the  force  and  moment  exerted  on  the  fuselage  by  the  jet. 

For  a  two-dimensional  body,  finding  such  an  image  system  is  possible,  as  long 
as  it  can  be  mapped  into  a  unit  circle  on  the  transformed  plane.  Thus,  the  power-effect 
forces  and  moments  may  be  calculated  for  a  general  class  of  bodies.  Unfortunately, 
there  is  no  corresponding  procedure  for  a  three-dimensional  body.  The  image  ox  a 
singularity  in  that  boundary  is  usually  a  complex  one.  Only  for  the  special  case  of  a 
sphere,  may  the  image  system  be  found  conveniently  and  exactly.  Since  a  fuselage  is 
a  three-dimensional  body  of  complicated  geometry,  it  is  virtually  a  hopeless  task  to 
determine  the  image  system  for  this  kind  of  body.  However,  if  we  can  approximate  the 
fuselage  by  a  prolate  spheroid  of  equal  volume,  the  possibilities  are  improved.  On 
the  strength  of  this  simplification,  we  propose  to  employ  Lagally's’ theorem  as  a 
possible  empirical  method  for  the  fuselage. 

Before  we  dwell  on  the  main  aspect,  some  explanation  seems  to  be  pertinent  in 
regard  to  th.  justification  of  this  simplification.  Lagally's  theorem  is  derived  by 


means  of  the  integral  equations  in  which  the  force  and  moment  are  expressed  in  terms 
of  momentum  flux  through  an  enclosing  control  surface.  The  process  of  integrating  the 
detailed  surface  pressure  has  in  effect  been  circumvented.  The  details  of  the  fuselage 
geometry,  which  reveals  in  the  form  of  velocity  potential,  become  increasingly  indistinct 
as  the  distance  from  the  surface  increases.  For  a  two-dimensional  body,  the  general 
complex  potential  at  large  distances  can  be  expressed  in  terms  of  a  Laurent  series. 

The  leading  term  is  proportional  to  l/(x  +  iy),  which  behaves  like  a  doublet.  The 
higher  order  terms  signify  more  refined  details.  By  the  same  token,  the  disturbance 
potential  of  a  three-dimensional  body  decays  like  l/(x  +  y  +  z  )  at  great  distances. 

This  is  recognized  to  be  a  doublet.  Thus,  the  total  force  and  moment  on  an  immersed 
body  by  all  singularities  will  be  insensitive  to  the  detailed  geometry,  provided  that  we 
can  show  that  a  substantial  portion  of  the  forces  and  moments  are  derived  from  the 
"induced"  velocities  and  their  gradients  at  positions  of  singularities  sufficiently  far 
from  the  fuselage.  By  examining  the  outputs  from  the  jet  model  in  which  the  strength 
of  singularities  increases  with  the  distance  away  from  the  exit  and  using  a  simple  body, 
say  a  sphere,  one  could  confirm  this  conjecture. 

b.  Lagally's  Theorem  and  Computation  of  a  Sphere  as  a  Test  Case 

Since  the  computation  is  based  on  inviscid  flow,  the  forces  acting  on  the  surface 
of  a  fuselage  arise  from  the  pressure  alone.  Thus,  the  components  of  the  force  ax-e 
given  by 

fi  =~  Pn:dS  ,  i=  / ,2,3  (135) 

(5) 

where  p  is  the  pressure,  and  n.  are  the  cartesian  components  of  a  unit  normally  drawn 
outward  to  the  element  dS  of  the  fuselage  surface  S.  This  equation  may  be  reduced 
by  Bernoulli's  equation  and  Green's  theorem  to 

Hi  r 

F;  =  2  f  T“jajn0dS,  (136) 

*=/  “(Sa) 

where  the  repeated  indices  j  refer  to  summation,  a.  are  the  velocity  components,  and 
Sa  are  the  surfaces  of  spheres  of  infinitesimal  radii  which  surrounds  singula  es 
at  the  positions  ra  (a  =  1,2,  —  If,).  This  ic  a  general  expression.  To  specify  tne 
nature  of  the  singularities,  we  shall  make  distinction  between  a  source  and  a  doublet. 
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For  a  source  of  strength  m  located  at  the  position  r  ,  the  velocity  components 

cL  ft 

on  the  surface  of  the  infinitesimal  sphere  a,  after  the  series  expansion,  may  be  written 
as 


“«=  % 


id 


Pic.  (Xj  -XL<1) 

Y 3 


+ 


where  the  subscript  a  indicates  the  position  of  singularity,  and  q.&  are  the  velocity 
components  at  the  point  r„  induced  by  all  causes  except  the  singularity  there.  The 

ft 

second  term  represents  the  velocity  components  by  the  singularity,  in  which  r  refers 
to  the  radius  of  the  infinitesimal  sphere.  The  third,  fourth  and  the  rest  of  the  terms  are 
members  of  the  Taylor  series.  Upon  substitution  of  this  expression  into  Equation  (136), 
we  obtain 

R<l=  '  (137) 


This  may  be  looked  upon  as  a  force  F„  acting  on  a  fuselage  in  a  line  through  the  source 

ft 

at  r  .  Thus,  the  moment  on  the  fuselage  due  to  in  is  the  vector  product  of  ro  and 
u  a.  ci  ci 

Fo  and  is  given  by 

a 

Mia.=  *‘2,3  (138) 


where  is  the  permutation  symbol,  «123  =  «231  =  «312  =  h  «213  »  e132=e321  =  311(1 

=  0,  if  i,  j,  k  are  not  all  different.  For  a  doublet  of  strength  ^  situated  at  the  point 
r^,  the  force  and  moment  may  be  obtained  by  a  similar  but  somewhat  more  complex 
procedure.  For  details,  one  is  referred  to  the  above  mentioned  references. 

The  total  force  and  moment  on  the  fuselage  are  the  summation  of  individual 
contributions  and  are  given  by 


Nt  tJz 

Fc  =  w[  l  +  2  1  -  ‘-h1-2-3 

a=\ !  b=i 


«;=  We,n[ 2 +2  =  3 

a=i  *--*  *  J 


-^2 
6=/ 


(139) 
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The  subscripts  a  and  b  have  been  dropped  for  convenience.  The  quantities  m,  Jl ,  r 
are  known  from  the  jet  flow  field  program,  but  the  velocity  components  q.  must  be 
determined  according  to  the  given  configuration.  Thus,  to  use  Lagally's  theorem,  the 
problem  is  that  of  determining  these  components. 

As  stated  above,  the  quantities  q.  at  the  point  in  question  are  the  velocity  com¬ 
ponents  induced  by  all  causes  except  the  singularity  there.  But  the  contributions  due 
to  other  singularities  outside  the  body  may  be  neglected  without  affecting  the  final 
result,  because  the  force  and  moment  induced  at  one  singularity  by  another  are  equal 
and  opposite  to  those  at  the  latter  singularity  by  the  former.  This  is  sometimes 
referred  to  as  the  reciprocal  relation.  Consequently,  the  relevant  induced  velocities 
q.  resulting  from  a  body  in  a  uniform  erf  ssflow  with  an  exhausting  jet  may  be  classified 
as:  (a)  the  velocity  induced  by  the  body,  (b)  the  velocity  induced  by  the  image  system 
of  sources,  and  (c)  the  velocity  induced  by  the  image  system  of  doublets.  These  veloc¬ 
ities  are  all  functions  of  body  geometry.  A  general  expression  for  an  arbitrary  body 
is  not  known  at  present.  The  velocities  of  type  (a)  can  be  determined,  if  we  assume 
that  a  fuselage  be  replaced  by  an  ellipsoid  of  equal  volume.  The  solution  to  this  is 
known  and  can  be  found  in  many  books.  Unfortunately,  it  is  still  difficult  to  obtain 
velocities  of  types  (b)  and  (c)  even  under  this  assumption.  Thus,  it  appears  that  some 
further  approximation  is  needed.  Before  exploring  possible  approaches,  it  seems 
worthwhile  to  know  the  magnitude  of  these  velocities  and  the  relative  importance  of 
each  type.  To  do  this,  we  shall  work  out  a  simple  example. 

Hie  only  known  three-dimensional  configuration  for  which  all  three  types  of  the 
induced  velocities  can  be  determined  exactly  is  a  sphere.  For  this  reason,  we  choose 
this  as  our  test  case.  The  power  effect  may  be  calculated  with  the  singularity  distribu¬ 
tion  given  by  the  jet  flow  field.  It  turns  out  that  type  (a)  contributes  roughly  80  percent 
of  the  total  effect.  For  an  elongated  body  such  as  a  fuselage,  it  is  conceivable  that 
this  percentage  may  be  somewhat  reduced.  Nevertheless,  the  evidence  seems  to  point 
out  the  fact  that  the  effect  of  the  image  system  is  relatively  small  as  compared  to  that 
due  to  the  body  and  may  be  neglected  for  the  first  approximation. 

Prior  to  this  finding,  attempts  were  made  by  the  trial-and-error  methods  to  get 
approximate  image  systems  for  a  source  outside  a  prolate  spheroid,  using  the  image 
system  of  an  ellipse  and  the  correspondence  between  a  circle  and  a  sphere  as  the 
guide.  This  effort  did  not  persist  long  enough  to  examine  various  approaches,  out  it 
appears  that  some  useful  approximation  should  result,  depending  on  the  accuracy 
required. 
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(c)  Application  cf  Method  to  Fuselage 


By  choosing  elliptic  coordinates  v,  £,  cj,  the  velocity  potential  at  a  field  point 
for  flows  past  a  prolate  spheroid  at  an  angle  of  attack  a  is  given  by 


(p  =-  U^Cos# -kv^  -  Umsinot4( c°$u>  -Av[jih  jzj  -  l] 


1 

?r  I 


- A'(i-y2f(t2-'iJ[j^f^T  "1^7 3to5tu, 


A  = 


Umkco  S0< 


fc  I  fl  fo+  I 
#*-/  “2^  &-I 


A<_  ~UC0'k 

/i  •“  i  x  J-j 


2  A- 1 


=  , 


e  = 


where  U^is  the  freestream  velocity,  a  and  c  are  the  polar  and  equatorial  radii  of  the 
given  spheroid  that  is  a  member  of  the  confocal  family  £  =  The  relationship  between 
this  and  the  cartesian  coordinate  system  is 


-  Vf/im  lu  ,  $  =  v?  cos  u>  . 

j.  £ 

I*  =  Mi-  pzf($2->) 


in  which  the  surfaces  $  =  constant  and  v  -  constant  are  confocal  ellipsoids  and  hyper¬ 
boloids  with  the  common  foci  of  (±k,  0,  0)  as  illustrated  in  Figure  79.  By  knowing 
the  velocity  potential  and  this  coordinate  transformation,  we  can  determine  q.  and 
dq./ax.  at  every  position  of  the  sources  and  doublets  in  the  jet  flow  field.  Application 
of  Equation  (139)  then  gives  the  force  and  moment  on  the  given  body  caused  by  the 
induced  velocities  of  type  (a).  These  are  assumed  to  equal  approximately  the  force 
and  moment  on  the  fuselage.  Calculations  were  subsequently  made  for  the  test  fuselage 
with  a  lift  jet  at  freestream  to  jet  velocity  ratios  equal  to  0. 2  and  0. 3.  The  results  are 
shown  in  Figure  80,  along  with  the  test  data  and  those  computer!  by  the  transformation 
method. 
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FIGURE  80.  POWER-EFFECT  LIFT  FOR  LIFT  JET  WITH  FUSELAGE 
CALCULATED  BY  LAGALLY’S  THEOREM  AND  ITS  COMPARISON 
WITH  TRANSFORMATION  METHOD  AND  TEST  DATA 


SECTION  IX 


HOVER  AERODYNAMICS 

The  hover  region  of  V/STOL  flight  requires  special  treatment  due  to  the  lack  of 
steady  freestream  conditions.  Two  areas  of  special  treatment  for  this  flight  regime 
have  been  studied.  The  first  deals  with  the  powei  effects  in  a  state  of  true  hover. 

The  second  treats  nonlinear  and  nonsteady  power  off  aerodynamics  at  near  hover  con¬ 
dition.  These  studies  provide  a  limit  point  for  the  transition  aerodynamics  and  permit 
an  approach  to  V/STOL  hover  flight  analysis. 
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1.  POWER  EFFECTS  IN  HOVER 


A  number  of  experimental,  empirical.,  and  theoretical  investigations  have  been 
devoted  to  the  problem  of  determining  the  forces  exerted  on  a  stationary  body  from 
which  a  jet  exhausts.  Particular  attention  has  been  devoted  to  thrust  losses  of  lifting 
jet  configurations  both  in  and  out  of  ground  effect;  a  good  summary  of  these  efforts 
may  be  found  in  Reference  (32). 

The  most  pertinent  results  to  date  are  contained  in  an  empirical  relation  for  lift 
loss  out  of  ground  effect, 


LT 

from  Reference  (33),  or 


X;  /  d 


(140) 
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from  Reference  (34),  and  in  an  empirical  correlation  for  the  change  in  force  due  to  the 
presence  of  the  ground  plane, 
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(35) 

,RAE  data 


(142) 


..  2.  OZ 


,NACA  data 


(36) 


(143) 


The  preceding  results  are  primarily  applicable  to  single  or  closely  spaced,  cen¬ 
trally  located  jets.  Although  derived  from  small-scale,  cold-jet  tests,  Equations  (140) 

(34  37) 

and  (141)  have  been  shown  to  apply  to  more  complicated  configurations.'  *  Equa¬ 
tion  (142)  has  been  apnlied  with  reasonable  success  to  full-scale  configurations  having 
closely  spaced  turbojet  propulsion. (37»38) 

Equations  (140)  and  (141)  differ  because  of  ivie  addition  of  a  jet  pressure  parame¬ 
ter  by  a  later  investigator. 
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Equations  (142>  and  (143)  vary  because  different  data  were  used  »n  determining 
the  constants.  A  cursory  examination  of  later  NASA  data^^  indicate!  a  preference 
for  Equation  (143)  and  a  brief  effort  was  made  to  determine  if  the  enr  :rical  constants 
of  Equations  (142)  and  (143)  could  be  related  to  the  area  ratio,  jet  dei  -y  rate,  or  out- 
of-ground  effect  loss. 


The  nature  of  Equations  (140)  thru  (143)  indicates  that  know’. e  of  jet  decay 
rate,  out-of-ground  effect  loss,  and  true  jet  thrust  for  simple  configurations  is 
required  for  such  an  effort.  A  requirement  for  this  data  group  gc  ^atly  reduces  the 
amount  of  useful  data  available.  As  can  be  seen  in  the  data  sum  i  ary  of  Table  3,  only 
data  from  References  (33)  and  (39)  remain,  and  Reference  (39)  doos  not  consider  ground 
effect.  References  (34)  and  (37)  contain  the  proper  data  group  but  involve  more  complex 
body  geometry. 

During  early  attempts  to  correlate  Equations  (142)  and  (143) ,  nonlinearities  of 
log-log  plots  of  G/T  -  (G/T)^  versus  h/d  led  to  the  plotting  of  G/T  versus  h/d  in  semi¬ 
log  form.  One  such  plot  is  shown  in  Figure  81.  As  the  quality  of  the  correlation 
seemed  remarkable,  additional  data  was  correlated  in  this  manner.  Figure  82  indi¬ 
cates  that  correlation  is  possible  for  circular  flat,  plates  with  an  error  band  of  20  per¬ 
cent.  Datum  outside  this  band  and  much  of  the  scatter  within  this  band  may  be 
attributed  to  errors  in  transcribing  information  from  Reference  (33), 


An  expression  of  the  form 


(144) 


is  therefore  suggested.  The  value  of  B  obtained  from  the  data  of  Figure  81  differed 
from  that  obtained  from  Figure  82;  B  was  found  to  be  0.97  and  1.23,  respectively. 
The  difference  may  be  attributable  to  the  different  jet  turbulence  levels  produced  by 
a  rectangular  plenum  (Figure  81)  and  a  cylindrical  plenum  (Figure  82).  For  the  two 
configurations  that  have  been  examined,  it  is  noted  that  the  constant  B  is  inversely 
proportional  to  the  square  root  of  the  out-of-ground  effect  loss. 

e'_ 

6" 


rc/n. 
L  ( G/  T )' 


O-  5 


(145) 
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TABLE  3  Jet  Hover  Data  Summary 


circular  rlate  R  -  rectangular  or  square  plate 

pod  or  fuselage  S  -  single  jet 

lift  loss  T  -  jet  thrust;  triangular  plate 

limited  data,  lift  T  -  calculated  jet  momentum  thrust 


OIO 


h/d 


FIGURE  81.  LIFT  LOSS  FOR  CIRCULAR  FLAT  PLATES,  RECTANGULAR  PLENUM  (Ref.  33) 
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Equation  (145)  in  combination  with  relations  (140)  and  (141)  may  be  used  to  relate 


the  constant  B  to  the  rate  of  jet  (impact)  pressure  decay. 

B_  _  p  )/  /(  b  (tjt/PTyp-  p)  /  'j 

6  M  >./d  )  A/b /  \  a(x/<J)  / x/J 


0.25 

(146) 


Relation  (146)  in  combination  with  Equation  (144)  suggests  Equation  (147). 


G 

T 


b  ( tji/frjp'  p)  /  b  U/d) 
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J 


(147)* 


However,  further  analysis  is  necessary  before  conclusions  can  be  drawn  regard¬ 
ing  the  validity  of  Equation  (147)  which  has  been  developed  using  limited  data  for  cen¬ 
trally  located  jets  in  circular  plates. 

It  should  be  noted  that  a  correlation  similar  to  that  of  Equation  (144)  may  be 
found  in  the  form 


0'"c(  g/t)  /(h/<in 

flu 


(148) 


This  suggests  that  (G/T)  is  proportional  to  /s/A.  - 1  and  not  to  \/s/Ai  as  indi- 
cated  in  Equations  (140)  and  (141).  This  seems  intuitively  correct  as  one  would  desire 
the  calculated  load  on  the  plate  to  go  to  zero  if  no  plate  is  pr.  _>ent,  S/A.  -  1.  A  com¬ 
parison  is  shown  in  Figure  83.  Thus,  \/s/A.  - 1  is  a  significant  parameter  both  in  and 
out  of  ground  effect. 

It  is  probable  that  relations  of  the  type  listed  above  will  not  be  useful  for  a  body 
very  near  to  the  surface  where  small  details  of  specific  body  and  jet  geometry  become 
significant.  Additional  terms  will  be  required  in  this  region. 


kNote:  For  circular  planforms,  / S/A^  -  D/d. 
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f.  PERTURBATION  METHOD 

To  attack  the  problem  of  nonlinear  and  nonsteady  aerodynamics  at  near  hover 
condition,  a  study  combining  theoretical  analysis  and  water  tank  tests  was  undertaken 
to  determine  the  feasibility  of  obtaining  a  method  for  predicting  stability  derivatives 
under  flight  conditions. 

A  small  existing  water  tank  was  used  to  try  out  possible  techniques  of  obtaining 
test  data  for  further  study.  This  tank  and  the  arrangements  used  for  these  tests  are 
shown  in  Figure  84.  The  model  was  suspended  from  a  beam  that  was  free  to  roll 
across  the  tank.  Motive  power  was  supplied  by  unbalancing  the  weights  on  the  two 
weight  pans.  Weights  were  shifted  from  the  aft  side  to  the  forward  side,  providing 
an  incremental  loading  while  keeping  the  total  mass  of  the  system  constant.  An  accel¬ 
erometer  attached  to  the  beam  supplied  a  force  trace  versus  time  plot.  Strain  gages 
were  mounted  on  the  supporting  rod  to  provide  moment  measurements. 

Acceleration  and  moment  plots  were  obtained  for  the  model,  and  measurements 
were  made  both  with  the  model  in  water  and  in  air.  Typical  plots  are  shown  in  Figure 
85  for  different  incremental  weight  loadings.  The  initial  accelerations  in  both  cases 
were  plotted  against  incremental  weight,  and  straight  lines  were  faired  through  the 
points.  From  these  straight  lines,  it  was  possible  to  determine  the  apparent  mass  of 
the  mode'  ?n  water.  (The  apparent  mass  in  air  was  assumed  to  be  negligible  compared 
to  that  in  water.) 

The  acceleration  curves  were  integrated,  and  the  values  of  force  and  moment 
versus  velocity  th  m  were  obtained  for  three  different  incremental  weights.  The  force 
due  to  the  acceleration  (apparent  mass  multiplied  by  acceleration)  was  subtracted  from 
the  total  force,  and  a  net  force  was  obtained  which  was  assumed  to  be  due  to  viscous 
effects.  The  results  of  this  data  reduction,  reduced  to  drag  coefficient  and  plotted 
against  velocity  multiplied  by  time,  are  shown  in  Figure  86.  The  shape  of  the  three 
incremental  weight  curves  shows  good  correlation  throughout  the  test  range.  The 
four-pound  curve  is  somewhat  low,  but  this  could  be  attributed  to  inaccuracies  in  the 
data  reduction. 

Since  both  accelerometers  and  strain  gages  were  not  only  measuring  the  linear 
motion  of  the  rigid  body  but  also  were  picking  up  elastic  transients  and  other  spurious 
inputs,  isolation  of  the  reduced  data  was  somewhat  difficult.  This  problem  was  fur¬ 
ther  increased  by  the  use  of  strip  chart  recorders  which  contributed  their  own  noise 
and  pen-a.'m  inertia  to  the  riots.  Of  the  models  tested  —  rectangular  plate,  cylindrical 
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In  Air 


In  Water 


4-lb  Load 


5-lb  Load 


FIGURE  85.  TYPICAL  TIME  HISTORIES 
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FIGURE  86.  VISCOUS  DRAG  COEFFICIENT  FOR  TRAPEZOID  PLATE 
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tube  (see  Figure  84),  and  trapezcidai  flat  plate  —  the  latter  offered  the  most  clearly 
defined  data,  especially  for  calculation  of  apparent  mass. 

The  recorded  data  were  surprising  in  two  more  ways:  First,  the  drag  coeffi¬ 
cients  obtained  near  the  end  of  the  plots  were  about  three  (based  on  the  frontal  area), 
which  seemed  to  be  too  high.  Second,  the  drag  coefficients  in  the  beginning  stages  of 
the  run  were  high,  with  indications  that  the  drag  was  more  nearly  proportional  to  V 
than  to  V^. 

The  data  were  analyzed  according  to  the  concept  of  indicial  functions.  Northrop 
assumed  that  the  curve  for  the  particular  derivative  or  coefficient  under  consideration 
could  be  built  up  from  a  series  of  step  functions.  For  example,  from  the  tests  made 
so  far  it  can  be  assumed  that  the  derivative  of  the  drag  with  respect  to  V  (or  perhaps 
V^)  is  a  function  of  time  alone.  Using  velocity  as  an  exampie.  this  can  be  written  as. 


(149) 


where  Dy(T)  represents  the  value  of  the  drag  for  a  unit  step  fraction  of  velocity.  This 
equation  can  be  evaluated  by  taking  the  Laplace  transform  of  both  sides  to  give 

0(4*)  ~  (150) 

and  hence, 

0,U >  -  Mii-  <15D 

Equation  (151)  was  derived  under  the  hypothesis  that  a  Dy(t)  could  be  defined 
which  was  independent  of  the  driving  force  and,  hence,  of  the  velocity  history.  It  was 
not  necessarily  expected  that  this  would  be  true,  but  that  by  examining  the  test  data 
and  seeing  if  the  resultant  curves  would  collapse,  a  proper  functional  relationship 
could  be  obtained  which  would  permit  a  step  function  to  be  defined  which  could  be  used 
to  reproduce  the  test  time  histories.  As  an  example  of  this,  two  other  definitions 
which  were  considered  to  define  the  step  function  were 

0(t)  -  jtff  (152) 
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and 


/c 

/>(VcJ  -  Jjzf  D** ir) 


Equations  (149),  (152),  and  (153)  together  with  alternative  formulations  were 
examined  to  determine  which  would  best  collapse  the  test  data  for  die  different  weight 
increments.  Equation  (133)  showed  the  most  premising  results  of  all  those  tested. 
However,  the  preliminary  tesi  data  indicated  a  discrepancy  for  the  various  weight  dif¬ 
ferences.  While  the  shapes  of  the  resulting  curves  were  similar,  the  magnitude  dif¬ 
fered.  In  order  to  bette  ~  correlate  the  data,  a  step  function  was  added  to  the  sedation 
to  account  for  data  beyond  the  test  period.  The  magnitude  of  toe  step  was  set  equal  to 
the  last  data  point,  so  that 


/  /  -41 ftr  r 

/  Of*  f**t>J-*  =  D.\-±J 

c  L 

and,  assuming  s>o 

X(°<(vx)}  "  b*  ^ - 


(154) 


(155) 


Through  the  use  of  this  modified  expression,  the  parameter  Dy(s)  was  evaluated 
for  0<  8  <5.0  for  incremental  weights  of  3,  4,  5,  6,  and  7  pounds.  Figure  87.  Excel¬ 
lent  agreement  is  shown  for  the  3-,  4-,  and  7-pound  test  runs,  while  the  5-  and  6 -pound 
runs  did  not  conform.  A  detailed  examination  of  the  testing  technique  j  and  tie  instru¬ 
mentation  arrangement  revealed  that  the  instrumentation  was  simply  not  sensitive 
enough  for  runs  with  incremental  weights  greater  than  7  pounds. 

In  spite  of  the  inadequacy  of  the  test  data  accuracy,  an  attempt  was  made  to 
determine  the  time-dependent  response  to  a  step  input.  For  this  purpose  the  results 
obtained  by  the  definition  of  Equation  (153)  were  used  and  the  inverse  Laplace  trans¬ 
form  was  obtained  to  see  if  any  additional  problem  areas  would  be  encountered. 

A  study  of  the  properties  of  the  Laplace  transforms  ina.cated  that  the  values  of 
2 

the  quantity  D(s)/V  (s),  Figure  88,  would  have  to  approach  a  constant  for  large  values 
of  8.  This  was  required  to  retain  a  finite  force  at  the  start  of  acceleration.  For  this 
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FIGURE  87.  TRANSFORMED  EXPRESSION  FOR  DRAG  FORCE  DUE  TC 

IMPULSIVE  VELOCITY  SQUARED 


FIGURE  88.  TRANSFORM!  D  DATA  FOR  DRAG  DUE 
TO  IMPULSIVE  VELOCITY  SQUARED 


reua»,  the  raises  of  tfc*  ratic.  yje  evaluated  at  large  values  of  s  as  shewn  is  Figure 
89.  Oner  this  larger  of  s  the  previously  obtained  correlaliot  between  the  test 

loads  was  lost.  This  is  not  surprising  since  the  increase  of  the  parameters  corre¬ 
sponds  to  a  decrease  o  the  rt  parameter  so  that  the  scatter  obtained  at  large  s  values 
caa  be  attributed  to  errors  in  measuring  the  very  lor  forces  at  the  very  early  part  of 
the  model  acceleration. 

A  typical  curve 

P(*)  _  Ao  a.  ^  2  Z  <156) 

a.-*- »« 

was  selected  within  the  range  of  the  test  data  scatter  and  the  drag  variation  with  time 

_2 

coe  to  a  step  input  of  V  was  obtained  as  shewn  in  Figure  99. 

This  study  has  shown  that  a  possibility  exists  of  obtaining  the  time-dependent 
cn 1 3ii  of  forces  and  moments  of  a  shape  accelerating  from  rest  using  a  combination 
of  testing  in  a  relatively  simple  ureter  tank  and  analysis  A  these  results  by  using  the 
indicia!  function  concepts  of  Reference  45.  Although  the  present  study  was  restricted 
to  the  drag  of  relatively  simple  models,  no  special  difficulties  should  be  encountered 
in  treating  any  force  or  moment  history  for  much  more  complicated  shapes. 

Although  no  practical  and  usable  results  hav  ?  been  obtained  to  enable  the  predic¬ 
tion  of  aerodynamics  of  V/STOL  2ircrait,  the  procedure  described  appears  very 
promising  for  future  refinement  and  ultimate  use  in  obtaining  prediction  techniques. 
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APPENDIX  I 


EFFECT  OF  ROTARY  VARIABLES 

In  Figcre  1  p,  q  and  r  represent  the  magnitudes  of  the  rotary  variables  defined 
in  the  coordinate  system  shown-  The  effect  of  the  rotary  variables  p,  q  and  r  on  the 
development  of  the  jet  may  be  represented  by  a  change  in  the  mainstream  velocity 
vector.  The  perturbed  mainstream  velocity  vector  at  a  given  point  along  the  jet  center- 
line  is  expressed  as 


Op  =  (p?+  +  ri}  *  & 


/K  ♦  &  _ 

where  i,  j,  at  are  unit  vectors  in  the  fixed  coordinate  system  and  R  ts  the  position 
vector  of  the  given  point  along  the  jet  centerline. 

In  evaluating  the  jet -induced  flow  field,  the  computations  are  carried  out  for  a 
number  of  segments  as  shown  in  Figure  2.  Each  segment  is  treated  as  a  sepc  .»te  jet, 
with  proper  initial  conditions  and  the  appropriate  mainstream  velocity.  The  number 
of  integration  steps  per  segment  may  be  arbitrarUy  specified.  The  perturbed  main¬ 
stream  velocity  vector  is  assumed  constant  over  the  extent  of  each  segment. 

In  the  following  discussion,  subscripts  1  and  2  refer  to  jet  properties  of  the 
segments  I  and  II  of  Figure  2. 

The  ie*  of  diameter  do^  is  located  at  Xo^,  Yo^,  Zo^.  Initial  conditions  for  seg¬ 
ment  I  are 

-  1  ,  Uj  -  Un  =-  m,  —  Uici 

*'3jc i  U 


The  initial  jet  exhaust  vector  is  utilized  as  an  approximation  to  locate  the  midpoint  A 

lUpi! 

of  segment  I.  At  this  point  —  is  determined  bv  vector  addition  as  shown,  where 

L  oC 


i  •*." 


FIGURE  1-1.  DEFINITION  OF  ROTARY  VARIABLES 


FIGURE  1-2.  COMPUTATIONAL  SCHEME  FOR  ROTARY  VARIABLES 


V*  = 
R*  - 


do  (pi  +  + 

U«o 

]L 

do 


The  direction  cosines  of  the  perturbed  mainstream  velocity  vector  Upj  are 
determined  and  a  local  coordinate  system  is  established,  centered  at  Xo^,  YOj,  Zo^, 
and  aligned  with  the  perturbed  mainstream  velocity  vector  and  die  jet  exhaust  vector 
(see  discussion  of  arbitrarily  directed  jet). 

The  effective  inverse  velocity  ratio  for  segment  I  is 


-  U»t  _y^  =5  m,/ 1  Up.? 

U®  lUpil  •  Goo 


The  equations  for  Uj,  d,  and  X  are  then  integrated  nu  me  “■rally  over  the  extent 
of  segment  I. 

Point  Xo2,  Yo2,  Zo2  then  becomes  the  origin  of  the  next  jet,  with  a  diameter  of 
<k>2  =  d*  dOj,  where  d*  is  the  last  computed  value  of  the  aondii.'.ensionalixed  jei  diam¬ 
eter  of  segment  I. 

Other  initial  conditions  for  segment  II  are 


1. 


=  Jin-  =  i 


u. 


m,  * 


— 

Dee 


=  U,  on, 


where  13*  is  the  last  computed  value  of  the  nondimensionalizod  jet  velocity  in  segment 

I. 

The  jet  velocity  direction  after  the  integration  overjhe  extern'  of  segment  1  is 

!UiV! 


used  to  approximately  locate  midpoint  B.  At  this  point 


U  oo 


-  is  evaluated  and  die 


direction  cosines  of  the  perturbed  mainstream  velocity  vector  Up2  are  determined. 
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A  new  local,  jet-oriented  coordinate  system  is  established,  centored  at  Xo2, 

Yo2,  Zo2.  The  effective  inverse  velocity  ratio  for  this  segment  is 

The  equations  for  U£,  d4,  X*  are  then  integrated  over  the  extent  of  segment  II. 

The  computations  described  above  continue  until  integration  over  the  extent  of 
the  entire  jet  has  been  accomplished.  Although  each  segment  is  treated  as  a  separate 
jet,  the  segments  are  linked  to  each  other  in  determining  the  degree  to  which  the  Jet 
has  deformed  from  its  initial  circular  cross  section.  This  is  done  by  setting 

m  •&!  ,  *&c>j  m  Z&z 

l.e.,  specifying  the  Initial  ratio  of  minor  to  major  axis  for  a  given  segment  to  bo  oqual 
to  the  last  computed  value  for  the  previous  segment. 

The  Induced  velocity  components  at  a  given  control  point  aro  determined  by  sum¬ 
ming  the  velocity  components  due  to  each  segment  of  the  jet.  Again  each  segment  is 
treated  as  a  separate  jet  for  these  computations. 

Figure  3  shows  the  effect  of  tho  rotary  variablo  r  on  the  conterlino  of  a  Jet 
exhausting  normally  into  a  crossflow  at  a  velocity  ratio  U^/ U^0  «  .  125.  One  curve 

(r*  «  -  0. 0)  shows  the  centerline  of  the  undisturbed  Jet  exhausting  into  a  crossflow 

.it  this  voloclty  ratio.  It  was  computod  by  considering  tho  jet  to  consist  of  a  number  of 
sogmonts,  as  described  in  dot&ll  In  the  preceding  discussion,  with  V*  •  0.0  for  oach 
segment.  Tho  contcrline  computed  4n  this  manner  deviates  less  than  1  percent  from 
previously  computod  results. 

Tho  sign  convention  for  tho  rotary  variable  r  was  established  in  Figure  1.  For 
r*  <  0.0  the  perturbed  mainstream  velocity  |Up|<U,o  and  the  Jet  centerline  shown  less 
defloctlon.  For  r*  >0.0  the  perturbed  mainstream  velocity  |U^J  >UW  which  results 
in  a  greater  deflection  of  the  jet  centerline. 

Induced  velocity  components  in  the  plane  of  the  jet  exit  were  evaluated  for  two 
values  of  r*.  Table  1  shows  the  incremental  change  in  volocity  components  at  three 
control  points  as  a  function  of  the  rotary  variable  r*.  It  can  be  observed  that  A  (u/U)/d  r*. 
4(v/U)/dr*  and  A{w/U)/Ar*  arc  not  constant  at  each  control  point.  Since  the  incremen¬ 
tal  changes  in  induced  velocities  are  directly  rotated  to  the  induced  forces  and  moments 
due  to  the  rotary  variable,  it  is  desirablo  to  establish  the  range  over  which  they  could 
be  considered  linear  with  r*. 
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Distance  in  Freestream  Direction,  X/d0 


FIGURE  1-3.  EFFECT  OF  ROTARY  VARIABLE  ON  A  JET  CENTERLINE 


^r^F-fT'  *"v 
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TABLE  I-I.  EFFECT  OF  ROTARY  VARIABLE  ON  INDUCED  VELOCITY  COMPONENTS 


q  =  0. 
q  =  .l 
C/C 


q  =  0. 
q  =  .l 
C/C 


X 

=  .95,  Z  = 

.7 

X 

=  .55,  Z  = 

7 

X 

=.15,  Z 

=  .7 

u/U 

v/U 

w/U 

u/U 

v/U 

w/U 

u/U 

v/U 

w/U 

.0180 

-.0307  -. 

0411 

.0337 

-.0162  - 

rH9 

.0273 

-.0016 

-.0095 

.0179 

-.0303  -. 

0407 

.0334 

-.0159  - 

.0269 

-.0015 

-.  0093 

.0177 

-.0299  - 

0410 

.0332 

-.0156  - 

0294 

.0267 

-.0013 

-.0094 

.95,  Z 
v/U 
.  1073 
.1070 
,  1075 


=  .l 
w/U 
-.0239 
-.0236 
-.0242 


X  =  .55,  Z  =.  1 
u/U  v/U  w/U 
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The  velocity  components  of  Table  1  were  obtained  by  considering  an  equal  number 

A 

of  integration  steps  along  the  Zf  -axis  in  solving  for  the  geometry  of  the  jet.  This 
means  that  for  every  value  of  r*  a  different  length  s  of  the  jet  was  considered  in  eval¬ 
uating  the  induced  velocity  components.  In  an  effort  to  separate  the  effects  due  to  a 
change  in  the  rotary  variable  from  the  effects  due  to  integrating  over  a  variable  length 
of  the  jet,  two  different  approaches  to  establishing  an  equivalent  length  s  were  utilized. 

One  approach  was  to  sum  contributions  to  each  induced  velocity  component  due 
to  the  singularity  distribution  along  the  jet  centerline  until  the  contribution  from  a  given 
segment  represented  less  than  .  1  percent  of  the  total  up  to  this  point.  Further  contri¬ 
butions  were  neglected.  The  other  method  utilized  a  straight  line  approximation  for 
the  jet  centerline  to  determine  an  equal  length  s  of  the  jet  for  each  value  of  r*. 

For  both  of  the  methods  described  above,  the  axis  of  rotation  was  taken  through 
the  jet  exit  and  the  behavior  of  the  incremental  changes  in  velocity  components  at  the 
various  control  points  was  examined.  For  the  range  of  r*  considered,  which  encom¬ 
passes  a  variation  of  r  from  .1  to  .4  radians/sec,  it  appears  that  the  incremental 
changes  in  velocity  components  are  small  enough,  relative  to  the  overall  accuracy  of 
the  method,  to  preclude  a  meaningful  determination  of  the  quantities  d(u/U>^r*, 
4(v/U)^r*,  4(w/U)ftr*. 

Computations  to  determine  the  effect  of  rotation  about  an  axis  other  than  the  jet 
exit  were  carried  out.  Rotation  about  the  center  of  gravity  in  a  reaction  control  con¬ 
figuration  for  the  XV-6  was  chosen  as  an  example.  Computations  with  the  moment  arm 
/=  11.2  ft,  dQ  =  .2  ft,  u.Q  =  400  ft/sec,  U»  =  50  ft/sec  and  q  =  .  1  rad/sec  again 
showed  small  changes  in  the  induced  velocities  evaluated  in  the  plane  of  the  jet  exit. 

A  check  case  in  which  the  freestream  velocity  was  assumed  to  be  uniformly  perturbed 

over  the  extent  of  the  jet,  U  =  dU,  where  dU  =  lx  q  =  1. 12  ft/sec,  was  also  com- 

P 

puted.  Table  3  shows  induced  velocity  components  at  various  control  points  for  q  =  0. , 
q  =  .  1  and  the  check  case.  Again,  indications  are  that  the  effects  due  to  the  rotary 
variables  in  a  reaction  control  configuration  are  negligible. 
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APPENDIX  D 


DETERMINATION  OF  TRAILING  VORTEX  INTEGRAL 


The  integral 


Urn- 

f  A 

JL  r* 


in  Equation  (121)  may  be 


written  in  the  form 


xf* 


Wktrc  X  =  ft  4-4  X,.  CX, 


using  the  same  notation  as  Equation  (123) 


An,  ,  *•*] 


(hJc 

(Wfr-i  ) 


ft,  ^  -tSSiA 


4  «  4n*i.  u<_ 


S  ° 


.2. 

«.  «  x*  +  (j-  t  -  a,  ])  1-  +*,  ^  ) 

i  =  -a  £*-•»,  Ly  i-o.i)  +(a  <-*,  j)i,  ] 


't-  I  w 

C  =  li-4,  +i, 
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Thu*  •  £»**.v 3 f  *-£.‘*<.7  -  *[t  t)‘. '*>.* 'V.iJ ; 

♦  *1>.  v«.  -V*  -C«J  f 

♦  HV  <-*/<*' +1?  +2*,*}  -2l,xi+24, ij 

Now  l  +  [l-tJ)  Ijj 

So  that  with 

A,  *  /+C 

A*.  -  «.*>0 
A3  r  J*, 

**/♦*■<  *««*.  *1 
As  *  -*£«,$  -i.*  +V* 

A4  =  J  •#*•#•*,  x.  +  *.  X  J  ♦<*,  £  +  .2<»,XJ  -2b,Xi: 

+  .a*,  *> 

We  obtain  AC  -  4  *  A(  £  *"  J  *’’  ♦  Aj  *  A*.  ?  f /fr  J  +■ 

and  then  cue- ty  -  A]'’  +1]  +  C 

where  A  =  *,  +  A*  fa^4  +  *Z  lH  U*4> 

*■  > 

B  *  ^  ih  -?;  /*  W.]  '1'  A*~<  +  *,[/.  -j>  /X  w.] 

■+■  A(.  +  A5  Ijff  fa*  4 
f 

C  -  4*  [/.  >K  +  4f  [  l  -Jr  /j/  U*f]i*c 

Now  Jcj+i 

where  p  -  [>/  -A, 

~  /-  hi  -x  -h*.  m  ~i,  3- 

*  s  * 
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Aftfy+C 


Fj+Gr 


Ho*  L*-0  +$-'iT+^  ~  ?  +E 
*»  B  =  -.?  r# » * » -(/.  -  7.  'i'w.)  'f 

and  *V  -J*fj[  -7;%1~t-]*-[t  -J.  % 

C  «  - - - 2 - * - 

!♦  A 


Thus  Equation  (121)  becomes 


*  ■*■  j  ^  - 7  +•  ^A  C*~1)l  f  p-^rr: — rr  ^«04,"1 

*  4r  U  7  ^r;4  37J  l  JZ&jb+he^  J 


Ay"’*??  iC 

Now  j-;  *  4t^.(l*.  (*'j)  -  +  X 


£*»*•*. 

with  /f  ^  I  +  4,  Uu^  | 

and  I  =}+*•*  -*i  [I  ~Z  Ijtf  ] 


so  that  finally  we  obtain 
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